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Versuche an Turbinenlaufschaufelgittern. 


Von H. Hausenblas. 


1. Einleitung. Uberblickt man die bisher in allgemein zugiinglicher Weise verdffentlichten 
MeBergebnisse iiber das Verhalten von Profilgittern!~8, so stellt man mit Ausnahme der neu- 
esten englischen Literatur ein fast vollstandiges Fehlen moderner Unterlagen fiir stark um- 
lenkende Beschleunigungsgitter fest, wie sie im Dampf- und Gasturbinenbau Verwendung fin- 
den. Soweit es sich um Waagenmessungen handelt, ist auBerdem einzuwenden, daB diese immer 
einen Integralwert der Verluste iiber der gesamten Schaufelhéhe einschlieBlich der héheren Rand- 
verluste® '° angeben, und eine Abtrennung der letzteren nicht méglich ist. Die Druckvertei- 
lungsmessungen in den Schaufelkandlen hingegen 
liefern keine Angaben iiber die aufgetretenen 
Verluste. Dasselbe gilt von den Interferometer- 
messungen E. Eckerts. Die Impulsschnittmethode 
hingegen gestattet die Vermessung der reinen 
Gitterverluste unter Ausschaltung der Rand- 
verluste, da man nur in einem einzigen mittleren 
Schnitt entsprechend langer Schaufeln miBt. 

Im folgenden soll eine Reihe von Versuchs- 
ergebnissen fiir stark umlenkende Beschleuni- 
sungsgitter nach der Impulsschnittmethode mit- 
eeteilt werden, die aus den Jahren 1942 bis 1944 
stammen'!, Obwohl es sich dabei nicht um 
Versuchsergebnisse aus systematischen Unter- 
suchungsreihen handelt, diirften diese doch von 
einigem Interesse sein. 

2. Die Versuchseinrichtung. Abb. 1 zeigt den 
Priifstand und zwar einen Blick auf die Ab- 
strémseite des Gitters mit der Durchkurbelein- 
richtung, die das Verfahren eines Prandtlrohres 
bzw. einer (in der Abb. nicht gezeigten) Zwei- 
fingersonde parallel zum Gitter gestattet. Damit 
konnte die Verteilung der Abstrémgeschwindig- 
keit hinter dem Gitter nach GréBe und Richtung 7% 
aufgenommen werden. Es wurde im Uberdruck- Abb. 1. Profilgitter-Versuchsstand 
betrieb gearbeitet, da bei Saugbetrieb und der 
iblichen Art der Versuchseinrichtung u. U. eine Fiihrung des aus dem Gitter austretenden Strah- 
les hatte auftreten kénnen, was durch Uberbestimmung der Umlenkung im Gitter zu einer Ver- 


1 4. Stodola, Dampf- und Gasturbinen. Berlin 1924. 

2 E. Knérnschild, Untersuchungen an Turbinenschaufelgittern mit Hilfe einer Schaufelwaage. Jb. 1941 
Itsch. Luftf.-Forschg., S. II 195/210. 

3 E. Knérnschild und K. Leist. Untersuchungen an Turbinenschaufelgittern, Jb. 1939 dtsch. Luftf.- 
Forschung., S. II 204/212. ; ra 

4 W. Traupel, Neue allgemeine Theorie der mehrstufigen axialen Turbomaschinen, Ziirich 1942. 

5 K. Christiani, Experimentelle Untersuchung eines Tragfliigelprofils bei Gitteranordnung. Luftf.- 
Forschg. 2 (1928), Heft 4. 

6 4A, R. Howell, J. Roy. Aeronaut. Soc. 52 (1948), 5.329. | 

7 K. W. Todd, Proc. Instn. Mechan. Engr. 157 (1947), 5. 482. 

8 D. G. Ainley, Proc. Instn. Mechan. Engr. 159 (1948), 5. 230. 

9 M. Sédille, Bull. techn, de l’Assoc. des Ingénieurs sorties de l'Université de Bruxelles, Nr. 3 (1948), 
>, le 

10 4, D. S. Carter, Proc. Instn. Mechan. Engr. 159 (1948), S. 255. aes 

11 An der Durchfiihrung und Auswertung der Versuche waren die Herren G. Rieger, S. Mielke, 
Vi. Lamblin und A. Pfleghaar beteiligt. 
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failschung der Gitterversuche gefithrt hatte. Ein Kinfluf des Abstandes der MeBebene hinter dem | 
Gitter auf die MeBergebnisse konnte fiir nicht zu kleine Abstinde nicht festgestellt werden. Die 
in Abb. 1 eingebaut gezeigten mittelsten Schaufeln des aus Holz gefertigten Gitters sind in Mes-_ 


sing hergestellte Sonderschaufeln, die im Mittelschnitt tiber den Umfang verteilte DruckmeB- 


stellen besitzen. Die verschiedenen Zustrémrichtungen wurden durch Auswechseln der Zustrém- | 
kisten erzeugt. Die untersuchten Gitterformen sind aus Abb. 2 und deren Hauptdaten aus Ta- | 


GitterB7 Gitter B2a 


Gitter B2b Gifter C7 Gitter C2 


Abb. 2. Untersuchte Profilgitter. 


Tabelle 1. Hauptdaten der untersuchten Gitter. 


Ronste. oncen : : Relative Relative | Relativer 
Profilgitter | Eintritts- iAciatrittse Staftelungs- Teilungs- Relative Kriim- Relative Dicken- | Eimtritts-] Strich- 
Sopa emicel winkel verhaltnis Pfeilhéhe mungs- Dicke micklace kanten- art 
riicklage : be radius 
= Bi Bs y 1/1 fil x4/l d/l xq/l ofl = 


54,5° 34,2° OU” 0,705 0,249 0,39 0,282 0,32 0,080 | ——— 

3020 33,0° 80,2° 0,733 0,395 0,40 0,233 0,35 0,033 }|——-—= 
1 39,0° 30,0° T° 0,790 0,336 0,42 0,248 0,26 0,063 aac uee 
2a 48 ,5° 32,0° HOS” 0,692 0,305 0,38 0,252 0,23 DADE: |) = o0— 
2b D0zom 30,0° 68,5° 0,800 0,305 0,38 0,252 0,23 050763) ee 
1 1350. 32,0- SA) 0,770 0,208 0,30 0,274 0,26 0,074, === 
2 42,5° 36,4° (Osoe 0,572 0,306 0,41 0,313 0,35 0.0765 | 


OOWWAeO 


belle 1 zu entnehmen, wo auch die fiir die einzelnen Gitter in den spateren Abbildungen an- 
gewendeten Stricharten angegeben sind. Die benutzten Bezeichnungen sind aus Abb. 3 zu er- 
sehen. Die Reynolds- und Machzahl, bezogen auf die Zustrémgeschwindigkeit zum Gitter und 
die Sehnenlinge der Schaufeln, waren bei den Versuchen im Mittel Re = 200000 und M = 0,128. 
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Die Wahl der Reynoldszahl der Versuche in der Nahe des kritischen Bereiches derselben war durch 
die Verhaltnisse in der betreffenden Turbine bedingt. 


3. Die Versuchsergebnisse. Die in den Gitterversuchen erhaltenen Auftriebsbeiwerte c, und 
Gleitzahlen ¢ sind in Abb. 4 wie bei Howell! in Abhangigkeit vom Winkel i zwischen Zustrém- 
richtung und Eintrittstangentenrichtung der Profilskelette zusammengestellt. Dabei wurde nicht 
c, selbst, sondern ¢, - I/t aufgetragen, da diese Wertekombination mit dem Umlenkungsdreieck des 
Gitters in direktem Zusammenhang steht. Ferner wurden an Stelle der Widerstandsbeiwerte 
die Gleitzahlen verwendet, da letztere in den Wirkungsgradformeln auftreten2,*. An Hand des 
Verlaufs der Gleitzahlen in Abb. 4 kann man deutlich folgende auch an den Profilumrissen in 
Abb.2 unterscheidbaren Profilfamilien fest- 
stellen: 


4. Die Profilfamilien. Die Profile der 1, 
ersten Familie (Profilgitter 0, C1, C2) sind ae 
nach der im alteren Dampfturbinenbau ein- ie 


gebiirgerten Methode aus méglichst wenigen a "8, 
v4 


Kreisbiégen zusammengesetzt und zeigen eintritstangente . 
dementsprechend — wie noch spater be- 
sprochen — nicht nur in der Nahe der Ein- 
trittskante, sondern auch infolge des Kriim- 
mungssprunges auf ihrem Riicken* starke 
Ubergeschwindigkeitsstellen. Sie zeigen da- 
her ihre besten Gleitzahlen bei Zustrém- 
winkeln 6,, die den Skelettlinien-Eintritts- 
winkel f{ wesentlich iiberschreiten (also bei 
negativem 1), wo laut Abb. 5 (in der das 
Verhaltnis der Abstrémgeschwindigkeit w, 
zur Zustrémgeschwindigkeit w, dargestellt 
ist) die Strémung im Gitter insgesamt derart 
beschleunigt wird, da8 die Ubergeschwindig- ee Nee uae b 
keitsstellenrelativ dazu geniigend verschwin- astride tandente 

den. Die Gitter dieser Familie sind also trotz 
der relativ starken Kriimmung ihrer Profile 
nur bei kleinen Umlenkungen giinstig ver- 
wendbar. Das Profilgitter C2 zeigt infolge 
seiner relativ zur Teilung zu groBen Profil- 
dicke und damit schlechten Kanalform 
héhere Gleitzahlen als die iibrigen. 

Die zweite Familie, der die Profilgitter 
Bl, B2a und B2b angehiéren, zeigen ver- 
feinerte Profilformen, aber ahnlich dicke 
Eintrittskanten wie Familie eins. Die Folge Abb. 3. Bezeichnungen im Profilgitter. 
dieser Verfeinerung ist eine Verminderung 
des Bestwertes der Gleitzahl und eine Verstaérkung der Empfindlichkeit der Profile gegeniiber 
Anderungen der Zustrémrichtung. Dabei sind die beiden mit relativ weiter Teilung versehenen 
Gitter B1 und B 2b dem enger geteilten Gitter B2a in beiden Hinsichten deutlich iiberlegen. 
Die von Friedrich® in demselben Bereich der relativen Teilungen festgestellte umgekehrte Ten- 
denz des Einflusses einer Teilungsinderung ist offensichtlich auf die mit scharfen Kanten be- 
haftete UmrifBform seiner Profile zuriickzufiihren. Diese Kanten kénnen nur bei durch enge 
Teilung bedingter straffer Strémungsfiihrung ohne starkere Ablésung umstrémt werden. Die 
mit steigender Teilung eintretende Verbesserung der Profilgitter wurde in dem hier betrach- 
teten Bereich der relativen Teilungen auch von Talbert und Smith® erhalten. 


Skelettlinie’ 


Umlenkungsdre/eck : 


1 Siehe FuBnote 6 von S. 75. 

2 C. Keller, Axialgeblase vom Standpunkt der Tragfliigeltheorie. Zirich 1934. 

3 W. Traupel, a.a.O. 

4 A. Betz, Luftf.-Forschg. 19 (1942), S. 129. 

5 H. Friedrich, Strémungsmechanik einer Dampfturbinenstufe bei wechselnden Betriebszustanden. 
Nurnberg 1948. 

6 J. E. Talbert und J. C. Smith, J. Aeron. Sciences 15 (1948), S. 556. 


6* 
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Die dritte Familie wird nur durch Gitter A vertreten. Die Profile desselben enthalten ahn- | 
liche Verfeinerungen der UmriSform wie jene der zweiten Familie, unterscheiden sich aber von _ 
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Abb. 4, Auftriebsbeiwerte und Gleitzahlen, Abb. 6. Verlustschleifen fiir den Gesamtdruck. 
letzteren durch eine relativ schlanke und scharf vorgezogene Eintrittskante. Man erkennt aus 
Abb. 4, daB diese scharfe Eintrittskante durchaus nicht so vorteilhaft ist, wie man zunichst 
glauben kénnte, indem die minimale Gleitzahl héher hegt als bei allen Gittern der zweiten Fa- 
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milie. Da auBerdem die Empfindlichkeit gegen Zustrémwinkelénderungen wesentlich gréBer 

ist als bei der zweiten Familie, wurde diese Profilfamilie nicht weiter untersucht. 
Bemerkenswert sind noch die gegeniiber Familie eins hohen Umlenkungen bzw. Auftriebs- 

beiwerte, die die Gitter der zweiten Familie bei besseren Gleitzahlen erreichen, was bedeutet, 


dafs man mit Familie zwei wesentlich gréfere Stufengefalle mit guten Wirkungsgraden verar- 
beiten kann. 


5. Die Gesamtdruckverteilung. Interessant war ferner der Verlauf der hinter dem Gitter 
vermessenen Gesamtdruckverteilungen, von denen drei charakteristische in Abb. 6 skizziert sind, 
Verteilung A entspricht dem Normalfall ohne stirkere Ablésungen, wahrend B und C Vertei- 
lungen mit starkeren Ablésungen an Schaufelbrust bzw. Schaufelriicken bei stark von der stof- 

ai freien KEintrittsrichtung abweichenden Zu- 

ja Be strémungen sind. 


| 


8 a a ; 


34° 


2 
7 
Ole a 4 by 
=r 27" =0 #9 =. 20° 30° Ga” a ae | 
Abb. 8. Mittlere Austrittswinkel der Strémung 
Abb. 7. Auftriebsbeiwerte. aus dem Gitter. 


6. Die Auftriebsbeiwerte. Um nun den Verlauf der erhaltenen Auftriebsbeiwerte naher analy- 
sieren zu kénnen, sind diese in Abb. 7 in der in der Gittertheorie gewohnten Weise in Abhingig- 
keit vom mittleren Anstellwinkel «, =f, — y aufgetragen und mit Hilfe der bekannten Be- 
ziehung 1~4 


Ca = a, Ky sin (4 ,, + 4) 


Tabelle 2. Umlenkverhalten der untersuchten Gitter. 


Aus den Auftriebsbeiwerten Aus den Austritts- 


nach a, K, und a, mit der Methode _ Profilgitter wohl Sie Ee inpetuecereayet 
der kleinsten Quadrate ausgewertet , a) Ky Oo Bo Boo CH 

(vgl. Tabelle 2), wobei die spater einer O 5,01 36,9° 30,6° Boat 0,9725 
Betrachtung unterzogenen stark ab- A 5,73 46,4° 33,8° 35,1° 1,0606 
weichenden Punkte bei a, > 30° ae ae Te ea eS eae 
(6, < 30°) nicht beriicksichtigt wur- B ob 5,93 AL1° 27.4° 95,2° 0.8679 
den. Dabei ist ao eine den Anstieg Cl 5,85 30,4° 26.6° 27.62 0.9159 
des Auftriebsbeiwertes angebende One 4,18 570° 153° DD BX° 0.8149 


1 F, Weinig, Die Strémung um die Schaufeln von Turbomaschinen. Leipzig 1935. 
2 G. Klingemann, Ing.-Arch. 11 (1940), 5. 151. 

2A Ce deh Wislie eres Trans. Amer, Soc. mechan. Engr. Aug. (1945), 5. 451. 

4G, F. Wislicenus, Fluid Mechanics in Turbomachinery. New York 1947. 
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Konstante, K, der GittereinfluB-Korrekturbeiwert, o der die Nullauftriebsrichtung des Profils 
im Gitterverband festlegende Winkel. Der Gitterkorrekturbeiwert Ko ist vom Teilungsverhalt- 
nis (t/l),, und Staffelungswinkel y,, des aquivalenten Streckenprofilgitters abhingig, welche 
Werte (insbesondere der Staffelungswinkel) von den entsprechenden Werten des wirklichen 
Gitters abweichen und am bequemsten durch konforme Abbildung mit Hilfe der elektrischen 
Analogie!~4 erhalten werden kénnen. Da bei den hier beschriebenen Versuchen (t/1),, und 


y nicht ermittelt wurden, ist die Trennung der beiden Faktoren von dg Ky nicht méglich. ve 
die in den Versuchen erhaltene Nullauftriebsrichtung mit jener fiir Potentialstromung identisch 
sein, SO ware Y= Y — %- \- 7 
7. Die Austrittswinkel. Trotzdem durch die 30%-“* l r r 4 
eben durchgefithrte Betrachtung der Auftriebs- 
beiwerte das Umlenkverhalten der Gitter bereits 
charakterisiert ist, wollen wir in Abb. 8 noch- 
mals die gemessenen Austrittswinkel der Gitter 25% 
direkt betrachten, da bei kompressiblem Stré- 
mungsmittel das Arbeiten mit dem Auftriebs- 
beiwert schwierig wird. Dieselben Uberlegungen, 
die zu dem oben erwahnten Gitterkorrekturbei- ig J | 
wert fiir den Auftriebsbeiwert fiihrten, ergeben ee / | 
fir den Abstrémwinkel f, vom Gitter bei j 
Potentialstrémung die Beziehung j/ 
ctg Bo — Cy ctg Boo — (1 > Cy) ctg Px ? 15% a 
wobei der Gitterkorrekturbeiwert Cy in Ab- /j | 
hangigkeit vom Teilungsverhaltnis und Staffe- /; 
lungswinkel des dquivalenten Streckenprofil- ji / 
gitters bekannt ist *+,° und immer unter eins liegt. at 7 Le 
Ve 7 
G5% =| 
6, Coat 
68 °o fo oO o ° oO ° ° tt 0 | ___ a 
Le AO EE be 0” a 59 2 3 4 5 6 7 8 9 
Abb. 9, Kontraktionsziffer der Strémung im Gitter. ; Abb. 10. Interpolation der Gleitzahlen. 


In Abb. 8 sind nun neben den wegen der Streuung durch geknickte Kurvenziige verbundenen MeB-! 
punkten, die mit Hilfe der Methode der kleinsten Quadrate erhaltenen Ausgleichskurven eingetragen. 
und in Tabelle2 die entsprechenden Werte von Cy und des Abstrémwinkels /.) bei ablenkungs-: 
freier Durchstrémung des Gitters zusammengestellt. Bei dieser Auswertung wurden die stark ab-: 
weichenden MeSipunkte fiir 6,;< 30° aus spiater besprochenen Griinden herausgelassen. Auer- 
dem ist zu beachten, daB die Gitter trotz ihrer groBen Schaufellangen nicht in ebener Strémung. 
arbeiteten (vgl. Kontraktionsziffer u = ¢,,,/¢m; in Abb. 9), was selbstverstandlich die Strémungs- 
austrittswinkel der Gitter und die GréBe der Werte Cy ebenfalls beeinfluBt. Dabei war inter- 
essant festzustellen, da® die Kontraktionsziffern nicht langs der Schaufellange konstant oder von 
Schaufelmitte nach den Schaufelenden zu einsinnig verlaufend waren. Die Gitter A und B1| 
haben ein dem theoretischen zuwiderlaufendes Verhalten und damit Werte Cy >1, doch werden 
auch die diesbeziiglichen MeBwerte durch obige Formel gut interpoliert. Dieses abnormale Ver- 
halten der Abstrémwinkel und die starke Verainderung der Abstrémwinkel fiir 6, < 30° ist dar-| 
auf zuriickzufiihren, da8 die Abstrémwinkel eines Gitters aufer vom potentialtheoretisch erfaB- 


1 P. de Haller, Application de l’analogie électrique a l’étude des grilles d’aubes. R | 
Cen aie q s grilles d’aubes. Rev. Techn. Sulzer 


K. H. Grofmann, Schweiz. Bauztg. 66 (1938), S. 429. 
L. Malavard, R. Siestrunck und P. Germain, Rech. Aéronautique, (1948), S. 33. 


W. Traupel, Zur Potentialtheorie des Schaufelgitters, Techn. Rundschau Sul | : 
5 Siehe FuBnoten 1, 3 u. 4 von S. 79. . Te LE GE e Ee came NA | 


mB Ww PD 
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baren Einflu® und den Abweichungen vom ebenen Strémungszustand noch von den im Gitter 
auftretenden Verlusten abhangen. Die bei wachsenden Verlusten sich verstarkenden Grenz- 
schichten und Abliésungsgebiete ergeben entsprechend wie die endliche Dicke der Profilaustritts- 
kanten eine Verminderung des Abstrémwinkels. Daher tritt ein dem theoretischen ahnliches Ver- 
halten der Abstrémwinkel und eine Veranderung der Auftriebsbeiwerte mit dem Sinus des Anstell- 
-winkels nur im Gebiet der besten Gleitzahlen auf, weshalb — wie schon erwihnt — bei den 
Auswertungen jene Versuchspunkte, bei denen die Gitter infolge Arbeitens als Verzégerungs- 
gitter (8, < 30°, 1. > 30°). hohe Verluste aufweisen, nicht beriicksichtigt wurden. Es sei noch 
erwahnt, daB die durch hohe Verluste bedingten Verringerungen der Abstrémwinkel die nach 
der Abschatzung von Traupel! berechneten Werte iibersteigen. Der merkwiirdige Verlauf der 
Auftriebsbeiwerte des Gitters C 2 im Bereich a. ~— 20 bis — 30° diirfte auf MeBungenauig- 
keit zuriickzufiihren sein. 


8. Die Gleitzahlen. Bei der Analyse der 


x Gleitzahlen ergab sich, da®B es nach Abb. 10 
10 in offenbar méglich ist, diese durch Parabeln fiir 
| \ den in der genannten Abbildung dargestellten 
49 r ; tT Wert e/(c,:l/t) von einer je nach Profilart 
j | verschiedenen Ordnung n entsprechend 
08 Ay ! \ E€ oN l l n 
iN ley ae erat 
; a fh \ = (Ca Ca — / opt 
a7 if ‘\ + \ t t / opt 
i? * hy 
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Abb. 11. StoBfaktor. Abb, 12. EinfluB der Rauhigkeit der Profile 
und des Turbulenzfaktors der Luft. 


anzunaihern, wobei C und n Konstanten sind und der Index ,,opt** die Werte fiir den Punkt mit 
dem niedrigsten ¢/(c, + 1/t) kennzeichnet. 


9. Die Geschwindigkeitsbeiwerte. Ausgehend von den im Dampfturbinenbau iiblichen Ge- 
dankengangen kann man aber auch den Verlauf der Kurve der Geschwindigkeitsbeiwerte iiber 
dem Zustrémwinkel /, als Superposition eines Grundverlustes y,,, (bei 4: ) und eines Kintritts- 
stoBverlustes zusammengesetzt denken, was fiir inkompressibles Medium auf die Beziehung: 


Yope [1+ x (ctg By ope — etg Py)? sin? B.]*/? 


fiihrt, wobei der StoBfaktor ~ ein Abminderungsfaktor gegeniiber dem Carnotschen StoBverlust 
ist. Die Auswertung einiger Gitter (Abb. 11) zeigte, daB® dieser StoBfaktor so stark variiert, daf 
er nicht als Konstante angenadhert werden kann. Das Maximum von liegt bei den Familien 
zwei und drei in der Nahe des Zustrémwinkels fiir Gleichdruckbetrieb, was sich auch bei Aus- 
wertung fremder Gittermessungen” mit wesentlich gréBerem 3, wiederholte. 


10. Einflu8 der Turbulenz. Da die vorstehend beschriebenen Gitterversuche knapp oberhalb 
jener Reynoldszahlen ausgefiihrt wurden, die dem kritischen Bereich fiir Profile entsprechen, 


1 Siehe FuGBnote 4 von §.75. 
2 Siehe FuBnote 3 von S. 75. 


Hausenblas: Versuche an Turbinenlaufschaufelgittern, Ingenieur-Archiy) 


wurde schlieBlich noch in einem Beispiel (Gitter 0) der Einflu8 einer erhéhten Turbulenz def 
Arbeitsluft und einer erhéhten Oberflaichenrauhigkeit der Profile festgestellt. Die dabei ethai 
tenen Geschwindigkeitsbeiwerte y des Gitters sind in Abb. 12 zusammengestellt. Beide Ein- 
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Abb. 13. Verteilung des Druckes um den Profilumfang. 


] 


fliisse ergaben merkliche Anderungen im cal 
warteten Sinne. 


11. Die Druckverteilung. Als letztes ist in 
Abb. 13 das Ergebnis von Druckverteilungs- 
messungen um den Profilumfang fiir Gitter B2a_ 
wiedergegeben. Aufgetragen ist der .,ideale“* 
Staudruck q’ (= Differenz von Gesamtdruck vor | 
Gitter und statischem Druck an der betreffen- 
den Stelle) dividiert durch den entsprechenden | 
Wert fir den Strémungszustand am Gitteraus- | 
tritt (q). Man erkennt au®er der Ubergeschwin-. 
digkeitsstelle in der Nahe der Eintrittskante | 
noch zwei weitere schwache Ubergeschwindig- | 
keitsstellen am Riicken entsprechend den Stellen | 
der Kriimmungsspriinge der Profilumri®Bform. | 
Bei entsprechenden Messungen am Gitter 0 war 
am Riicken in Ubereinstimmung mit dem dort 
vorhandenen staérkeren Kriimmungssprung eine 


wesentlich starkere Ubergeschwindigkeitsstelle | 


vorhanden, die zu einer entsprechend starken Ab- 
lésung und hohen Verlusten fiihrte. Zusammen- 


setzung von Profilformen aus einzelnen Kreis- 
bégen erscheint also trotz der Verwaschung der | 
Kriimmungsspriinge bei der Fertigung und der 
ausgleichenden Einfliisse der Grenzschicht als zu 
grob, weshalb vom Verfasser eine Profilsyste- 
matik mit kontinuierlichen Kriimmungsiiber- 


gangen ausgearbeitet wurde, die sich gut bewahrt. 

Ferner wird in Abb. 13 infolge der Bezugnahme auf q, verdeutlicht, in wie starkem MaBe die 
Strémung in der Nahe der Austrittskanten der Profile von der Zustrémrichtung f, unbeeinfluBbt 
bleibt. 


(Eingegangen am 28. Juli 1950.) 


Anschrift des Verfassers: Dr. H. Hausenblas, Saarbriicken 1, Francoisstr. 30. 
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Systematische Druckverteilungsmessungen an Pfeilfliigeln konstanter Tiefe 
bei symmetrischer und unsymmetrischer Anstrémung. 


Il. Mitteilune*. 
Von W. Jacobs. 


1. Einleitung. In der I. Mitteilung ist iiber umfangreiche systematische Sechskomponenten- 
messungen an Pfeilfliigeln verschiedener Pfeilung und Zuspitzung berichtet worden. Fiir diese 
Systematik war die wichtige Klasse der Trapezfliigel herausgegriffen worden. Bei festem Sei- 
tenverhaltnis (A = 5) wurden die Pfeilung (~ = —30°, 0, +15, 30, 45°) und die Zuspitzung 
(Z = 1,0; 0,6; 0,2) variiert. Besonderer Wert wurde hierbei auf die Bestimmung der Stabili- 
tatsbeiwerte gelegt. Die Ergebnisse stimmten z. T. recht gut mit theoretischen Werten iiberein. 
Um genauen Aufschlu8 iiber die Krafteverteilung im einzelnen zu gewinnen, wurden im An- 
schlu8 daran Druckverteilungsmessungen durchgefiihrt, und zwar bei symmetrischer und un- 
symmetrischer Anstrémung. Aus der bei den Sechskomponentenmessungen durchgemessenen 
umfangreichen Systematik wurden hierfiir zunachst die Fliigel mit konstanter Tiefe (Z = 1) 
herausgegriffen. Spater sollten die Ergebnisse auch auf einige andere Fliigel mit anderer Zu- 
spitzung erweitert werden. Die Ergebnisse dieser Zusammenstellung sind bereits frither in ein- 
zelnen Teilberichten mitgeteilt worden!??. Wahrend iiber Druckverteilungsmessungen an Pfeil- 
fliigeln bei symmetrischer Anstrémung heute bereits eine Anzahl von Ergebnissen vorliegt 4~7®® 
ist das Material iiber unsymmetrische Anstrémung aufBerst sparlich. An systematischen Unter- 
suchungen, die ja besonders wichtig sind, fehlt es noch vollkommen. Theoretisch ist die Berech- 
nung der Auftriebsverteilung von Pfeilfliigeln bei symmetrischer Anstrémung mit ausreichender 
Genauigkeit méglich. Es sind hier die Methoden von Weissinger, Falkner! und Mutterperl'? zu 
nennen. Fiir gerade Fliigel mit nicht zu kleinem Seitenverhaltnis ist die Methode von Multhopp™ 
die heute am weitaus gebrauchlichste. In neuester Zeit sind umfangreiche Rechnungen nach den 


* Die vorliegenden Messungen wurden zwischen 1942 und 1945 im fritheren Aerodynamischen Institut 
der Technischen Hochschule Braunschweig ausgefiihrt. (Leitung: Prof. Dr. H. Schlichting.) Sie wurden 
wihrend des Krieges in mehreren FB- und UM-Berichten der ZWB (Zentrale fiir wissenschaftliches Berichts- 
wesen) mitgeteilt, die heute nicht mehr zuganglich sind. Hier wird abschliefiend (in zwei Mitteilungen) ein 
zusammenfassender Bericht gegeben. Die I. Mitteilung erschien im Ing.-Arch. 18 (1950), S. 344. 

1 W. Jacobs, Druckverteilungsmessungen an Pfeilfliigeln konstanter Tiefe bei symmetrischer An- 
strémung. Untersuchungen u. Mitteilungen der D. V.L. Berlin-Adlershof, U. M. 2059 (1943). Bibliography 
of Technical Reports. Washington, P. B. 39148. 

2 W. Jacobs, Druckverteilungsmessungen an einem Pfleilfliigel konstanter Tiefe (y = 15) bei unsymme- 
trischer Anstromung. Untersuchungen u. Mitteilungen der D.V.L. Berlin-Adlershof, U. M. 2083. Biblio- 
graphy of Technical Reports. Washington, P. B. 39139. 

3 W. Jacobs, Druckverteilungsmessungen an zwei Pfeilfliigeln konstanter Tiefe (py = 30°, 45°) bei un- 
symmetrischer Anstrémung. Untersuchungen u. Mitteilungen der D.V.L. Berlin-Adlershof, U. M. 2110. 
Bibliography of Technical Reports. Washington, P. B. 39141. 

4 H. Luetgebrune, Druckverteilungsmessungen an Pfeilfliigeln. Forsch.-Ber. der Zentrale fiir wissenschaftl. 
Berichtswesen. Berlin-Adlershof, F. B. 1501 (1941). 

5 R. Géthert und H. Scholkemeier, Sechskomponentenmessungen und Druckverteilungsmessungen an 
einem Rechieckfliigel und einem Ellipsenfliigel. Forsch.-Ber. der Zentrale fiir wissenschaftl. Berichtswesen. 
Berlin-Adlershof, F. B. 1451 (1941). 

6 Th. Schwenk, Messungen an einem trapezférmigen Pfeilfliigel. Forsch.-Ber. der Zentrale fiir wissen- 
schaftl. Berichtswesen. Berlin-Adlershof, F. B. 1781 (1943). 

7 Thiel-Weissinger, Druckverteliungsmessungen an einem geraden und einem 35° riickgepfeilten Trapez- 
fliigel. Untersuchungen u. Mitteilungen der D.V.L. Berlin-Adlershof, U. M. 1293 (1944). 

8 Nicholac, van Dorn, Dr. Young, A comparison of three theoretical methods of calculating span load 
distribution on swept wings. N. A.C. A. Techn. Note No. 1476. 

9 M. Ingelman-Sundberg, Experimental determination of pressure distribution on a plane wing with 40° 
sweepback at low speed. Technische Hochschule. Stockholm, A. EK, R. O. TN 8. 

10 JT, Weissinger, Math. Nachr. 2 (1949), 5. 45. 

11 VY, M. Falkner, The calculation of aerodynamic loading on surfaces of any shaps. A. R. C. Rep. a. 
Memoranda 1910 (1943). 

12 W, Mutterperl, The calculation of span load distributions on swept back wings. N. A.C. A. Techn. 
Note 834 (1941). 

13 H, Multhopp, Luftf.-Forschg. 15 (1938), S. 153. 


Jacobs: Druckverteilungsmessungen an Pfeilfliigeln konstanter Tiefe usw. Ingenieur- Archiv 


7 Le 

& bit DT 
With 
Wa i} 


I ie: 
| —ili : 
Ha Tat He ! 
HAE Marit HE ig 


PILL ] 
{LL 7 L, Ll L/ L / L/ 


Hi | Hh 
LMT ip 
LY, 


(tie | 
inf ee Ba 
Al oy 
£7 


40 
G5 
0 — 


ee Wii 
“tt yi UiElh Li 
Lee 


Ul) 
GZ [LAs 
8 B 7g L } ty 


o 
723 
msi O=+H3° B=0° 


ppb ber ichtsskizze der  Fligel u und Lage der Meb- Abb. 2. Normalkraftverteilungen am Pfei eifltigel 10/15. 
und MeGanbohru ngen, Profil NACA 23 012. &=-+ 1,9; 5,7; 8,5; aa: 14,3°; 6 = 


|b oe AAA 


| fp HE Hy, Hila AM: > 
We ie Ui re tl fe 


tas ff 
g/l 


to 
Hy. fey 
OMB 7S” rete 


ey - ew 


AL en I Hy 


Ming ema 


Vi 


Vy, tee ob 


”y A 
Le 2 easel Oe 


ps, ge 
#57° ae 70° ‘erg 4 o 073 Zp 


sit 2 3 4 HeLa 


n 


| 
rTM) Ma 
hae 


i] / 

st i], Air 15 
L, ML Hy LT 4 
iy Vie ri Ley ye. YP 
7 Vag Bs LE he eke. ue8s soa? 


= 485° hes 
aT — 


We 


wan | es 
ie “heal Wy 
ALT Kor aA 
ALL TY. WE 


se er a AT 
it fs De My if WT) 
i hy My 
ETHYL 


fl4-s 
i123 


4 rae 


| 
k rn 
pay 


EL [Lb bh 7 


Mn | 


ag 
WEY. 


Wy thin 
Ve HME UL) HA] wy 


By 
Hy) 

2) ee a) 
1 a ae agate 


WE, Hd Li 
gat: 


wa, 
M7 a rhe 
iff: 129 yor of eg nee wg 


+m fo=10 
\, 


LEW 


ss 
Hh Zé 


Ht 
| 


Hae 


et Hh 
| al ay 


ITT 
i ai 


b 
Uf ra 
I 
e wit a=+3° B=20° 
A / Zo ; 
a 7 Wn PH 
ye 
r ve ed Wey Z 
itt: 124 e it th 70 ri 
a=4143° =I Wa $576 94; 1B 1g 
lo 


MC 
ae Wy ie 


LLf LY 
Ute LY 
=+57° f=20° leery Mf y 


Wal 
nna Tiny} 


85 


'Y Z 
Ge 
Uf x 
Yo’ a 
ied 
Hh Zo 


N. 
\ 
IX 
X. 
~. 


iy 


el eA 


86 


Jacobs: Druckverteilungsmessungen an Pfeilfliigeln konstanter Tiefe usw. Ingenieur- Archiv 


Mabistab 


00 
io\= 


GS 


Wp_p, A 
O- W 
MeP a. 
St 


holst 
60 [> | Us 


CMa «SAN 
_ Hers 
Oy Min 
\e lj 
Wh Yi Mian jie 
a Ee LN 
1 


yf 
Aiea 


Malistab AY 
15 i, I - | 
10F p,- My 9 Ge 

ae Ae : 
MN Ge ih fl 17 y 
. es Al 7 4 Te fe ° f @ Wy) LY by) ye jj t Vf f “; 16 Py WW 3 
ee qo Q-+t Cae = Ly) me 0 vor 9 Ce B-10° ome i Ze 
aif vi a chy Zp 4 A & oe 


e Hil eas Hi | 
ot ee aly | LT 
a Ns WN 
ee Yip ae y Droit yall Wy) ie ; ayes 
; 4 0 | aS a fe ay LW 
ee i Ee 


Abb. 3. Normalkraftverteilungen am Pfeilfligel 10/30. 6 =+1,9; 5,7; 8,53; 11,4; 14,3° B= 0, 10, 20° 
9 99'5 O05 oS 9 > =U, ? . 


X. Band 1951. 


Jacobs: Druckverteilungsmessungen an Pfeilfliigeln konstanter Tiefe usw. 87 


igen Methoden durchgefiihrt und miteinander verglichen worden!. Hierbei wurden Pfeil- 
inkel, Zuspitzung und Seitenverhaltnis in weiten Grenzen variiert. Als Ergebnis ist zu nennen, 
1B Falkners Methode die besten Resultate lieferte, allerdings mit groBem Rechenaufwand. Weis- 
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Zu Abb, 3. 


singers Methode ergab praktisch gleiche Ge- 
nauigkeit und ist daher infolge des geringeren 
Rechenaufwandes zu empfehlen. Mutterperls 
Methode stellt keine Verbesserung dar. 

Bei unsymmetrischer Anstrémung ist die 
Berechnung der Auftriebsverteilung bisher 
nur fiir gerade Fliigel méglich?. Wir wollen 
im folgenden eine einfache Naherungslésung 
hierfiir angeben, die jedoch recht brauchbare 
Ergebnisse liefert. 


2. Bezeichnungen und Bezugsgréfen. Es 
wurden die Bezeichnungen nach dem Norm- 
blatt Din L100 verwendet: 

x, y Flugzeugfeste Koordinaten, 

€, 7 Dimensionslose Koordinaten, § = 2 x/b; 
i= 24/b, 

p Druck am Profil, 

Pp = Druck auf Druckseite, 
ps = Druck auf Saugseite, 
Cn ry Ortlicher Normalkraftbeiwert, 
(Ac, értl) 2 == (c, értl — Cn srtly_9) Pf 
a) (es, értl —_ Cn drt y)p=0 
= zusatzlicher von der Pfeilung her- 
rihrender értlicher Auftriebsbeiwert, 

x, ¥; Koordinaten des geometrischen Neutral- 
punktes (= Flachenschwerpunkt einer 
Fiiigelhalfte), 

xx, ¥n Koordinaten des aerodynamischen 
Neutralpunktes (Abb. 12), 

yn, y-Koordinate des  aerodynamischen 
Neutralpunktes beim geraden Filiigel, 

| Abstandskoordinaten des 
ni fe | aerodynamischen Neutral- 
yn = ¥s— ¥n_ | punktes vom Flachenschwer- 

punkt (Abb. 12), 

(x), Neutralpunktsverschiebung des Pfeil- 
fliigels mit der Auftriebsverteilung des 
geraden Fliigels (Abb. 12), 

Neutralpunktsverschiebung 

des Pfeilfliigels infolgeAnde- 

rung derAuftriebsverteilung 

II: | gegeniiber dem ungepfeilten 

Fliigel (Abb. 12), 

Momentenachse: Achse durch 1/4-Punkt des 
Profilschnittes im Flachenschwerpunkt 
einer Fliigelhalfte = geometrische Neu- 
tralachse. 


A xN = XN— XNg 


1 Dr. Young, Theoretical additional span-loa- 
ding characteristics of wings with arbitrary sweep, 
aspect ratio and taper ratio. N. A.C. A. Techn. 
Note No. 1491. 

2 I.Weissinger, Jb. 1940 dtsch. Luftf.-Forschg., 
5.1145. 
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Profilbezeichnung: Die ersten beiden Ziffern 
geben die Zuspitzung in Zehnerprozenten 
und die beiden letzten die Pfeilstellung an. 
Bakes 10/308; L== 1.0% 10 =a 0°, 

Weitere Bezeichnungen siehe I. Mitteilung. 


3. Modelle. Fiir die Untersuchungen wurden 
vier Fliigel konstanter Tiefe mit den Pfeil- 
winkeln y= 0, 15, 30 und 45° benutzt (Abb. 1). 
Samtliche Fliigel haben gleiche Tiefe, gleiche 
Spannweite und gleiches Seitenverhaltnis (A = 
5,18). Als Profil in Schnitten parallel zur Sym- 
metrieebene wurde NACA 23012 gewahlt. Die 
Fliigelenden waren mit der iiblichen Abrundung 
versehen, wonach der jeweilige Abrundungs- 
radius gleich der halben értlichen Profildicke 
ist. V-Stellung und Verwindung waren nicht 
vorhanden. 

Die genauen Abmessungen der Fliigel sind: 

b’—0,770m, b6=—0,750m, F= 0,1146 m2, 
l= 0,150 m, Ain OpL Os 
Die Lage der MeBschnitte und Anbohrungen ist 
aus Abb. 1] zu ersehen. 


4, Durchfitthrung und Ergebnisse der Messung. 
Die Messungen wurden im Windkanal des Aero- 
dynamischen Instituts der Technischen Hoch- 
schule zu Braunschweig durchgefiihrt. Bei einer 
Windgeschwindigkeit von 40 m/sec betrug die 
Reynoldsche Zah] Re = vl/y = 4,2 - 10°. Neben 
den Druckverteilungsmessungen wurden auch 
stets Kraftmessungen durchgefihrt. 

a) Drei- und Sechskomponentenmes- 
sungen. Auf die Kraftmessungen wollen wir 
hier im einzelnen nicht eingehen, sondern nur 
auf die I. Mitteilung verweisen, in der die Ergeb- 
nisse einer groBen Systematik von trapezférmigen 
Pfeilfligeln mitgeteilt wurden. Die Neutral- 
punktlagen sind in Tabelle 1 und 2 enthalten. 

b) Druckverteilungsmessungen. Die 
Ergebnisse der Druckverteilungsmessungen in 
den verschiedenen Schnitten sollen hier im ein- 
zelnen nicht wiedergegeben werden. Bis auf die 
von der Randumstrémung beeinfluBten auBeren 
Schnitte zeigen die Kurven den normalen Ver- 
lauf. Fiir statische Berechnungen interessieren 
vor allem die Gesamtlastverteilungen [Normal- 
kraftverteilungen [(pp — ps)/q] iiber Fliigeltiefe 
und die integrierten Werte lings Spannweite, 
also die Auftriebsverteilungen, die natiirlich 
auch fiir Stabilitatsbetrachtungen wichtig sind. 
Die ersteren Werte sind in einer Schragprojek- 
tion itber dem Fliigel aufgetragen und zwar fiir 
die Schiebewinkel 6 = 0,10 und 20° (Abb. 2 
bis 4)1. Die Punkte gleicher Riicklage sind durch 


1 Die Normalkraftverteilungen fiir den Rechteck- 
fliigel siehe bei Méller: Druckverteilungsmessungen 
an einem Hoch- und Tiefdecker bei symmetrischer 
und unsymmetrischer Anstrémung. Techn, Berichte, 


Bd.11, Heft 5 (1944), Abb. 35 bis 38. 
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diinne Kurven miteinander verbunden, ebenfalls sind die maximalen Werte durch Strichelung 
verbunden. Bei allen Fliigeln ist in symmetrischer Anstrémung ein starker lokaler Kinbruch 
in der Auftriebsverteilung im Mittelschnitt vorhanden. Ferner ist aus dieser Auftragung deut- 


lich zu entnehmen, daB bei den gepfeilten Fliigeln mit zunehmendem Schiebewinkel die vorgehende | 


Fligelhalfte stirker tragt als die riickliegende. Dieser Unterschied ist um so starker, je gréBer 
der Pfeilwinkel ist. Beim geraden Fliigel dagegen erhalt der zuriickgehende Fliigel etwas mehr 
Auftrieb als der vorgehende, abgesehen von den Randgebieten, an denen infolge der Rand- 
umstrémung an der vorderen Kcke eine starke positive und an der hinteren eine starke negative 
Auftriebsspitze auftritt. 
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Abb. 5. Auftriehsverteilungen am Rechteckfligel 10/00. « =+ 1,8; 5,6; 8,5; 11,4; 14,3°; 8 = 0, 5, 10, 15, 20°. 


Die aus den Druckverteilungsmessungen erhaltenen Lastverteilungen langs Spannweite sind 
in Form der 6értlichen Normalkraftbeiwerte ¢, 5, in Abb. 5 bis 8 zusammengestellt. 

Betrachten wir zunachst einmal die Auftriebsverteilungen bei symmetrischer An- 
strémung (Abb. 14 bis 17), so sehen wir, da die Auftriebsverteilung mit wachsender Pfeil- 
stellung des Fliigels eine wesentliche Anderung erfahrt, und zwar tritt eine Verlagerung des Auf- 
triebes nach den Fliigelenden hin auf. Wahrend beim Rechteckfliigel eine kontinuierliche Ab- 
nahme des Auftriebes nach auBen hin stattfindet, wird bei gréBeren Pfeilwinkeln der Auftrieb 
in der Mitte so stark abgebaut, da nach aufBen hin sogar eine Zunahme festzustellen ist. Im 
Camax-Bereich wird dies zur Folge haben, dafs der Beginn des AbreiBens der Strémung sich mit 
wachsender Pfeilstellung nach auBen verschiebt. Unterstiitzt wird die Ablésung im auBeren Teil 
noch durch den ,,Grenzschichteffekt‘*. Infolge der Pfeilstellung des Fligels tritt ein Druck- 
gradient in Querrichtung auf, dem vor allem das Grenzschichtmaterial infolge seiner geringen 
Tragheitskrafte folgt. Beim riickwarts gepfeilten Fliigel ist auf der Saugseite infolge der Druck- 
abnahme nach aufen eine Anhaufung des Grenzschichtmaterials im auBeren Gebiet festzustellen, 
was die Abrei®gefahr in diesem Teil erhéht. Die von der Pfeilung herrithrende Auftriebsanderung 
und der Grenzschichteffekt wirken also im gleichen Sinne. 
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Infolge der durch Pfeilung hervorgerufenen értlichen Auftriebsinderungen wird eine Ver- 
schiebung des Lastschwerpunktes einer Fliigelhalfte und damit der Neutralpunktlage gegeniiber 
dem geraden Fliigel bewirkt. Wir wollen diese Neutralpunktverschiebung durch Pfeilstellung 
aus unseren Druckverteilungsmessungen bestimmen. 

Im ersten Teilbericht hatten wir die Beziehung abgeleitet [I. Mitteilung, Gleichung (5) | 


Axy Ayn b 
Siew pare ee (1) 


wobei Axy/l die Neutralpunktverschiebung infolge Auftriebsinderung des Pfeilfliigels gegeniiber 


T ir 
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09-9=19°. | g "orth 
a =0=49 06: See 
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Abb. 6. Auftriebsverteilungen am Pfeilfliigel 10/15. «0 =-+1,9; 5,7; 8,5; 11,4 14,3°; B = 0, 5, 10, 15, 20°. 


dem ungepfeilten Fliigel bedeutet (Abb. 12). Als gesamte Neutralpunktverschiebung gegeniiber 
dem geraden Fliigel ergibt sich, wenn wir die einzelnen Kinfliisse getrennt auffiihren, 


oa ied? YN VN \ 8 
Fee ner Se (ia aT, eee oe: 21 '8P- (2) 
In dieser Beziehung stellt der erste Term die Verschiebung des Flachenschwerpunktes infolge 
der Pfeilstellung dar, was gleichbedeutend ist mit der Neutralpunktsverschiebung bei konstanter 
Cn ory Verteilung; der zweite Term umfaft die Verschiebung, die dadurch entsteht, daf die Auf- 
triebsverteilung des Rechteckfligels von der konstanten ¢,5,,-Verteilung abweicht; der dritte 
Term schlieBlich gibt die weitere zusatzliche Verschiebung infolge der Abweichung der Auftriebs- 
verteilung des Pfeilfliigels von derjenigen des ungepfeilten Fligels an. 

Aus der Lage der Lastschwerpunkte einer Fliigelhalfte kénnen wir samtliche oben abgeleiteten 
GréBen bestimmen. Den Abstand des Lastschwerpunktes von der Symmetrieebene (yy) erhalten 
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wir aus der gemessenen Auftriebsverteilung auf graphischem Wege, indem wir bilden 


1 
nnd 
= Yn wk if pales (3) 
Cues prowa ai i. | 
Sen dy 4|| 
M4 . 
Diese fiir alle vier Fliigel ermittelten Werte sind in Tabelle 2 enthalten. Hierbei wurde der Mittel- | 
wert gebildet fiir die Anstellwinkel « = 1,9; 5,7; 8,5 und 11,4°. Der gréBte Anstellwinkel wurde | 
fortgelassen, weil in diesem Falle schon értliche Abrei®erscheinungen auftreten, und infolge- 
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Abb. 7. Auftriebsverteilungen am Pfeilfliigel 10/30. B =+ 1,9; 5,7; 8,5; 11,4; 14,3°; 6B = 0, 5, 10, 15, 20°. 


dessen das Ergebnis falschen kénnen. Aus der Gri®e yn/56 und aus der Pfeilung lassen sich 
dann nach obigen Bezichungen die Neutralpunktverschiebungen berechnen. Diese sind in Ta- 
belle 2, 3 aufgefiihrt. Fiir ¥n/4 6 ergibt sich ein Anwachsen mit der Pfeilstellung, der Lastschwer- 
punkt wandert also nach auBen. Gegeniiber dem Rechteckfliigel betragt diese Wanderung beim 
Pfeilfliigel p = 45° ungefahr 3% der Halbspannweite des Fliigels (Aynx/$b = 0,031). Dies ent- 
spricht einer zusatzlichen Neutralpunktverschiebung infolge Anderung der Auftriebsverteilung 
von 8% der Fliigeltiefe (Axy/l= 0,08). Die Auftriebsverteilung des Rechteckfliigels darf daher 
fiir die Berechnung der Neutralpunktlage des Pfeilfliigels unter keinen Umstinden zugrunde 
gelegt werden. 

Bei unsymmetrischer Anstrémung, also beim Schieben, erfahren die Auftriebsvertei- 
lungen der Fliigel recht betrachtliche Anderungen, die stark von der Pfeilung abhangen (Abb. 5 
bis 8). Abgesehen von den Fliigelenden zeigt der Rechteckfliigel an der vorgehenden Fliigel- 
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halfte eine starkere Abnahme der Auftriebsverteilung als an der zuriickgehenden, wie Abb. 5 


zeigt, auf dem die Verteilun 
p=, 10, 15 und 20° zusammengetragen sind. Anders | 
Fliigeln. Bei dem Pfeilfliigel gy = 15° (Abb. 6) 


gen jeweils fiir einen festen Anstellwinkel bei den Schiebewinkeln 
iegen die Verhaltnisse bei gepfeilten 
ist schon deutlich sichtbar, da®B der Auftriebs- 


verlust beim Schieben auf der riickgehenden Hiilfte stirker ist als auf der vorgehenden. Noch 
stiirker tritt der Einflu8 der Pfeilung bei dem Fliigel gp = 30° hervor (Abb. 7). Hier ist der Auf- 
trieb an der vorgehenden Fliigelhalfte fast unabhangig vom Schiebewinkel §, an der riickgehenden 
tritt eine betrachtliche Abnahme auf. Setzen wir den GesamteinfluB8 am schiebenden Pfeilfliigel 
zusammen aus dem des schiebenden Rechteckfliigels und einem zusitzlichen Pfeileffekt, so ist 


bei dem 30°-gepfeilten Fliigel die Zunahme des Auftriebs an der vorgehenden Fliigelhalfte in- 
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Abb. 8. Auftriebsverteilungen am Pfeilfliigel 10/45. 
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folge des Pfeileffektes ungefahr gleich der Abnahme des Auftriebs beim geraden Fliigel, so daB 
sich diese beiden Einfliisse nahezu aufheben. An der riickgehenden Halfte tritt eine besonders 
starke Abnahme auf, da sich hier die Einfliisse unterstiitzen, denn auch der Rechteckfliigel zeigt 
an der riickgehenden Halfte eine Abnahme, die allerdings schwicher ist als an der vorgehenden. 

Der Fliigel mit der Pfeilung g — 45° zeigt schon ein betrachtliches Uberwiegen des Pfeil- 
effektes, d. h. es tritt auch an der vorgehenden Halfte eine merkliche Zunahme des Auftriebs mit 
dem Schiebewinkel ein (Abb. 8). Die Abnahme auf der riickgehenden Halfte wird noch starker 


als bisher. 


Den reinen EinfluB der Pfeilung gegeniiber dem Rechteckfliigel erhalten wir, wenn wir die 
Differenzen der Ac, 5,-Werte des Pfeilfliigels gegeniitber denen des Rechteckfliigels bilden: 


(Ach s)p¢ = (Cn drtl ka ag Cn drtig 9) Pf — (Cn set aia Cn értty_»)p=0 : 


(4) 


Diese ausschlieBlich von der Pfeilung herrithrenden zusatzlichen Normalkraftheiwerte sind fiir 
die Pfeilfliigel ermittelt worden, wobei gleiche geometrische Anstellwinkel bei Pfeilfliigel und 
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Abb. 9. Zusiitzliche Auftriebsiinderungen beim Schieben, herriithrend von der Pfeilstellung. Pfeilfliigel 10/15. 0 =+1,9; 5,7; 8,5; 11,4; 1 
B = 5, 10, 15, 20°. 


Messung , heorie 


G8 


= 


vorgehend rickgehend vorgehend Pickgehend—~ 
Abb. 10. Zusitzliche Auftriebsinderungen beim Schieben, herrithrend von der Pfeilstellung. Pfeilfliigel 10/30. ~« = +1,9; 5,7; 8,5; 11,4; 14 


B = 5, 10, 15, 20°. 
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- Rechteckfliigel zugrunde gelegt wurden. Die Ergebnisse sind fiir die drei Pfeilfliigel] bei 6 = 5, 10, 
15 und 20° in Abb. 9 bis 11 aufgetragen. Der Pfeileffekt wirkt sich auf beiden Fliigelhalften un- 

| gefahr gleich aus, die Kurven verlaufen nahezu antimetrisch. Sie sind proportional dem Schiebe- 
} winkel 6 und dem Anstellwinkel. AuBerdem wachst der EinfluB natiirlich mit der Pfeilung 
) selber. Eine ziemlich sprunghafte Anderung des Pfeileffektes ist in Fliigelmitte zu beobachten. 
i Die Werte (4 Cn ort) Pf gehen hier in einem ziemlich schmalen Gebiet von positiven Werten des 
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B11. Zusatzliche Auftriebsianderungen beim Schieben, herriithrend von der Pfeilstellung. Pfeilfliigel 10/45. ~ =-+1,9; 5,7; 8,65 11,55 14,4°; 
| = Oy 105 os 20. 


vorgehenden Fligels zu negativen des zuriickgehenden Fliigels ttber. Im héheren Anstellwinkel- 
bereich verlaufen die Kurven nicht mehr antimetrisch. Die Abnahme des Auftriebs an der zuriick- 
gehenden Halfte ist wesentlich geringer als die Zunahme an der vorgehenden Fliigelhalfte. Be- 
merkenswert an der zusatzlichen Auftriebsverteilung infolge der Pfeilung ist die ziemlich kon- 
stante Verteilung in einem groBen Bereich der Halbspannweite, ahnlich dem Einfluss der Ka 
Stellung eines Fliigels auf die Auftriebsverteilung beim Schieben. Diese Tatsache wird uns, wie 
spater naher ausgefiihrt, zu einer naherungsweisen Berechnung der Auftriebsverteilung eines 


schiebenden Pfeilfliigels fiihren’. 


1 Da die groBen Auftriebsspitzen an den Fliigelenden nach der Messung nur ungenau erfalit werden 
konnen, beim Rollmoment aber wegen des grofien Hebelarmes stark eingehen, wurde auf die Ermittlung 
2 . . 
des Rollmomentes aus den Druckverteilungen verzichtet, 
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Das durch die antimetrische zusatzliche Auftriebsverteilung erzeugte positive Rollmoment 
kann im Gebiet des Auftriebsmaximums durch AbreiBerscheinungen ins Gegenteil umschlagen. 
Diese Gefahr besteht vor allem bei starker gepfeilten Fliigeln, an denen mit zunehmendem Schiebe- 
winkel der vorgehende Fliigel erhéhten Auftrieb bekommt und daher leicht den maximalen 
Auftriebsbeiwert iiberschreiten kann. Es ist hierbei jedoch auch noch der Grenzschichteffekt 
zu beriicksichtigen, der vor allem bei gréBeren c, -Werten infolge des starken Druckgradienten 
in Querrichtung des Fliigels ein starkes AbflieBen von Grenzschichtmaterial bewirkt, was in ge- 
wissen Teilen des Fligels erhéhte Ablésegefahr zur Folge hat. Betrachten wir einen schiebenden 
riickwarts gepfeilten Fliigel, so wird die Schraganstrémung der vorgehenden Fliigelhalfte ver- 
mindert, die der zuriickgehenden dagegen verstarkt. Der Grenzschichteffekt wirkt also an der 
riickgehenden Fliigelhalfte mit wachsendem Schiebewinkel im Sinne einer Verstarkung der Ab- 
lésung, wahrend er an der vorgehenden im Sinne einer Verminderung wirkt. 


5. Vergleich von Messung und Theorie. a) Auf- 
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Abb. 13. Dimensionslose Zirkulationsverteilungen fiir die Pfeilfliigel 10/00, 10/15, 
Abb. 12. Skizze zur Nentralpunktberechnung. 10/30 und 10/45 nach den Verfahren von Weissinger und Multhopp. 


nungen vergleichen. Fiir die Berechnung der Auftriebsverteilung von Pfeilfliigeln wollen wir 
folgende Verfahren betrachten. 


1. Verfahren von Multhopp. Dies baut auf der Prandtlschen Traglinientheorie auf unter 
Anbringung eines bestimmten Korrekturfaktors. 


2. ,,Tragflachenverfahren‘‘ von Weissinger. Hierbei wird eine flachenhafte Zirkulations- 
verteilung am Fligel zugrunde gelegt. 


3. ,,Traglinienverfahren‘* von Weissinger. In diesem Falle wird wieder die tragende Linie 
beniitzt. Sie wird in '/, der Fliigeltiefe gelegt, und die Zirkulationsverteilung wird nun so be- 
stimmt, daB der vom Wirbelsystem induzierte Abwind auf der 3 1/4-Linie des Fliigels gleich der 
Normalkomponente der Anblasgeschwindigkeit ist. 

An weiteren Verfahren waren zu nennen das von Falkner und das von Mutterperl. Das Ver- 
fahren von Falkner liefert ungefahr dieselbe Genauigkeit wie Weissinger, ist jedoch wesentlich 
miihsamer, wahrend das Verfahren von Mutterperl keine Verbesserung darstellt. 

Weissinger hat gezeigt, daB zwischen dem Tragflachenverfahren und dem Traglinienverfahren 
im Ergebnis praktisch kein Unterschied besteht, und man daher ohne Bedenken das einfachere 
Traglinienverfahren anwenden kann. Am einfachsten ist das Verfahren von Multhopp. Es zeigt 
sich jedoch, daB dieses Verfahren fiir gréBere Pfeilwinkel unzureichende Ergebnisse liefert. 

Die dimensionslosen Zirkulationsverteilungen y (7) fiir die verschiedenen Pfeilfliigel, berechnet 
nach den Verfahren von Weissinger und Multhopp, sind in Abb. 13 enthalten. Dabei bedeutet 
y=T'/bv. Die Rechnung von Weissinger gilt fiir cj. = 22, wahrend nach Multhopp mit 
Ca~ — 5,20 gerechnet wurde. Bei dem Verfahren von Multhopp wurde der Korrekturfaktor 
x = 1 gestezt. Aus dem Vergleich von Theorie und Messung ist folgendes zu entnehmen: 
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1, Auftriebsanstieg. Nach Multhopp ist der Proportionalitatsfaktor zwischen Zirkulation 
und effektivem Anstellwinkel unabhangig vom Pfeilwinkel, also dc,/0%—= konst. Die Theorie der 

tragenden Flache nach Weissinger liefert jedoch einen Abfall des Gesamtauftriebes und damit 

des Auftriebsanstieges mit wachsender Pfeilung. Die nach Weissinger ermittelten Werte dc,/00 

sind mit in Tabelle 1 enthalten. Es ergibt sich nach der Messung eine Kleinere Abnahme von 

) Oc,/0% als nach der Theorie, und zwar betragt sie ungefahr 60°, des theoretischen Wertes. 

2. Auftriebsverteilung. Wichtiger als die Frage : 

) nach dem Gesamtauftrieb ist die ae nach der eas Cee oe ete eae 

: : : ; es ezugspunkt: 1/4-Punkt 
| der Auftriebsverteilung, denn diese ist mafigebend fiir des Profilschnittes im Flichenschwer- 
» die Lage des Neutralpunktes und fiir das AbreiBverhalten punkt einer Fliigelhalfte. 


| bei groBen Anstellwinkeln. Der Vergleich der theoretischen Eos ie i Dia 

| Auftriebsverteilung mit der aus der Druckverteilung er- a ak as ar. 5 ey 

/ mittelten ist in Abb. 14 bis 17 enthalten. Es sind hierbei Ee | mes es 

' die nach Weissinger und Multhopp errechneten Ver- 0 0 3,89 3,90 

_teilungen jeweils auf die entsprechenden c,-Werte der ¥ ae oe ae 
Messung umgerechnet worden. Wie aus Abb. 14 ersicht- 45 0.040 333 3.00 


lich, zeigt schon der Verlauf der Kurven beim Rechteck- 
 fliigel merkliche Unterschiede. Die Auftriebsverteilung nach Multhopp ist vélliger als die nach 
Weissinger. Die Theorie von Weissinger liefert eindeutig eine bessere Ubereinstimmung mit den 
Druckverteilungsmessungen als die von Multhopp. Dieser Unterschied kann also durch Nicht- 
| beriicksichtigung der flachenhaften Zirkulationsverteilung erklart werden. Er ist natiirlich auch 
| in der Berechnung der Auftriebsverteilung am Pfeilfliige] enthalten. Hinzu kommen noch andere 
_ Unterschiede der beiden Berechnungsverfahren. 

Der Vergleich fiir die Pfeilfliigel ist in Abb. 15 bis 17 enthalten. Die Auftriebsverteilung 
nach Weissinger weist bei allen Pfeilwinkeln gute Ubereinstimmung mit der Messung auf. Nur 
in der Fliigelmitte treten fiir die groBen Pfeilungen m = 30° und 45° starkere Abweichungen auf; 
| die Theorie liefert hier zu kleine Werte. Bei groBen Anstellwinkeln ist natiirlich infolge von Ab- 
reiBerscheinungen keine gute Ubereinstimmung zu erwarten. Die Theorie von Multhopp liefert 
in allen Fallen weniger gute Ubereinstimmung mit der Messung als die von Weissinger. Nach 
der Messung ist die Verlagerung des Auftriebes nach auBen mit der Pfeilung nicht so stark wie 
nach der Theorie von Multhopp. 


Tabelle 2. Lage des Lastschwerpunktes einer Fliigelhdlfte und Neutralpunktverschiebung, hervorgerufen durch 
Anderung der Auftriebsverteilung; Vergleich von Theorie und Messung. 


Dene entel : Mittelwert von recht- 
ruckverteilungs- 


Kraftmessung Theorie Weissinger? Theorie Multhopp? eckiger und elliptischer 
Ne Ata) Auftriebsverteilung 
Oe a eS 
Paar aes, | Love ods oy 1) ae | AI Ae IN| Ya Ae We) Soe) aay 
b/2 =| 6/2 l b/2 b/2 ji b/2 b/2 l b/2 “b/2 et! b/2 b/2 1 
0 | 0,454 | 0 0 0,451 | 0 0 0,463 | 0 0 0,463 0 0 
15 | 0,462 | 0,008 | 0,006 | 0.470 | — | — | 0,461 | 0,010 | 0,007 | 0,476 | 0,013 | 0,009 | 0,463 0 0 
30 | 0,472 | 0,018 | 0,027] 0,478 | — | — | 0,470 | 0,019 | 0,028] 0,491 | 0,027] 0,039 | 0,463 0 0 
45 | 0,485 | 0,031 | 0,078 | 0,483 | — | — | 0,484 | 0,034) 0,085 | 0,512 | 0,047 | 0,118 | 0,463 0 0 


3. Neutralpunktlage. Die genaue Kenntnis der Auftriebsverteilung ist besonders wichtig 
fiir die Vorausbestimmung der Lage des Neutralpunktes. Wir wollen untersuchen, wieweit die 
nach der Messung und nach der Theorie ermittelten Neutralpunktlagen iibereinstimmen. Die 
Ermittlung geschieht auf die Weise, daB wir zunachst die Abstiande der Lastschwerpunkte der 
Fliigelhalfte bestimmen und hieraus dann die Neutralpunkte unter der Annahme, daf} diese Last- 
schwerpunkte auf der 1/4-Linie des Fliigels liegen. Die Lage der Lastschwerpunkte ist aus Ta- 
belle 2 zu entnehmen. Hierbei ist nach der Messung der Mittelwert der Auftriebsschwerpunkte 
fiir «= +1,9; 5,7; 8,5 und 11,4° gebildet worden. Die Anstellwinkel « — — 3,9 und ~+- 14,3° 
sind nicht beriicksichtigt worden, um das Ergebnis nicht zu verfalschen. Die Auswertung wurde 
graphisch vorgenommen. Wir sehen aus der Tabelle 2, daB die Ubereinstimmung mit der Theorie 
von Weissinger recht befriedigend ist. Die Theorie von Multhopp liefert in allen Fallen zu hohe 
Werte. Die zum Vergleich noch herangezogenen Ergebnisse der Kraftmessungen, die aus unserer 
I. Mitteilung entnommen sind, zeigen geringe Unterschiede gegenitber den Druckverteilungs- 
messungen. Dies kann wohl darauf zuriickgefiihrt werden, daB der Verlauf der Momentenkurven 
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sind noch die Werte fiir konstante Cnou-Verteilung und fiir den haufig beniitzten Mittel- 


nicht ganz geradlinig ist und daher die Ableitun 
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Aus der Lage des Lastschwerpunktes kénnen wir jetzt sofort die Lage des Neutralpunktes be- 
rechnen. Die Ergebnisse fiir die verschiedenen theoretischen Auftriebsverteilungen sind in Tat 
belle 3 den Druckverteilungsmessungen gegeniibergestellt. Der Neutralpunkt ist hierbei gerechaal 
von 1/4 in Fliigelmitte nach Abb.12. Zunachst ist als einfachster Fall wieder konstante 
Casrn7 Werteilung angenommen. Die hieraus folgenden Neutralpunktlagen sind identisch mit 
den Riicklagen der Flichenschwerpunkte. Als hiufig beniitzte Naherung ist auch der Mittel- 
wert zwischen rechteckiger und elliptischer Auftriebsverteilung herangezogen 
worden. Weiter sind die Neutralpunktlagen aufgefiihrt, die sich unter Zugrundelegung der ge-) 


y 


ergeben. Die Ergebnisse nach der Theorie von Multhopp und nach der Theorie von Weissinger 
sind neben den experimentellen Werten nach der Druckverteilungsmessung und Kraftmessung) 
ebenfalls in Tabelle 3 enthalten. Fordert man eine Genauigkeit von 1% der Fliigeltiefe fiir die} 


Tabelle 3. Vergleich der Neutralpunktlagen nach Theorie und Messung. Werte xy/l 
(xn gerechnet von 1/4 in Fliigelmitte nach Abb. 12). 


Rechnung Messung 
; mo (2) (8) (4) (5) (6) (7) 

Mittelwert umgerechnet mit nach Druck- | 

Gp” konst. cy g,yJ. |von rechteckiger| der gemessenen nach Theorie nach Theorie nach verteilungs- 

Verteilung und elliptischer | Verteilung des Multhopp? Weissinger! Kraftmessung messung 
Verteilung Rechteckfliigels 
I 0 0 0 0 0 0 0 0 

15 0,335 0,310 0,304 0,319 0,309 0,315 0,310 

30 0,722 0,669 0,657 0,709 0,678 0,690 0,683 

45 1,250 1,158 1,136 1,279 1,210 1,210 1215 


Werte Axy/1 = I—I = Neutralpunktverschiebung infolge Abweichung der Auftriebsverteilung 


von der konstanten ¢n éril- Verteilung. 


(1)' 


LEO 0 0 0 0 0 0 
15 —0,025 —0,031 —0,016 —0,026 —0,020 —0,025 
30 —0,053 —0.065 —0,013 —0,044 —0,031 —0,039 
45 —0,092 —0,114 +0,029 —0,040 —(0,040 —0,035 


Werte Axy/l = I — T(3): = Neutralpunktverschiebung infolge Abweichung der Auftriebsverteilung 
von derjenigen des ungepfeilten Rechteckfliigels. 


lil 0 0 0 0 0 0 0 
15 0,031 0,006 0,015 0,005 0,011 0,006 
30 0,065 0,012 0,052 0,021 0,033 0,026 
45 0,114 | 0,022 0,143 0,074 0,074 0,079 


Vorausberechnung der Neutralpunktlage, so zeigt der Vergleich zwischen Theorie und Messung, 
daB diese Forderung nur von dem Berechnungsverfahren von Weissinger erfiillt wird. Das Ver- | 
fahren von Multhopp liefert schon zu groBe Werte. Erfreulicherweise ergibt der Mittelwert zwi-| 
schen rechteckiger und elliptischer Auftriebsverteilung bis zu 30° Pfeilwinkel brauchbare 
Werte, die in den meisten Fallen fiir Uberschlagsrechnungen ausreichen®. In Tabelle 3 II 
ist weiter noch die Verschiebung des Neutralpunktes infolge der Abweichung der Auftriebs- | 
verteilung von der konstanten ¢,,,-Verteilung enthalten (I —I,1). Den Einflu8 der Auf- 
triebsinderung des Pfeilfliigels gegeniiber dem Recheckfliigel zeigt Tabelle 3111 (I — I, 3). 

Zusammenfassend kénnen wir auf Grund dieser Ergebnisse feststellen, daB® die Theorie von 
Weissinger die beste Annaherung an die Messung liefert. Die Unterschiede bewegen sich hinsicht- 
lich der Neutralpunktlage innerhalb von 1% der Fliigeltiefe. Die Theorie von Multhopp liefert 
groBe Werte fiir die Abstainde der Lastschwerpunkte und damit bei gré®eren Pfeilwinkeln auch 
fiir die Neutralpunktlagen. 


1 Vgl. Note 10 von S. 45. 

® H, Multhopp, Die Anwendung der Tragfliigeltheorie auf Fragen der Flugmechanik bei uns ymme- 
trischer Anstromung. Bericht S, 2 der Lilienthalgesellschaft fiir Luftf.-Forschg. 1939. 

* DaB der Mittelwert von rechteckiger und elliptischer Auftriebsverteilung trotz der groBen Abwei- 
chungen beim Lastschwerpunkt (Tabelle 2) fiir den Neutralpunkt doch so gute Werte liefert, liegt daran, 
dai einerseits die grofen Unterschiede im Lastschwerpunkt bei kleinen Pfeilwinkeln wenig ausmachen, 
andererseits bei grofen Pfeilwinkeln der Lastschwerpunkt nahezu richtig ist. 
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b) Theoretische Naherung fiir die Berechnung der Auftriebsverteilung eines 
schiebenden Pfeilfliigels. Bei unsymmetrischer Anstrémung ist die Berechnung der 
Auftriebsverteilung bisher nur fiir den geraden Fliigel méglich. Die Methode von Weissinger 
liefert hierfiir recht brauchbare Ergebnisse!. Fiir die theoretische Berechnung der Auftriebs- 
verteilung eines schiebenden Pfeilfliigels liegt noch kein brauchbares Verfahren vor. Es 
soll im folgenden der Versuch einer naherungsweisen Ermittlung der Auftriebsverteilung eines 
schiebenden Pfeilfliigels gemacht werden. 

Wir gehen bei der Ableitung der Naherung aus von der Voraussetzung, daB jede der beiden 
Fligelhalften eines Pfeilfliigels einen schiebenden Fliigel darstellt, und daher die Anblasung auf 
beiden Halften in eine wirkungslose Tangentialkomponente und eine wirksame Normalkompo- 
nente zerlegt werden kann. 

Bezeichnen wir die Grundgeschwindigkeit mit V, und die wirksame Normalgeschwindigkeit 
mit V,, so gilt fiir einen Fliigel mit dem Pfeilwinkel ~ (Abb. 18) fiir den Schiebewinkel 6 = 0° 

Vi,== Voces B-. vy) 
Betrachten wir jetzt einen schiebenden Pfeilfliigel, so vorgehende Hilfte | guriickgehende Maifte 
erhalten wir fiir die Normalkomponenten der beiden 
Fliigelhalften 

Vorgehende Fliigelhalfte: V,,—= V, cos (py — f) , Vn 

Zuriickgehende a5 a Ve bcos M .) , I= 00s (pf) Nc 
Als Verhaltnis der Anderung der wirksamen Anblas- pe 
geschwindigkeit, herriihrend vom Schieben, gegeniiber 
der Normalkomponente bei 6 = 0 folgt fiir die vor- 
gehende Fligelhalfte 


4,=Grundgeschwindigkelt 


5 Pin 5 (Vn)p=o = AY, = V cos (p — B) — Vo cos p Y4,=wirksame Normalkomponente 
) B=0 ( y) p=0 V, cos p aer Geschwintlgkeit 
a Bee (p =f) = Ag Abb. 18. Skizze zur Bestimmung der wirk- 
an cos samen Anblasgeschwindigkeit am schiebenden 
Pfeilfliigel. 


oder entwickelt nach kleinen Werten von [ 


AVn 2 
cae - +B tgp. (5a) 
Ganz entsprechend ergibt sich fiir die zuriickgehende Fligelhalfte: 
AVn : Sb 
ER or ea (5b) 


‘Die Anderung der wirksamen Normalkomponente kann man nun in erster Naherung ersetzen 
durch eine aquivalente Anstellwinkelanderung. 
Die Prandtlsche Integralgleichung als Ausgangsgleichung fiir die Berechnung der Auftriebs- 


verteilung des geraden Fliigels lautet im allgemeinen Fall: 


+b/2 
1 LN Ci le cea ZEN) a oe 6 
4 Vo i Che me gee al Cho Vol(y) as(¥) - (6) 
—b/2 


Diese Gleichung gilt streng genommen nur fiir gleichmafige Anstrémungsgeschwindigkeit V 
langs der ganzen Spannweite ; man kann jedoch auch noch bei schwach inhomogenen Strémungen 
mit ihr rechnen, wie Bausch? gezeigt hat. Dies wurde auch durch spatere theoretische und 
-experimentelle Untersuchungen bestatigt’. 

Nehmen wir also V als veranderlich an langs der Spannweite, so gilt 


| +6/2 

| we / dn er a EEG). a 

| Ce Or yy ce Vy) tg) 
—b/2 


| 1 W. Jacobs, Sechskomponentenmessungen an drei Pfeilfliigeln. Forsch.-Ber. der Zentrale fiir wissen- 
ischaftl. Berichtswesen. Berlin-Adlershof, F. B. 1629 (1942). Bibliography of Technical Reports. Washington 
IP. B. 36166. 

2 K. Bausch, Luftf.-Forschg. 16 (1939), S. 29. 

3 H. Schlichting und W. Jacobs, Jb. 1940 dtsch. Luftf.-Forschg., 5. I 81. 
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oder 
1 +P aT ) dy’ 2I'(y) V(y) AV | 
= ( y y [te x a — ¥) —— 1 Se . fl | 
4 Vo | dy’ y—y' | C0 Vo! (Y) "8 Vo 25(0)( iE r,) ( ) 
—b/2 


Es wird also die Auftriebsverteilung eines Fliigels mit der geometrischen Anstellwinkelver-, 
teilung a,(y) in einer inhomogenen Geschwindigkeitsverteilung gleichgesetzt einer Auftriebsver- 
teilung descelben Fliigels in homogener Strémung mit der geometrischen Anstellwinkelverteilung, 

a,(y)* (1+ AV/V,). Diesen Winkel nennen wir den ,,Aquivalenten geometrischen Anstellwin-, 
kel, Die GréBe «, 4V/V, kénnen wir als eine zusatzliche iquivalente geometrische Verwindung} 
ansehen. Analog kénnen wir die Auftriebsverteilung erhalten durch Uberlagerung der Grund- 
auftriebsverteilung und der infolge der aiquivalenten Verwindung erzeugten Auftriebsverteilung. . 

Zur Berechnung der Auftriebsverteilung eines schiebenden Pfeilfliigels setzen wir jetzt fiir - 
AV|V, die Zusatzgeschwindigkeit A V,,/ Vng-o an der vorgehenden bzw. zuriickgehenden Flii-} 
gelhalfte ein. Wir erhalten auf diese Weise eine ,,fiquivalente Verwindung“ des schiebenden Pfeil- 
fliigels. Auf Grund dieser aquivalenten Verwindung 1a8t sich dann die zusatzliche Auftriebs-} 
verteilung des schiebenden Pfeilfliigels bestimmen, z. B. nach der Methode von Multhopp. Eine | 
Verfeinerung wiirden wir erhalten, wenn die Auftriebsverteilung fiir die obige Verwindung nach | 
der Methode von Weissinger fiir Pfeilfliigel berechnet wird, da der Zusammenhang zwischen Ge-| 
schwindigkeitsinderung und dquivalentem Anstellwinkel derselbe bleibt. Die gesamte Auftriebs- | 
verteilung wiirde dann aus der Grundauftriebsverteilung und der zusatzlichen Auftriebsverteilung | 
folgen. 

In unserem Falle wurde fiir die Pfeilfliigel nur die zusatzliche Auftriebsverteilung nach der 
einfachen Methode von Multhopp berechnet. Beriicksichtigt wurden in den Gleichungen (5a) 
und (5b) nur die linearen Glieder in f. Es wurde also gerechnet mit 


AVn . 
(es, Se (8) 


Der Vergleich mit den Messungen ist in Abb.10u 11] enthalten. Es ergibt sich im groBen und gan- 
zen eine recht brauchbare Naherung. Bei grofen Anstellwinkeln treten am riickgehenden Fliige! 
stirkere Abweichungen auf, wahrend die vorgehende Fliigelhalfte immer noch gute Uberein- | 
stimmung zeigt. Inwieweit die Ergebnisse auch fiir andere Fliigelformen zutreffen, 1aBt sich 
nicht mit Sicherheit sagen; aber es ist zu erwarten, daB die Ubereinstimmung ahnlich sein wird. 


6. Zusammenfassung. Es wird iiber Druckverteilungsmessungen an Pfeilfliigeln konstanter 
Tiefe mit den Pfeilungen m= 0, 15, 30 und 45° bei symmetrischer und unsymmetrischer An- 
strémung berichtet. Hieraus werden die Auftriebsverteilungen am Fliigel und die zusatzlichen 
von der Pfeilung herrithrenden Auftriebsénderungen ermittelt. 

Die Auftriebsverteilungen bei / — 0 und die hiermit gewonnenen Neutralpunktlagen werden 
mit theoretischen Rechnungen verglichen. Es ergibt das Verfahren von Weissinger eine recht 
gute Ubereinstimmung mit der Messung. Die Unterschiede zwischen den gemessenen und nach 
diesem Verfahren gerechneten Neutralpunktlagen sind etwa 1% der Fliigeltiefe. 

Fitr die Berechnung der Auftriebsverteilung an einem schiebenden Pfeilfliigel wird ein Nahe- | 
aT oeemete oe angegeben, das beim Vergleich mit der Messung recht brauchbare Ergebnisse 
liefert. . | 


| 
| 
(Eingegangen am 16. August 1950.) | 


Anschrift des Verfassers: Dr. habil. W. Jacobs, Stockholm (Schweden), Flygtekniska férséksanstalten, 
Ulvsunda. 
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Kin Integrationsverfahren fiir die Berechnung der Biegespannungen 
achsensymmetrischer Schalen unter achsensymmetrischer Belastung*. 


I. Mitteilung. 
Von H. Miinz. 


1. Einleitung. In der Reifner-Meifnerschen Theorie der achsensymmetrischen Schalen unter 
achsensymmetrischer Belastung erfordert die Berechnung der Biegespannungen die Kenntnis 
der allgemeinen Lisung einer gewodhnlichen linearen Differentialgleichung vierter Ordnung. 
H. Reiner’ hat erstmalig bei der Berechnung der Kugelschale konstanter Wandstirke erkannt, 
da®B sich die schwer zu iiberblickende Fiille von Gleichungen dieses Elastizitatsproblems itber- 
sichtlich zusammenfassen 1aBt in zwei gekoppelte, lineare Differentialgleichungen zweiter Ord- 
nung fiir die Querkraft und die Verdrehung der Meridiantangente. E. Meifner? gelang es dann, 
auch fiir Schalen beliebiger Meridiankriimmung und veranderlicher Schalendicke entsprechende 
Differentialgleichungen aufzustellen, die in Sonderfallen durch eine komplexe Differential- 
gleichung zweiter Ordnung ersetzt werden kénnen. F’. Télke® hat spater gezeigt, daB bei gewissen, 
mit der Theorie vertraglichen Vernachlassigungen, das Integrationsproblem immer auf die Lésung 
einer komplexen Differentialgleichung zweiter Ordnung hingefiihrt werden kann. 

Im allgemeinen Fall miissen diese Differentialgleichungen natiirlich schrittweise numerisch 
integriert werden. Fiir die Reifner-Meifnerschen Differentialgleichungen hat W. Lohmann* 
ein solches Integrationsverfahren entwickelt. 

Eine Schalenrechnung erfordert gewoéhnlich einen sehr erheblichen Rechenaufwand. Im 

Bestreben, diesen Aufwand auf ein Minimum herabzudriicken, wird man sich die Frage stellen, 
ob sich die von Reifner-Meifner und Télke aufgestellten Differentialgleichungen als die fiir 
eine numerische Rechnung am besten geeignete Formulierung des Schalenproblems darbieten. 
Die Koeffizienten dieser Differentialgleichungen sind sehr verwickelt gebaut. enthalten erste 
und zweite Ableitungen der Meridianform und der Schalendicke und sind, wie die Differential- 
sleichungen selbst, durch eine Reihe von Differentiations- und Eliminationsprozessen entstanden, 
deren mechanischer Sinn wohl kaum noch durchschaubar sein diirfte. 
_ Der Versuch, ein der Rechnung besser angepaBtes Gleichungssystem aufzufinden, gab den 
Anstof& zu dieser Arbeit. Der Leitfaden, der uns hierbei durch die Reifner-Meifnersche Theorie 
fiithren wird, ist die Variationsrechnung und die Integrationstheorie der kanonischen Differential- 
gleichungen, wie sie in der analytischen Mechanik von Lagrange, Jacobi und Poincaré entwickelt 
worden ist. Die beste Orientierung itber diese Methoden hietet der Enzyclopadieartikel von 
G. Prange®. Prange selbst hat in seiner Habilitationsschrift® die Theorie der kanonischen Systeme 
auf das Fachwerk und den elastischen Kérper angewendet und in einer besonderen Arbeit’ 
auf den Balken. Prange zeigt die engen Analogien auf, die zwischen den Variationsproblemen 
der analytischen Mechanik und denen der Festigkeitslehre bestehen. Neuerdings hat W. Giinther* 
diese Verfahren der Variationsrechnung fiir die Kreisringplatte nutzbar gemacht. 


2. Die Grundgleichungen der Schalentheorie. Wir wollen vorweg die Grundgleichungen der 
Reifner-Meifnerschen Schalentheorie zusammenstellen. Dies kann in aller Kiirze geschehen, 


* Dissertation Braunschweig 1950 (Hauptberichter: Prof. Dr.-Ing. H. Schaefer, Mitberichter: Prof. 
Dr.-Ing. O. F6ppl). 

1 H. Reifner, Miiller-Breslau-Festschrift 1912 S. 181, Leipzig 1912. 

2 EF. Meifner, Physik. Z. 14 (1913) 5. 343. 

3 F, Télke, Ing.-Arch. 9 (1938) S. 282. 

4 W. Lohmann, Ing.-Arch. 6 (1935) S. 338. Vgl. auch R. Grammel, Ing.-Arch. 10 (1939) 5. 407. 

5 G. Prange, Enzyclopadie der math. Wissensch., IV, 2; Snp0 ps 

6 G. Prange, Das Extremum der Formanderungsarbeit, Habilitationsschrift Hannover 1916, unver- 
offentlicht. ¥, 

7G. Prange, Die Theorie des Balkens in der technischen Elastizitatslehre. Z. Arch. Ing.-Wesen 65, 
N. F. 24 (1919) S. 83—96, 121—150. 

8 W. Giinther, Die Biegung kreissymmetrischer Ringplatten veranderlicher Dicke als Problem der 
Variationsrechnung. Dissertation Braunschweig 1946. 
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da ausfiihrlichere und sehr itbersichtliche Darstellungen in den bekannten Biichern von W. Fliigge* 


und Biezeno-Grammel? zu finden sind. 
Das Element einer achsensymmetrischen Schale ist begrenzt durch zwei benachbarte Meri- 


dianschnitte und zwei benachbarte Breitenkreisschnitte (Abb. 1 und 2). 


geometrischen Daten und der Beziehungen zwischen den Verformungen und Verschiebungen aus-; 
| 


Q. Ardimmungsmitteypunkt & 
\\ . des Meridianelements NS | 
~ 1 


EX pow 
ae ag-~ \AN 


No ONS Ss VES 
NG ae acs mS 
\ \ WN x Hy Se 
\ See — Ludvede be [pt 4(rbp)/a3 
¢ Me ae Mittelflche ee \ bs [Myra Mp/fad 
\ y 


x 
; ao 
ra aly Ng 
we 


<P Scholenachse \\i ff s f 7My Be “A is / 
Mittelpunkt aes Breitenkreises — Lweiten Krdimmungs- rgd eS ss aw WN »ap 
mittelounkt des ? Mi 


Schalenelements 


Abb. 1. Geometrische Daten. 


. pm = Winkel zwischen Meridiantangente und Schalen- 
achse. 
r, = Kriimmungsradius des Meridians. 
r,dp = Bogenelement des Meridians. 
r = Radius eines Breitenkreises. 
rd = Bogenelement des Breitenkreises. 
x, y,z — Achsen eines rechtshindigen, kartesischen Koordi- 


a 
ay ee 


\] 
me _Schalenachse, 


Abb. 2. Schnittkrafte und -momente. 
No = Langskraft (in Richtung der Meridiantangente). 


Q = Querkraft (in Richtung der Normalen der Schalenmittei- 


flache). 
No» = Ringkraft (in Richtung der Breitenkreistangente). 
My = Biegemoment (um die Breitenkreistangente). 
My = Ringmoment (um die Meridiantangente). 


natensystems, und zwar 

x-Achse in Richtung der Breitenkreistangente, 

y-Achsein Richtung der Meridiantangente, 

z-Achse in Richtung der Normalen der Schalen- 
mittelfliche. 


reicht. Die Schnittkrafte und -momente werden 
ohnehin auf den Mittelpunkt jeder Schnittflache 
bezogen. 

Es seien o, und oy die Normalspannungen (Zug positiv, Druck negativ), die am Breitenkreis- 
schnitt bzw. Meridianschnitt angreifen. Wegen der Achsensymmetrie von Schale und Belastung | 
treten nur am Breitenkreisschnitt Schubspannungen t auf. Die Spannungen o, normal zur Mittel- 
flache seien im Vergleich zu den anderen Spannungen so klein, da® wir sie und ihren EinfluB 
auf die Formanderungen vernachlassigen kénnen. Die aus den Spannungen resultierenden 
Schnittkrifte und -momente werden je Langeneinheit des zugehérigen Bogenelementes der 


Schalenmittelflache angegeben. Es folgt, wenn z= -- h die duBere bzw. innere Schalenlaibung | 
bedeutet, 
+h | 
z Cos 
“= [o,(1+ 4 ) dz, an 
aay 
+h 


Th 


Vp = [ o»(1 ane —) 2dz, 


fe 
f(t E22) as, (3) 


=i 


Os = 


* W. Fliigge, Statik und Dynamik der Schalen. Berlin 1934. 
* Biezeno-Grammel, Technische Dynamik, S. 450. Berlin 1939. 
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gts | a tee dz, (4) 
ve 
| ao(1—) edz. (5) 


Die Kriimmungshalbmesser rund r, sind nach Voraussetzung sehr viel gréBer als die Schalen- 
dicke 2h, so dab die Anderung von r, und r mit z zwischen den Schalenlaibungen vernachlassig- 
bar ist. In die Integrale (1) bis (5) werden daher fiir r, und r die Werte der Schalenmittelflache 
eingesetzt. Fiir die Normalspannungen nehmen wir eine lineare 


o (2) =o +2 x. (6) 
Aus (1), (2), (4), (5) folgt mit (6) 


h cos 
WN, = 2h (on +—— Pea h cos 
ro ot Seon) Mieke cs a lala ae 0): 
| h 
Ne = 2h(os—+— 0 Cmegisa x 
0 3 Ts O1] ? Ms == a oy h O91. = = O90] - Abb. 3. Verteilung der Normal- 
1 Spannung iiber die Schalendicke. 


Aus diesen Gleichungen werden die an der auBeren und inneren Schalenlaibung auftretenden 
) Normalspannungen berechnet. Da h <r, und r ist, wird 


1 Se 
Opi h— oF (Ne 3 m) (7) 
1 meV 
CoLhk ar (No See | . (8) 
_ Die Gleichgewichtsbedingungen am Schalenelement lauten (Abb. 2) 
d 
—(r N,) —rQ,—17, Nesing=—rr,Y, (9) aN 
dy \ 
d \ 
rN. —— (rQ,) — r, Necos y= —rr, Z, (10 tee 
: e+ dy (7 Qe) et 1 Z, (10) aN 
gt aS? hr Up 0. (11) mits ay 


Dabei sind Y und Z die auf das Schalenelement 
einwirkenden dufBeren Krafte in Richtung der 
Meridiantangente bzw. der Normalen der Schalen- 
mittelflache (Dimension Kraft/Flaiche). 

_ Fir die vorliegende Theorie werden die auf- 
tretenden Verformungen und Verschiebungen als 
infinitesimal angenommen. Der Formanderungs- 
austand diimner Schalen wird durch die Ver- 


formung der Schalenmittelflache angegeben. on r| | r(teEy) 
‘AuBerdem sei vorausgesetzt, daB die Normale ~~~ __ “ _|__ Schalenachse 
der Schalenmittelflache auch nach der Verfor- Riving Reh on eee ee 


mung noch normal zur neuen Schalenmittelf lache — Verschiebung i Richtung der Meridiantangente. 
‘ist. Diese Bedingung ist nur dann erfillt, wenn — Verschiebung in Richtung der Normalen der Schalenmittel- 
Z flache. 
wir keine Querkraftverformungen zulassen. 1% = Verdrebung der Meridiantangente. 
‘ p ' ds =r, dp = Bogenelement des Meridians, 
Dieses setzt wiederum nur die Betrachtung oe 
; 5 . eg = Meridiandehnung. 
idiinner Schalen voraus. Die Verschiebungen v_ ,, — Ringdehnung. 
und w (Abb. 4) sind positiv in Richtung der 
jpositiven y, z-Achsen; die Tangentenverdrehung y ist positiv, wenn sie, in Richtung der 
positiven x-Achse gesehen, im Uhrzeigersinn erfolgt. 
' Die Beziehungen zwischen den Verformungen und Verschiebungen, die an Hand von Abb. 4 
armittelt werden kénnten, sollen hier mit Hilfe der virtuellen Arbeiten der inneren und auweren 


ec 


q 
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den virtuellen Verriickungen A, = dv, A, = dw, A, = 6x multipliziert, dann addiert und itber eine: 
Schalenzone integriert. Wir erhalten so fiir die gesamte virtuelle Arbeit der duBeren und innerent 


Krafte die Gleichung 


TIN) ; d 
x [A ie -rQ,-4 ri Nosing +7, | rN, oe 4 rs Ny cos | + 


| As| eM) 1 ri Mo sin p—r7sQy| ip taf ruta dpe=0e 


Werden die Glieder mit der ersten Ableitung nach ¢y einmal partiell integriert, so erscheint dies 
virtuelle Arbeit der inneren Krafte in der Form 


2a | InN, (Te 22) + No Ay sing 10089) +109 (A + a ns) | | 


7M, S++ Mo lysing dp. | 


Die ausintegrierten Glieder geben die virtuelle Arbeit der Randkrafte an. 

Wenn de, und dey die virtuellen Verzerrungen des Meridianelementes bzw. des Breitenkreml 
elementes sind und 6x, und 6%, die entsprechenden virtuellen Verkriimmungen, so ist die virtuelley 
Arbeit der inneren Krafte ohne Querkraftverformung 


27 [ {No 6€)+ Ny de,+ M, 0x, + M,dx»} rr,dq. (13) 


Vergleichen wir die beiden Formen der virtuellen Arbeit der inneren Krafte (12) und (13) miteiml) 
ander, so ergeben sich zwischen den virtuellen Verformungen und den virtuellen Verschiebungen} 
Beziehungen, die auch fiir die wirklich auftretenden infinitesimalen Verformungen und Verschiag 
bungen gelten: 


=. 


(12 


dv w ; 
Ted a ee (14) 
jp (15) 

dy 
Kop rd (16) 
%s= (17) 


dw v dw 
Mecha aie are ho io 
d.h. der Verschiebungszustand ist bereits durch zwei der drei GréBen v, w, y bestimmt. 
Zu den Gleichgewichtsbedingungen (9) bis (11) und den hier abgeleiteten Beziehungen (14)) 
bis (18) tritt fiir die 12 unbekannten GréBen No, No, Qo; Mo, Mo; v, w, 13 Eq, E93 Ho» % NOCH deat 
Zusammenhang zwischen den SchnittgréBen und den Wertienitnigen! Fiir homogene, isotrope } 
Schalen bei infinitesimalen Verformungen gilt das Hookesche Elastizitatsgesetz. Daraus ergeben 
sich folgende Beziehungen, falls die Abhangigkeit der Verformungen vom Abstand von der 
Schalenmittelflache vernachlassigbar ist: 


2uN, = D(e, +689), (19) 
24%Ny = D(és + v¢,) , (20) | 
20M, = K(x, + m5) ; (21) ) 
2% Ms = K(x» + %,), (22) : 


worin D = 27 a die Dehnungssteifigkeit ist und K — 20k . een die Biegesteifigkeit | 
v 
mit dem cee E, der Querkontraktionsziffer » und ae ania Schalendicke h. 


3. Das Prinzip vom Minimum der potentiellen Energie. Die potentielle Energie der belasteten | 
Schale ist | 


[hes 27 ie [Ny €o-+ No eo-+ Mg %o-+ Myo x9] rr, dy — [ e420) nT dp} : 
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wenn wir zur Vereinfachung der Schreibweise von Randkraften und -momenten absehen (die 
Grenzen der Integrale sind durch die Begrenzung der Schalen durch zwei Breitenkreise gegeben). 
Hierin sind mit Hilfe des Hookeschen Gesetzes die Schnittkrifte und -momente durch die Ver- 
formungen auszudriicken und diese aus den Verformungs-Verschiebungs-Gleichungen durch die 
Verschiebungen, so daB die potentielle Energie schlieflich in folgender Form erscheint : 


ji 1D dv ee easy dy w\vsinp+ weosp , /vsing + w cos p\2 
2 aa KD py r,dy iy r r 


dy \2 dy ysing x sin ~\2 
| KI, a) ms ats : H( : J] rnd — 22. f (Yo Ze) rr dp. 


Das Prinzip vom Minimum der potentiellen Energie fordert 

OL ==... 

eater den méglichen Verschiebungszustanden ist derjenige herauszugreifen, der bei Beachtung 
ider geometrischen Randbedingungen die potentielle Energie zum Minimum macht. Nun besteht 
‘zwischen den Verschiebungen v, w, 7 die Gleichung (2, 18), die wir als Nebenbedingung mit einem 
|Lagrangeschen Faktor A versehen dem Integranden von // hinzufiigen wollen, um die drei Ver- 
sschiebungen unabhangig voneinander variieren zu kénnen. Die potentielle Energie erhalten wir 
‘somit in der Gestalt 


P dv dw d 
T= [E (rms da’ Ip? ap) PA) ap 


= : ie & Ww) 20 (Tp w) (v sin p+ wos g) + 3 (v sin p+ w cos a)? |+ | 
r [dy dy 7 - 
+kK Fe(ae) +27gezsine +7 sing) 420A 2 i —1y)—20 (Yv+Zw)rr, | 


oie 


(1) 


Die Eulerschen Differentialgleichungen unseres Variationsproblems, die man beziiglich v, w und ¥ 
zu bilden hat, sind die Gleichgewichtsbedingungen, ausgedriickt in diesen VerschiebungsgréBen. 
Stellt man diese Gleichungen auf, so zeigt ein Vergleich mit (2, 2) und (2, 9) bis (2, 11), daB der 
\Lagrangesche Faktor die Bedeutung rQ,, hat. Dies war zu erwarten, da A die der geometrischen 
‘Nebenbedingung entsprechende Reaktionskraft ist. Diese Nebenbedingung (2, 18) folgte aber 
jaus der Voraussetzung, dafs eine Schubverformung ausgeschlossen ist, wodurch aber die Quer- 
tkraft zu einer Reaktionskraft entartet. 

Da diese drei Eulerschen Gleichungen im folgenden nicht benétigt werden, verzichten wir 
fauf ihre Wiedergabe, bemerken jedoch, da man aus ihnen durch Elimination von A und x bei 
‘Benutzung der Nebenbedingungen (2, 18) zwei Differentialgleichungen dritter Ordnung fiir v 
jund w von allerdings sehr uniibersichtlichem Aufbau erhalten kann. 


_ 4, Das kanonische Variationsproblem und die kanonischen Gleichungen. Nach dem Beispiel 
ler analytischen Mechanik ersetzen wir nun in der ,,Lagrangeschen Funktion“ L (3, 1) die ersten 


‘Ableitungen der Verschiebungen nach der unabhiangigen Verdnderlichen p durch die zugehérigen 


oL oL oL 
eee aac Pa a i? Pi Gy 2 
p= DIF (v' — w)+-v (using +-weos 9)| = 2nrQN,, (1a) 
} Pw=20A Jar), (1b) 
pa= K [Fo +ozsing) ==2arM, . (1e) 


iDie Ableitung nach g wird durch den Strich bezeichnet.) Durch eine Legendresche Transfor- 
mation wird die Hamiltonsche Funktion eingefihrt: 


UD (Dope Pes © Us) = Pol + p,w' + p, x4 — Liv’, w', x5 v, w, 43 An) 


ind das Variationsproblem vom Minimum der potentiellen Energie geht iitber in das kanonische 


i 


1 Diejenigen Glieder von L, die in v’, w’, x’ linear sind, liefern keinen Beitrag zu H; w’ tritt iiber- 
q . . 
faupt nur linear in L auf. A 
8 
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| 
Variationsproblem 
bf {p,v’ pe py — iH, dg 0), (2a)} 


Jetzt werden die ,,Impulse“* und Verschiebungen, die ,,Koordinaten” im Sinne der analytischen| 
Mechanik, unabhangig voneinander variiert, und es ergeben sich daraus die kanonischen Glei-| 


chungen 
he oH be Ort 
Pee Pe ye Ps 
; oH F OH 
enh) PaaS 
ye OH yy — OH 
Pe Gr Me ay 
Die Hamiltonsche Funktion ist 
A= — pov + pypW +1 PwX + 
a aly: a 2 y Pe (v sin y + woos g) — (1 — »”) (vsin p + weos @)?|4- 
Der AND D (2) 


+ Kz) 27K xsing— (1 — v2) 72 sin? p 4 2nrr,(Yo+Zw), | 


und demnach lauten die kanonischen Gleichungen 


Pvp sin y = Pot D(1 — »?) 2 (v sin p+ w cos 7) sin p — 2a r7r,Y, (3a) | 
Po—? 2p, cosp=— p,+D(1- vw) (v sin p + w cos g) cos p— 2ar7,Z, (3) | 
Py —? 2D, sing = — r, py + K (1 P) 2 x sin? @ , (3c) 
v + y2ovsing = w (1—» cos g)+ 22 ; (3a)| 

w=—v+trny, (3e} 
L+y2ysing=ie. (3f) 


Die kanonischen Gleichungen bilden cin System von sechs linearen Differentialgleichungen erster 
Ordnung fiir die Verschiebungen v, w, y und die drei Impulse, die nach (la) bis (1c) die an einem 
Breitenkreise angreifende Langskraft und Querkraft angeben und das zugehérige Biegemoment. 

Die erste Gruppe der kanonischen Gleichungen (3a) bis (3c) wiirde man direkt aus den Gleich- | 
gewichtsbedingungen (2, 9) bis (2, 11) erhalten, wenn man die Ringkraft Ny und das Ringmoment 
My, durch die Verschiebungen ausdriickt. Es ist naémlich 


vsin y + w cos p 
9 
r 


2aNos=v2a0N,+D (1—??) 


2 My =v 2.0M,+ K (1—p2) 2288 


Die Gleichungen (3d) und (3f) der zweiten Gruppe der kanonischen Gleichungen ergeben die 
Definitionen der Impulse p, und p, nach (1a) und (1c) — das sind also die Beziehungen zwischen 
Langskraft N, und den Verschiebungen v und w bzw. zwischen dem Moment M, und der Ver- 
drehung 7. Die kanonische Gleichung (3) ist die Nebenbedingung zwischen den Verschiebungen 
(2, 18). 


5. Reduktion des kanonischen Gleichungssystems mit Hilfe eines bekannten Integrals. Unter 
einem Integral des kanonischen Gleichungssystems (4, 3) versteht man einen Ausdruck von der 
Form G (p,, Pw. Py; V, W, 3 P), der fiir jede Lésung p, (g),... v(p),..- des kanonischen Systems 
unabhangig von ¢, also konstant wird. Wie in der Hamilton-Jacobischen Integrationstheorie der 
analytischen Mechanik gezeigt wird, lassen sich aus der Kenntnis eines Integrals Vorteile fiir 
die Integration eines kanonischen Systems gewinnen. Ware z. B. ein Integral bekannt, das linear 
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in den Impulsen ist, also die Gestalt 
a4 (2, 1, 13 9) Pot a (0, 1, 13 P) P+ (0, 5 13 9) Py + a4(0, 10, 43 9) = konst. 


hat, so lieBe sich unser kanonisches System von sechs linearen Differentialgleichungen erster Ord- 
nung auf ein solches mit nur vier Differentialgleichungen reduzieren!, 

Ein Integral des vorliegenden kanonischen Systems laBt sich aber tatsichlich angeben. Wie 
man leicht sieht, erhalt man bei der Elimination von N, aus den Gleiéhapwiohtobedtupudzen (2, 9) 
und (2,10) die sofort integrierbare Gleichung 


a ie Ne cos gy — rQ, sin g] = — rr, (Y cosy — Z sin P) ; (1) 


da auf der linken Seite die Ableitung der Komponente der resultierenden Schnittkraft in Rich- 
tung der Schalenachse steht und auf der rechten die Komponente der resultierenden Belastung 
in Richtung der Schalenachse. Daf dieses Integral angegeben werden kann, liegt darin begriin- 
det, dafs bei einer Verformung der Schale die Verschiebung der Schalenpunkte in axialer Rich- 
tung immer von einer beliebigen Verschiebung der starren Schale in axialer Richtung iiberlagert 
werden darf. 


Die erste Gruppe der kanonischen Gleichungen (4,3) stellen die etwas umgeschriebenen 
Gleichgewichtshedingungen dar. Es muf sich deshalb die der Gleichung (1) analoge Beziehung 
n den Impulsen ergeben, wenn (4, 3a) mit cos m und (4, 3b) mit sin y multipliziert werden und 
man beide voneinander abzieht: 


d 
ae (P, Cos YP — Py sin p) = — 2arr, (Y cosy — Zsing). (la) 


Wir fiihren als neuen Impuls die bei gegebener Belastung bekannte Funktion 


Pi = Pv COS P — Pw Sin P (2) 
pin. 

Damit man sich bei der Einfithrung neuer Impulse und Koordinaten nicht dér Integrations- 
vorteile eines kanonischen Gleichungssystems begibt, mu das urspriingliche Variationsproblem 
4, 2a), wieder in ein kanonisches Variationsproblem iibergehen; die Impulse p,, py, p, und die 
Xoordinaten v, w, ¥ miissen also durch eine ,,kanonische Transformation‘ in die neuen Impulse 
1, Po, p, und die zugehorigen Koordinaten q,, qo, q3 itberfiihrt werden. Die Bedingung hierfitr ist, 
laB sich die Integranden des alten und des neuen Variationsproblems héchstens um ein totales 
Differential unterscheiden diirfen, da ein solches Zusatzglied die Euler-Lagrangeschen Gleichungen 
nicht beeinfluBt. Die ,,erzeugende Funktion“* V der kanonischen Transformation ist hierbei eine 
unachst beliebige Funktion der alten Impulse und der neuen Koordinaten, V— V(p,. Pw. P, 
hy I> 33 gp), fiir die wir folgenden Zusammenhang mit den neuen Impulsen und den alten 
<oordinaten vorschreiben: 


_av Wek EN 
Ly ae FR Px Fa,’ Pa qa | ’ 
este eer)” i | (3) 

Col ho OP, OPy 


Der Zusammenhang der neuen Hamiltonschen Funktion H, (p,, P2: P33 Vis 2» 93) mit H (4, 2) 
iwerde gegeben durch 


OV 
= at 4, 
HH (4) 


Der Integrand des neuen kanonischen Variationsproblems muf wieder die Gestalt haben 

} Pit +P2% + P3943 — Ay. 

nter Beriicksichtigung der Transformationsgleichungen (3) und (4) laBt sich folgende Umrech- 
mung angeben: 


2 / / d if / / 
Pit, + Pah + Pd — Hy = G5 LV — (pe + wPw +X P2)| + Po? + Po + Pit — Ht. 


/ 1 Zu diesem in der Hamilton-Jacobischen Integrationstheorie bekannten Satz vgl. G. Prange, Encyclo- 
hadie, S. 735. 
9x 
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Die Integranden des alten und des neuen kanonischen Variationsproblems unterscheiden sich 
also wie verlangt um ein totales Differential. 


Da im vorliegenden Fall nach Gl. (2) und (3) 


aV : 
ih = Pp, COS Y — Py sin P 
sein mu}, wahlen wir die ,,erzeugende Funktion*‘ V zu | 
V = (py cos p— pw sin p) H+ (Pefi(P) + Pua (P)) 42+ Pr I (5). 
und erhalten aus den allgemeinen Gleichungen (3) der kanonischen Transformation 
Pi = Pv.00s PW — Py Sin” , (6a) | 
P2= Po fi(%) + Pw fo (Q) » (6b) ] 
P3= Px» (6c) 
Ube 1 OS P+ fi (¥) > (7a) 
w=—qsing+® fog)» (7b) | 
io EG (Ze) | 


Hierbei sollen f,; und f, zwei zunachst beliebige Funktionen von y sein. Aus (7) ergeben sich die? 
neuen Koordinaten zu 


oe v f,— wf, 
Lime cosp+f,sing’ (Bay I 
__ vsing + w cos 9 | 
a cosp+f, sing’ (8b) 
INS Gian (8¢e) 
Die neue Hamiltonsche Funktion ist 
OV vy Pw, > > Woo > 
H, = H (py, Par Py %, % 139) (Pe. P Me i> I2> Iss P) (9) | 
mit 
OV : ij j 
Bp a (Po Sin P+ Pw 008 Y) H-+( Poi (P) + Pos (P)) (9a) 


und H nach (4, 2). In (9) werden die urspriinglichen Koordinaten v, w, y nach (7) durch die neuen 
di > 42> 93 ersetzt und die urspriinglichen Impulse p,, p,, p, durch die neuen p,, p,, p3, indem 
man (6) nach p,, py, p, auflést: . 


_ _Pifa+ p,sin p ‘ 
Po Ff cosp+fsin®’ (10a) | 
__ _P, cos p— pif 
Pw Ff, cos p+ f, sin y’ (10b) 
Py ae Pea (10c) | 


Der hier gemachte Ansatz (5) fiir die erzeugende Funktion V, fiir die lediglich die Gleichungen (2) 
und (3) vorgeschrieben sind, ist nicht der allgemeinste Ansatz, doch ist er wegen der vorhandenen | 
Symmetrie zwischen den transformierten GréBen sehr augenfallig. Offensichtlich wird die Trans- | 
formation einfach, wenn die Koeffizientendeterminante | 


fx(p) cose +fi(y) snp=1 (11) jj 
gesetzt wird. Aus dieser Bedingung folgt 
92 = VsiNy-+weosg, 
d. h. q,ist unter der Bedingung (11) immer die aus v und w resultierende Verschiebung senkrecht | 


zur Schalenachse. Um (11) zu erfiillen, sind drei Kombinationen von f,(p~) und fy (~) sehr ein- || 
fach zu iibersehen, auf die wir uns beschranken wollen: . 


(a) fi) =sing, f2(p) = cosy, 
1 


(b) A@= fry) = 0, 


sing’ 


(c) Aly)=9, fa(¢) = 


i 
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Zu (a): Nach (6) und (8) wird 


Pi = Pr COS P — p, sing , q = veosp— wsing, | (13a) 
P2= Pv SIN Y + Pw cosy, q92=vsing+weosy. | 


Unter Beriicksichtigung der Bedeutung von p, und p, erkennt man, da durch diese Transfor- 
‘mation an Stelle der Langskraft N, und der Querkraft Q, die Resultierenden dieser Schnittkrafte 
§ parallel und senkrecht zur Schalenachse treten, und entsprechend an Stelle der Verschiebungen v 
jund w deren Resultierende parallel und senkrecht zur Schalenachse (Abb. 5a und 5b) 


Zu (b): Aus (6) und (8) folgt 


° Ww 
Pi => Pv ©OS GY — Py sin ®@ , Ui sing 2 | 


(13b) 


93 = vsing -+weosp. | 


wa 
a ee Ze 
i 
y “9\ y a\ . 2 ; 
AN Lad ; Ze A Jernachse 
fo 7 Schalenachse__ * Z Schalenachse __ ee 7 _ Schalen oe 
Abb. 5a. Fiir Fall a) Zusammenhang Abb. 6a, b. Abb. 7a, b. 
zwischen q,, qg, und v,w. Abb. 5b. Fiir Fall b) Deutung von p, und q,. Fiir Fall c) Deutung von p, und q,- 


* Zusammenhang zwischen p,, p, und Pro Pry 


In diesem Fall ist p, im wesentlichen die Langskraft N,,, doch kann man nach Abb. 6a p, deuten 
jalseine Kraft senkrecht zur Schalenachse, deren Komponente in Richtung der Meridiantangente 
idie Langskraft ist. Ebenso ist die Koordinate q, zu deuten als eine Verschiebung parallel zur 
'Schalenachse, deren Komponente in Richtung der Normalen der Schalenmittelflache die 
fnegative Verschiebung w ergibt (Abb. 6b). 


Zu (c): Aus (6) und (8) folgt 
Pi = Pv COSY — Py sing, q = : | 


~ cos p’ 
i sa (13) 
cos Pp” 


q,—= vsing+weosg. 


Poa 


Wegen Da 227Q, ist in diesem Fall p, eine Kraft senkrecht zur Schalenachse, deren Kompo- 
‘nente in Richtung der Normalen der Schalenmittelflaiche die Querkraft ergibt. Diese Variabele 
‘und die Verdrehung y haben Reifner und Meifner! zur Aufstellung ihrer beiden Differential- 
'gleichungen zweiter Ordnung benutzt. Die Koordinate q, ist in diesem Fall eine Verschiebung 
parallel zur Schalenachse, deren Komponente in Richtung der Meridiantangente die Verschie- 
bung v ist (Abb. 7a und 7b). 

; Die Hamiltonsche Funktion fiir die neuen Impulse und Koordinaten lautet, falls wir die bis 
auf die Bedingung (11) beliebigen Funktionen f, und f, beibehalten, 


| 1 H, Reifner, a. a. O., E. Meifner, a. a. O. und E. Meifner, Vjschr. Naturforsch.-Ges, Ziirich, 60 (1951), 
S23. 


mh . oe . 1 Moh hivl 
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A,= q.Pi (f2-+-F3-44 Sife fife) 114 (p, cos p — Py fi) + 


$5 [p[(efat patina? _ 9 yg, Pale tae? — (1 — st) g§| + | 
a7 | D ' D (14)) 


+ K|(e) — 2 vg, sin 9 — (1 v9 sing + 


+ 2arry [q (Y cos p — Zsin gy) +o (Vf +Zfa)| - 
Die Koordinate q, tritt nur linear als Koeffizient der Belastungskomponente in Richtung der) 
Schalenachse auf. Die rechte Seite der kanonischen Gleichungen kann also keine Glieder mit 4; | 
enthalten. Aus den kanonischen Gleichungen erhalten wir daher die Liésung 4 (¥) durch eine | 
Quadratur, falls die Lisungen der restlichen Koordinaten und Impulse aus den linearen Diffe-. 


; ; / 0H, 
rentialgleichungen erster Ordnung ermittelt worden sind. Die kanonische Gleichung q; = ap, | 
wird hier 

/ BY 2 / / r rT, Pi fo + P2 sin Y 
=o (FR fe fe fale) — 9 ee i (15) 


Man bezeichnet q, als eine ,,zyklische Koordinate‘‘. Wir haben also an Stelle des linearen Glei- | 
chungssystems von sechs Differentilagleichungen erster Ordnung (4,3) in den kanonischen 
Gleichungen, die zu der obigen Hamiltonschen Funktion gehéren, ein lineares Gleichungssystem | 
von nur vier Differentialgleichungen erster Ordnung erhalten fiir p., p; und q,, q3. Zu diesen treten 
noch die beiden Quadraturen fiir p, und q, hinzu nach (la), (2) und (15). | 

Da in der Hamiltonschen Funktion (14) die beliebigen Funktionen f, und f, nur als Koeffizien- 
ten von p, und der Belastungskomponente parallel zur Schalenachse auftreten, stimmen die homo- 
genen Gleichungen fiir p,, p3; und q,, gz in allen drei Fallen (a), (b) und (c) itberein. Wie wir im fol- 
genden sehen werden, besteht die Hauptarbeit in der Integration der homogenen Gleichungen; 
daher wiire es zunachst gleichgiiltig, fiir welchen der drei Falle wir uns entscheiden. Wenn wir den | 
Fall (a) vorziehen werden, dann deshalb, weil p, und p, statisch anschaulich zu deuten sind, nam- | 
lich als die Komponenten der resultierenden Schnittkraft parallel und senkrecht zur Schalen- 
achse, und weil q, und q, als Verschiebungen parallel und senkrecht zur Schalenachse unmittel<_ 
bar praktisch interessieren. 

Die homogenen kanonischen Gleichungen lauten in allen drei Fallen 


2° (16a) 


pi—v" py sing = D(1—»*)* 


1 
= 


P3— yp, sin p = K(1— 2”) . qs; sin? ~ — 1, pp cos @ , (16b) 


tc Ps 
Rp 


q2 +94 q2 Sin Y = sin? p +-1,q3 cos, (16c) | 


D 

g +y2gsingp= 28, (Lodi 
denn fitr p; = 0 werden die Querkraft und die Lingskraft die Komponenten der resultierenden _ 
Schnittkraft senkrecht zur Schalenachse. Falls also die Schale nur durch Randkrafte senkrecht _ 
zur Schalenachse und durch Randmomente belastet ist, | 
hat p, in allen drei Fallen dieselbe statische Bedeutung — 
(Abb. 8). 

Wenn man nun das Moment p, und die Verschiebung 
senkrecht zur Schalenachse q, aus diesen vier Differential- 
gleichungen erster Ordnung eliminiert, erhalt man fiir die 
Verdrehung q, und die Schnittkraft senkrecht zur Schalen- 
achse p, das lineare Differentialgleichungssystem zweiter 

4 % Ordnung, das Reiner und Meifner aufgestellt haben. 
= \ -——1 = aC Und ee liefert i Fall (c) ae ee Gleichungen, 
wahrend die beiden anderen Faille, wie schon erwahnt, 
andere inhomogene Glieder aufweisen. Der symmetrische 
Aufbau des Gleichungssystems (16) 1a8t sofort erkennen, daB es leichter ist, mit (16a) und (16d) 
p; und q, zu eliminieren, als q,; und p, mit (16b) und (16c), und so wird man zwangslaufig auf 
die homogenen Reifner-Meifnerschen Gleichungen fiir die Verdrehung q; und die Schnittkraft 


Abb 8: Zusammenhang zwischen p, und Qy> Ny 
fiir die unbelastete Schale. 
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_ senkrecht zur Schalenachse p, gefiihrt, die z. B. fiir konstante Wandstirke (D und K konstant) 
} lauten 


(F- pt) — pa(“ sin? g +r cos p) = D (Lv) r gy cos p., 
1 

ie ae Tis r 

(7:08) — 49("2 sin? p— 208 9} = — nr, Reeos py 


6. Die Grundgleichungen und die kanonischen Gleichungen, aufgestellt in den achsenfesten Ko- 

f ordinaten. Wenn nun einmal bekannt ist, daB die Komponenten senkrecht und parallel zur Schalen 

| achse der resultierenden Schnittkraft und die entsprechenden Verschiebungen die zu bevor- 

' zugenden Variabelen der Schalentheorie sind, so kann man gleich von einem neuen Koordinaten- 

| system ausgehen, dessen Achsen in der Schalenachse und 

) senkrecht hierzu liegen. Wir fiihren also folgendes &, 1, ¢- 

/ Koordinatensystem (Abb. 9) ein: U) 

| &-Achse in Richtung der Breitenkreistangente (negative 
x-Achse), 

y-Achse senkrecht zur Schalenachse, 

¢€-Achse in der Schalenachse. 


( Der Ursprung dieses Koordinatensystems wird fest auf der 
( Schalenachse gewahlt. Die Meridiankurve der Schale ist 
|} dann gegeben durch r=/f(¢). Das Bogenelement des “9 

a6 : Trey -: COUY. | Sohalenachse G 

) Meridians ist d s = jl +f’? d¢1 und der Kriimmungshalb- : ee ae 2S 
_messer des Meridians r, = (1-+f’?)3/f"". Fir flache Schalen Abb 04 Dadtachvontoster Koordinntonayetces 

s ist es zweckmaBiger, den Meridian zu geben durch € = €(r) 

i (vgl. spater II. Mitteilung), wahrend die hier eingefiihrte Bezeichnung giinstig fiir rohrférmige 
' Schalen ist, fiir die also f’ beschrankt bleibt. Da 


: he 1 | € n G 
sin @ = =» cos Y = ————— a 
”)  Wdere out (see wish a0 0 
+ ist, hangen das (x, y, z)- und das (&, 7, €)-Koordinatensystem nach ie 1 
) nebenstehendem Schema der Richtungscosinus zusammen. Durch 0 Vi+f? | Vl+f2 
| Umrechnung auf die neue unabhangige Variabele ¢ statt m und , ; 
| unter Beriicksichtigung der Transformationsformeln (5,13a), er-  , | 9 7 mae 
} halten wir statt (2, 9) bis (2,11) und (2, 14) bis (2, 22) folgende Vi+f? | Vi+f? 


| Beziehungen: 
1. Gleichgewichtsbedingungen am Schalenelement (K, und K, sind die Komponenten der 
| Belastung parallel und senkrecht zur Schalenachse) 


seers Nesee | = fT? K,, 1 
dt Vie fe ( p | i Q,) fi 1 ie 1 ( ) 
— fla e+ @)] FITEP = TEP Ks, . 


— £(fM,) +1 Mo—f \TFF? 0, =0- (3) 


2. Hookesches Gesetz (ohne Querkraftverformung) 
2% Ny =D (&+ ves) » 
20 No == D (es-\- 925) 5 
) 
) 


296M, == K (4 +-0 Hs) » 
2p Mg == IS (Hp 9 X,) - 


3. Beziehungen zwischen den Verformungen und Verschiebungen' 


as Se 


1 Durch den Strich wird jetzt immer die Ableitung nach ¢ bezeichnet. 
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Cl — 4 ? (9) 
q 

4%o— F 5 (10) 

TH ere | 

A= 2S . (11) | 

fyi+f? 

4. Geometrische Nebenbedingung zwischen den Verschiebungen 
TG et ef) ds — 0 (12) | 


5. Der Zusammenhang zwischen den Verschiebungskomponenten q, parallel und q, senkrecht | 
zur Schalenachse und der Verdrehung q, einerseits und den natiirlichen Verschiebungen vin Rich- | 
tung der Meridiantangente, win Richtung der Normalen der Schalenmittelflache und der Ver- 
drehung der Meridiantangente andererseits nach (5, 13a) 


v —f'w 
= —_—_., Sra 
Nh jit f? ( ) 
fiv+w 
Se l4a 
A igpaucar (14a) | 
UE ay fs (15a) 


6. Der Zusammenhang zwischen den Schnittkraftkomponenten p, parallel und p, senkrecht 
zur Schalenachse, dem Biegemoment p; einerseits und der Querkraft Q,, der Langskraft NV, und | 
dem Biegemoment M, andererseits nach (5, 13a) und (4, 1) 


Ny —f'Q: | 
EAD g ak 13b 
Pi mf yl+f? +f? ( ) 
f' No+ Qe 
P2 i yi+f® ( ) 
P27 (15b) 
Die Hamiltonsche Funktion (5, 14) erhalt in den hier gewahlten Variabelen die Gestalt 
of / 1 {nf (Pitf'P2)? Pi+f'P, VI+f? . | 
t= (oF pdt PDAS ayy hilt yy EF a], 
2 = (P2—SP1) % 2\ (pefyitf2 Y I Df ( y?) f a (16) 
ess Ps\? Ps oy fai I | 
Fie an ee ATP al| ! 
+22 fyi+f? (Kyq+ K.q) 
und die kanonischen Gleichungen lauten 
Be ou Pore x 
Lay FE loans (o= 1,2,-3)< 
also 
Pi = —2nfyl +f? K,, (17a) 
/ ye" oe J1- fe P ———— 
Dawe) pe Bama) ee) era ta ve — 2ufli+f? Ky, (17b) 
P= yh P— (Ey) = qs — (P2—f Pi) » (17c) 
f Rie 
Posy DORI / q2 
Sie i ae Bhs | 
q Ef q f DfVi+f? r 43 » (17e) 
ene peed taal 2 


dp 
1 Damit oy Podqo —H,dqp in es podqs — H, d= tibergeht, mu® H, = H, eins wobei noch in 
H, (5, 14) i ne und f, = cos 5 zu setzen sind, 
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Die ersten drei kanonischen Gleichungen sind — entsprechend den kanonischen Gleichungen (4, 3) 
— identisch mit den Gleichgewichtsbedingungen (1) bis (3), in denen nach dem Hookeschen Ge- 
setz fiir die Ringkraft 


2% No=v Pi CSP a 1 D (1—»2) 2 18 

. fVI+F2 ( i uy) 
and fiir das Ringmoment 

2 My = v®2 + K (1—?)- Pgs 19 

ipa FV lf 1) 


gesetzt ist. Die zweite Gruppe den kanonischen Gleichungen gibt eine Verkniipfung des Hooke- 
ischen Gesetzes (4) mit (8) und (9) und der Nebenbedingung (12) an und das Hookesche Gesetz 
fiir das Biegemoment (6) mit (10) und (11). 

Die kanonischen Gleichungen (17) bilden ein lineares Differentialgleichungssystem erster 
)Ordnung von vier Gleichungen — q, ist zyklische Koordinate, und daher werden q, und p, durch 
| Quadraturen ermittelt —, fiir die Komponenten der Schnittkraft und der Verschiebung senkrecht 
zur Schalenachse, das Biegemoment um die Breitenkreistangente und die Verdrehung der 
jMeridiantangente. Sie sind die Eulerschen Gleichungen des kanonischen Variationsproblems 


lS Paarl = 
a | Dire gra == 0% 


a | 


7. Dimensionslose Gro®en. In den homogenen kanonischen Differentialgleichungen 


mF p= Daw LIF a, (1a) 
pi—» 5 p= Kv) ae oe (1b) 
g r= 5 Pe b Qs » (Le) 
i +95 a= Zl tty, (1d) 


stéren die Koeffizienten D, K und (1—»?) (vgl. Abschn. 2) den symmetrischen Aufbau. Nun 
1aBt sich eine kanonische Transformation angeben, bei der die transformierten Impulse und Ko- 
ordinaten die mit dem Symmetrisierungsfaktor 13 (1—»)/h multiplizierten p, und q, sind, 
so da das transformierte kanonische Gleichungssystem eine véllig symmetrische Gestalt hat. 
‘Mit dieser Symmetrisierung wollen wir gleichzeitig dimensionslose G1éBen einfithren. Da die 
‘dimensionslosen Differentialgleichungen wieder kanonische Gestalt haben sollen, sind die Bezugs- 
igriBen fiir die Krafte, das Moment und die Verschiebungen nicht unabhangig voneinander; die 
‘Einfiihrung dimensionsloser GréBen mu demnach durch eine kanonische Transformation ge- 
_ Die dimensionslosen Impulse und Koordinaten seien P,, P,, P; und Q,, Q,, Q;. AnStelle der 
‘Schalendicke h, des Breitenkreisradius f und der Koordinate auf der Schalenachse € werden 
die Dimensionslosen 7, 9, € wie folgt eingefithrt: 

FAG ase z__ 

wobei die Bezugsgréfen 7, und hy ein mittlerer Breitenkreisradius baw. cine mittlere Schalen- 


'dicke sein sollen. 
| Die ,,erzeugende Funktion‘‘ V, entsprechend Abschn. 5, sei 


3 1 — v2 i i 
Ve oo c Pi Q+- we Ca P3 0, . (3) 


Die Transformationsgleichung (5, 4a) lautet fiir diesen Fall zweckmafig umgeschrieben 


( 3 P.dQ, — Hal) + sora ( 


wobei der Faktor 1/4aEh2r, auch den zweiten Summanden dimensionslos macht. Dann 
ergibt sich aus den allgemeinen Gleichungen einer kanonischen Transformation, entsprechend 


3 
4 4p,— Hat) = aV, 
Bf 


o= 


4a Eh? 
y3(1—»*)_ Pp 
n 4a Eh? ” 
p= v3 —»?) Ps 
” 4a Eh?r, 
und 
Oe ee 
42 Eh?r, do aa ODPo i 
also 


Th eee eee, 
q = 3( i y”) To Q;: 


3 1—»?) 
(Cars Se n 7 Q >» 


Beare) 
1h AT 3° 


Die dimensionslose Hamiltonsche Funktion ergibt sich entsprechend aus 


H H. av 
Ei ee eee 
4a Eh? ac 
mit 
aV 3(1— 2 
Lo Hees (7) P2 Q2 + Ps Qs) » 
indem man in (5, 16) und in (7) die p,, q, durch die P, und Q, nach (4) und (5) ersetzt: 
H s | 3 ae 2/1 


Hepa oa oe (P,Q, + P3 Qs) — 


1—v An P, F 2 ae 
5 3 2 ef pd 


— 770,47 BO oyi fe? (x04 Ko): 


Der Faktor A= a ist sehr klein gegen eins. Die kanonischen Gleichungen lauten 
0 


Vo (ae wel: 
Pe 1 /3- —y2) ao leee 2K 


2 E}? Le 


pluie e’\ p> _ V30—) Vl+e?, | » 1 ¥3(1—) 0 .-—— 
Pit (% —72) P= | 10 Qo 10 P, oS Ew Tt 0? K, 


, : : 4 /1—r? Or 
P; (” oe oO, (P, ETA 


i 
ire 5 . 
v= 5 A—— Q, 7 Of 
/ n’ Qo’ 1 /1— no’ P, 
Go = 78 | yaa ee ee eee Pe 


' ‘ 0’ —v?) /1+ 9’ 
es eae ie aE 
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Dies sei die endgiiltige Gestalt der Schalengleichungen, deren Integration in der II. Mitteilung 
besprochen werden soll. 

Riickblickend kénnen wir feststellen, da uns die Variationsrechnung, insbesondere der 
Ubergang zum kanonischen Variationsproblem, zwangslaufig auf die bei unserem Schalen- 
problem zu bevorzugenden Variabelen gefiihrt hat, die unmittelbar anschauliche Bedeutung 
haben. Durch Verwertung eines bekannten Integrals gelang es, das Integrationsproblem von 
‘der sechsten auf die vierte Ordnung zu erniedrigen, wodurch nun auch geklart ist, weshalb 
Reifner und Meifner ihre beiden simultanen Differentialgleichungen zweiter Ordnung aufstellen 
-konnten. Die kanonische Gestalt des Differentialgleichungssystems ist, wie wir im folgenden 
zeigen werden, fiir die numerische Integration vorteilhaft. Vor allen Dingen lat sich sehr ein- 
fach durch Variation der Konstanten eine Partikularlésung ermitteln, so da® wir nicht mehr auf 
die Kenntnis des zugehérigen Membranspannungszustands angewiesen sind, der ja keineswegs 
bei allen Schalen vorhanden ist. 


(Eingegangen am 11]. August 1950.) 


Anschrift des Verfassers: Dr.-Ing. H. Miinz, Technische Hochschule Braunschweig, 
Institut fiir Technische Mechanik. 
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Neuartige Behandlung der Poissonschen und der inhomogenen 
Bipotentialgleichung bei rechteckigen Bereichen mit Anwendung auf 
Probleme der Torsion und der Piattenbiegung. 


Von Th. Schade. 


1. Einleitung. Zur Behandlung potentialtheoretischer Randwertprobleme geht man be-4 
kannterweise so vor, da®B man eine solche Koordinatentransformation der kartesischen Koordi-+ 
naten ansetzt, da eine der Variablen auf den Rand konstant wird. Gelingt nach Transfor--+) 
mation der Laplaceschen Differentialgleichung A1® — 0 die Trennung der Variablen, so versucht |} 
man durch Reihenansatz der so gewonnenen Potentialfunktionen die Randbedingungen zuij 
befriedigen. Analoges gilt fiir die Bipotentialgleichung J A® = 0. Hat man es mit der inhomo- |}, 
genen Differentialgleichung 

A®= f(x, y) 
zu tun, so itberlagert man eine Partikularlésung der inhomogenen Differentialgleichung mit) 
Lésungen der homogenen Differentialgleichung so, daB wiederum den Randbedingungen geniigt | 
wird. | 

Diese beiden Bedingungen, Erfiillung der Differentialgleichung und Erfiillung der Rand- | 
bedingungen lassen sich in vielen Fallen nicht exakt durchfiihren, so da man gezwungen ist, 
irgendwelche Naherungen zuzulassen (wie z. B. das Beispiel der vollstandig eingespannten recht- 
eckigen Platte mittels eines Ansatzes von Fourier-Doppelreihen zeigt). 

Es soll nun hier versucht werden, durch eine Koordinatentransformation in schiefwinkligen | 
Koordinaten die Differentialgleichung und die Randbedingungen ganz oder teilweise exakt zu 
erfiillen und dafiir Unstetigkeiten héherer Ableitungen der Lésungen, die physikalisch fiir das | 
Problem nur eine untergeordnete Rolle spielen, an gewissen Stellen in Kauf zu nehmen. 

Diese Methode wird vornehmlich an Plattenproblemen durchgefihrt, wobei zunachst bekannte 
Aufgabenstellungen untersucht werden, um einen Vergleich der Ergebnisse zu erméglichen. 

Dazu sei noch folgendes bemerkt: Die bisher ttblichen Methoden der Fourier-Doppelreihen | 
oder der anderen numerischen Verfahren auf der Grundlage der Differenzenrechnungen haben | 
gerade beziiglich der Berechnung der physikalisch interessierenden GréBen gewisse Nachteile. 
Bei den Reihendarstellungen sind es Konvergenzschwierigkeiten, und bei den numerischen Ver- 
fahren sind es Fehler in der numerischen Differentiation, abgesehen von den meist recht kompli- 
zierten Berechnungsverfahren. Es soll hier dagegen versucht werden, geschlossene analytische | 
Ausdriicke zu finden, die das physikalische Problem zumindestens fiir die Praxis genau genug 
darstellen und nicht die oben erwahnten Nachteile in sich tragen. 

In den zur Diskussion stehenden Fallen soll es sich um eine rechteckige Berandung handeln, 
wobei speziell die quadratische Berandung untersucht werden soll. 

Um zu einer entsprechenden Parameterdarstellung zu kommen, kann man sich der Ko- | 
ordinatentransformation bedienen, die R. Grammel! in einer Arbeit iiber ,,Erweiterte Kreisfunk- | 
tionen® angibt. 

Ausgehend von der algebraischen Gleichung 

ay (n ganzzahlig) 
definiert R. Grammel die erweiterten Kreisfunktionen durch 
== Cnii(n) os || 
y= Sin(n)v | 
(Abb. 1), wofiir in der erwahnten Arbeit Reihenentwicklungen angegeben worden sind. Gleichung 
(1, 1) stellt fiir gerades n eine geschlossene Kurve dar, die fiir n—> oo in ein Quadrat ausartet; 
die Periode 2 7, liegt fiir n—> oo in den Grenzen 2% <2, <8. 


Ks liegt nun nahe, diese Funktionen zu benutzen, um eine Parameterdarstellung der kar- 
tesischen Koordinaten in der eingangs geforderten Form aufzuschreiben, was schon durch die 


1 R. Grammel, Ing.-Arch, 16 (1948), 5. 188 sowie 18 (1950), S. 250. 


(1, 1) 
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einfache schiefwinklige Koordinatentransformation 
x= opj(n)v, | 
y=oe Gin(n)v | (1, 2) 
erreicht wird (Abb. 2). Mit n—> co erhalt man fiir o = konst. eine Sch 


Quadrate und fiir v = konst. die Geradenschar durch den Nullpunkt. 
Fiir eine rechteckige Berandung kénnte man analog die Parameterdarstellung 


ar ineinanderliegender 


x—= aon] (n)v, 
y= bo Sin (n)v 
benutzen, wobei dann, wenn man die geometrischen Bezeichnungen der Ellipse iibertragend an- 
wendet, a und 6 die Halbachsen der geschlossenen Kurve 
@ = 1 darstellen. Der Einfachheit halber soll zunachst das 


~Quadrat untersucht werden. 


g 
Abb. 1. Graphische D <ti i * Abb. 2: 
raphische Darstellung der Funktionen Cos (n)v, Sin (n)v. Das schiefwinklige 0, v-Koordinatensystem. 


2. Transformation der Potentialgleichung und ihre Lésung. Geht man den in der Einleitung 
erwahnten iiblichen Weg, so hat man (1,2) auf die Potentialgleichung 


ed PP 


AP=Fa ae: (2, 1) 
anzuwenden. 
a, O@ 6 p= Ov 1 
Mit 2—-Goj(n)""v, j=—— Sin(n)», | 
a) O : n—] Ov 1 . (2, 2) 
fae Sin (n) ~ v, ree | Co} (n) v, 
erhalt man fiir die Differentialgleichung (2,1) 
FE [Coj fn)?" ? w+ Sin (my? 0] FG oy [Sin (mp4 Cof (m)? 0] + 
o®@ 1 ’ Tl ~ al Sli 
= Aa [Sin (n)"v Co} (n) v — Coj (rn) v Gin (n) vo] + | ae 
+5 = Sin (n)"—* v Coj (n)"~* v [Cof (n)? v + Gin (n)? v] + 
“Loe x [Gof (n)"—2v Sin (n) v — Gin (n)”~ 10 Cof (n) oJ . 
Mit dem Produktenansatz 
@, = o* F, (v) 
zur Trennung der Variablen und der Substitution ade 
U;, (v) =F, (v) (Gin (n)? v + Co (n)2») 7 
ergibt sich die Differentialgleichung fiir U, (v) 
PUL) 1 yy, (v) [(n—1) Gin (n)"~? v Cof (n)"~* 0+ | re 


+ (42 — 1) (Gin (n)? v + Coj (nm)? oy] = 0, 
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mit der allgemeinen Lésung 


U, (v) = Gin (n)? v + Cof (n)? v ay Aare tg (Tg (n) v). 

Fiihrt man in (2,4) die Substitution Tg (n) v = tg u ein, so erhalt man fiir ganzzahliges 1. 

und n—> co unter Einfithrung der komplexen GréBe z= x +-iy, baw. == x — iy die bekann- 
ten Potentialfunktionen in kartesischen Koordinaten 


A _ Re| 1 
a) SROs coi ee 


b) A= 0:@) = arc tg (Zg (n) v) = arc tg ae 


ce) A> 0:0) = Fe | aA 


Fiir die spateren Beispielrechnungen ist es sehr viel praktischer, die Darstellung der Funk- | 
tionen im (0, v)-Koordinatensystem | 


GA = beizubehalten. 

Aus (2,5) erkennt man, daB die | 
unter a) gegebenen Potentialfunk- | 
tionen einen Pol bei 9 = 0 besitzen, | 

+} + wahrend die unter c) angegebenen | 
AL einen Pol bei 90 =oo besitzen; diese | 
TD | Tatsache bestimmt ihre weitere Ver- | 
4 3 vp wendung bei vorgegebenen Rand- | 
0 7 y wertproblemen der Potentialtheorie. 
Aaé 4 Unter b) ist die von 9 unabhangige 
AB Potentialfunktion angegeben. 
A=6| | | 
el al =| 


‘Diagonal- 
Nien 


Abb. 4. Einteilung der Bereiche fiir das 
Randwertproblem, 


Man erhalt allgemein fiir alle in 
Polarkoordinaten r, qm bekannten 
Potentialfunktionen durch folgende 
Zuordnung die entsprechenden im 
(0, v)-Koordinatensystem: 


r= 9 | Sinn) v Cop %(n)v, | | 
v 


lye y— arc tg Tq (n)v. | 
Le 
Abb. 3. Graphische Darstellung der Funktionen F,. a; Fy Me Die Bipotentialfunktionen ue erhalt | 


i man durch Multiplikation der Po- 
tentialfunktionen mit 2? + y?= zz, womit sich analog zu (2,5) ergibt 


a) Neat res he) zB | 


Jm { Ac1” 
b)EAa= 0; = 22 are tg = = (x2 + 9") are tg, 


(2.7) 
c) A= 05 Pam faa | 
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Fir 2, = 2 erhalt man die von 9 unabhangige Bipotentialfunktion, fiir 1 ~ 2 diejenigen mit 
einem Pol bei 9 = oo und fiir A, < 2 solche mit einem Pol im Nullpunkt. Die unter (2,6) ge- 
gebene Zuordnung behalt natiirlich auch hier ihre volle Giiltigkeit. 


Mit ®, (0,0) = o F, &) = o" (Fa +i Fy) (2,8) 
erhalt man die in Abb.3 dargestellten Funktionen fiir n—> oo. Teilt man weiterhin, wie in Abb. 4 
angegeben, das Quadrat in 4 Bereiche I bis IV, so interessieren fiir doppelsymmetrische Pro- 
bleme nur Funktionen, die zu den Mittellinien v — 0; 4; 2; 6 symmetrisch sind; analog erhalt 
man fiir rein unsymmetrische Probleme die hierzu unsymmetrischen Funktionen. 

Diese kann man sich aber leicht aus den in Abb. 3 gegebenen Funktionen konstruieren, 
wenn man sie stiickweise bei v = 1; 3;5; 7 zusammensetzt. Der absolute Betrag der Ordinaten 
an der Flickstelle ist fiir beide Kurventeile mit dem gleichen Parameter A von vornherein gleich. 

Abb. 5a und 5b zeigen die somit geschaffenen Funktionen fiir eine Viertelsperiode. 

Diese Zusammensetzung bedingt, daB man bei dem entsprechenden Randwertproblem sich 
nur um die Bereiche I und II zu kiimmern braucht, wobei man zusitzlich gewisse Stetigkeits- 
bedingungen der Ableitungen an den Diagonallinien zu erfiillen hat. Es wird sich zunachst 
um doppelsymmetrische Probleme handeln. 

Fiir die Potential- und Bipotentialfunktionen im Bereich I ergeben sich folgende Polynom- 
ausdriicke im (0, v)-Koordinatensystem: 


ee all e2 Ge saga ee 4 eh ar Ay (= 10% ' 


Oe ie 0 ee Pry ; v2 eed tae -+(— 1)? (2n+-1) a 5 


(2n ile ar 4 oo 1h one ? 
V, = 02 (1 +0) O. 
Die Funktionstafeln! kénnen rekursiv errechnet werden. Wie man leicht zeigen kann, gilt 


ees as I) UP emcee rate re | 


YS tana iy sav oe a || 


Besonders wertvoll ist die Tatsache, daB sich die Ableitungen der Lésungsfunktionen, die 
zur zahlenmaBigen Berechnung der physikalischen GréBen notwendig sind, wieder als Linear- 
kombinationen der Potential- und Bipotentialfunktion angeben lassen, so daf man sich der 
gleichen Tabelle bedienen kann. 

Es ergeben sich im einzelnen folgende Ausdriicke, wobei es sich auch in Zukunft von selbst 
verstehen soll, daB man Real- und Imaginarteil fiir sich zu nehmen hat: 


Ib 


(2,10) 


oot = 10, =i, a= da 1) Q—2) Oa =i SF 

PAAR 1) Os =—— FH, geg = ARNO 2) 94s = iggy 
$e th a= Oi 5 

ae = A+) Pas Ort. ao ot Met mee oS aly 

=i [A +1) Pir — Pin), Fe AED AY) Fs +r] 

aa (+1) DY +20). a= iA(A +1) [(A—1) Yas +O], 

= A+Y[-AM 2 420), gga = 40 +D[-O-) Ma +O], 
Japs KIA FY Ys, 


A U4 (A +1), 


1 Diese Funktionstafeln sind im Originalexemplar der Habilitationsschrift angegeben. 
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Abb.5b. Symmetrische Potentialfunktionen. 
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Fiir die nun folgenden Beispielrechnungen soll allgemein festgesetzt werden, daB Potential* 
funktionen mit dem Koeffizienten a, und die Bipotentialfunktion mit den Koeffizienten b, 
versehen werden. 


3. Beispielrechnungen. a) Torsion eines prismatischen Stabes mit 
muadratischem Querschnitt (Abb. 6). Fir das Problem des tordierten prisma- 
ischen Stabes lautet die Differentialgleichung 

AD=—— K=— 26 (3,1) 


mit G als Gleitmodul und w als Verdrehwinkel, woraus sich 
lie Schubspannungen ergeben zu 
__ a@ aD 


ee — Be Ea) =| eee 0 (3,2) 


Die Randbedingung lautet 


P= fiir O= 0- (3,3) Abb. 6. Querschnitt des prismatischen Stabes. 


Dazu kommt die zusatzliche Randbedingung aus der Schubspannungsfreiheit der Mantelflache 


0® - 
ee eae 980 iby Weil 
und die Stetigkeitsbedingung der ersten Ableitung an den Diagonalen 
ob Aa®@ 
Cee oP a ae 3,4 
ar ie -Qe= ib. (3,4) 


Man bedient sich bei der numerischen Rechnung der im vorigen Abschnitt definierten symmetri- 
schen Potentialfunktionen. Da man die Stetigkeitsbedingung (3,4) nicht exakt erfiillen kann, 
oll naherungsweise die Bedingung 


od . o® . 
eS fir 9 =a/2, v= 1 
sestellt werden. 
Mit dem Ansatz 
K a; a ay 
P= F-Pt ag OF + Des +25 Pea +5t Pas +5i Poo — 795 Per] 

rhalt man zur Erfillung der Bedingungen (3,2), (3,3) und (3,4) durch Koeffizientenvergleich 
n den v-Potenzen 5 Gleichungen mit 5 Unbekannten, woraus sich fiir die Koeffizienten ergeben: 


208 3 1 ie ae ee: 
COTE TLIO ST a a reat me ana SSM IS 
Jieraus erhalt man fiir die maximale Schubspannung und das Torsionsmoment 
Malt 
Tmax 40 Koy ~) 
s 588 ; 
M, = Er Ko; 
M, 
der Tapp = GY Fs 
ait F als Querschnittsflache und b als Seitenlange des Quadrats . 
egeniiber einem Faktor von 4,81, der sich aus der Fourter- Abb. 7. Definition der Spannungen am 


Reihenmethode ergibt. 
b) Die quadratische diinne Platte mit konstanter Belastung 
Abb. 7). Die Differentialgleichung der diinnen Platte mit konstanter Belastung p, lautet 


AAw="2—K, (3,5) 
venn w die Durchbiegungen bedeutet, N die Plattenkonstante 
cet EO! 
= 12 (1 —) 


ait h als Plattendicke, E Elastizitaétsmodul und y als Querkontraktionszahl. 
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Fiir die Berechnung der Spannungen und Querkrafte gelten die Beziehungen 


Belity: io ew. aw ie a ) 

6, = My. Mm, = Nie | a) Pr=— Nw Aw, | 
122 0? w 0? w (a) 
Cy= My, My = N53 Ly aa) | cpp aN a Aw, 3.0) 
12 02 w 
Sy peso Gh pee | 
Txy = Mxy Fa» Mey (1 vy) N PCE: } | 


a) Fir die an allen Randern frei aufliegende diinne Platte konstanter Belastung gelten die 
Randbedingungen (Bereich I) 


0 fiir O=o- (3,7)) 


0? w 

Hj — Ny eas 

Als Zusatzbedingung soll noch die Stetigkeit der Spannungen an den Diagonalen, was mit der) 

Forderung (6,),—1 = (0y),—1 gleichbedeutend ist, erzielt werden, oder auch | 
ew dw 
Ax2 Oy? 


fir aoe. (3,8)) 


Weiterhin sei zusatzlich das Verschwinden der Querkrafte im Nullpunkt gefordert, d. h.| 


d= 9, == Otro — 0. Mit, deme Anzarz 


K 4 b. as 
ah of (1 (Sia | ay 04 | bo 05 Peo t5 Pest Pea to Pxs 


erhalt man zur Erfiillung von (3,7) und (3,8) durch Koeffizientenvergleich in v 6 Gleichungen 
mit 5 Unbekannten. Die iiberzahlige Gleichung wird dabei von selbst erfiillt. Das Ergebnis lautet 


24 16 1 1 
Oo = exo i al Done a, =a b= => 
und man erhalt fiir die maximale Durchbiegung 


a Po 5* 
Wmax — 9,00469 N 


(statt des Faktors 0,00406 nach der Rechnung mit Fourier-Reihen). 
B) Analog erhalt man fiir die an allen Randern eingespannte Platte mit den Randbedingungen | 
Ow 


w== 0! cpates Oe tire =) — 0, (3,9) 
und den Zusatzbedingungen 
e GE % 
en a A fir v=1, q=q=0 fir p=0 (3,10) | 


und dem Ansatz 


K f p 
w= o* (1 = v*)2 = a, 05 4 by 05 Moot ot Past ts Pes +7 Pes 
0 


die Koeffizienten 


8 8 2 
Oy = Sar iy arses bee a——l1, a 
und 


— Po bt 
Ue = 000156 “WO 


(statt des Faktors 0,00128 nach der Naherungsrechnung mit Fourier-Reihen). 


Die hier auftretenden Differenzen sind auf die Unstetigkeit der Querkrafte q,, qy an den | 
Diagonalen zuriickzufiihren, womit auch das Gleichgewicht gestért wird. Hieriber an aim| 
SchluB der Arbeit noch einiges gesagt. | 

Ks sei hier besonders erwahnt, daB die angegebenen Lésungen eindeutig sind, da die Auflésung 
des Koeffizientenschemas fiir eine Erweiterung des Ansatzes mit Potential- und Bipotential- | 
funktionen die Koeffizienten dieser Funktionen exakt zu Null macht. Weiterhin sei bemerkt, | 
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das, wie man leicht allgemein zeigen kann, die Stetigkeitsforderung der zweiten Ableitungen 
die Stetigkeit der ersten Ableitungen nach sich zieht. Mit Hilfe der Differentiationsformeln 
(2,11) lassen sich leicht die Spannungen und nach Mohr die Hauptspannungen nach GréBe und 
Richtung bestimmen. Es ergeben sich die Trajektorienbilder Abb. 8 und 9. 


Lugspannungen 
——— —Drickspannungen 


— Zugspannungen 
— — — — Druck Sspanningen 


Y, © 
Ba ip 
Lf) LIE 
LARIT AES 


Abb. 8. Spannungstrajektorien der an allen Rindern frei Abb. 9. Spannungstrajektorien der an allen Réandern 
aufliegenden Platte fiir z > 0. eingespannten Platte fiir z> 0. 


4. Die allseitig eingespannte diinne Platte mit einer Lastfunktion als Potenzreihe achten Gra- 
des in 9. Die Lastfunktion sei gegeben in der Form 


8 o\n 
P= To = ea() . (4,1) 
Dann lautet der Ansatz fiir die Lisung der Differentialgleichung 
K 0? 0° 6 Q 3 0 5 2 
== Se! RO abi Se Does! = c ee ee 
b! 0 Q4 we 252° 7 “8 op TE A etsy 2 peel 
eee ey! 2 go | Ki Leer 
7 “6 oF 11 Cr 01 2) “808 61)" 04 Sets Wo 
an 8 
11 
bn | K | 
+ by 02 Wg ee ee Pg calle oie ri Sra ean 
i 0 


Das sich nach Einsetzen der Stetigkeitsbedingung (3,10) ergebende Gleichungssystem ist in 
Tabelle la und 1b angegeben. 


Tabelle la. Koeffizientenschema der linken Seite des Gleichungssystems. 


a by a, | b; | by b; os b, bs b, A2 3 by 
a 1 1 iL 8/5 il 2/7 1 16/9 1 2 1 1 Il 
II %, 4 8 6 2 8 16 10 2 12 13 13 
Til 1 6 | — 8 5 6 | —28 —48 —21 10| —66 —78 —p4, 
IV 119? 24 20 30 168 336 216 90| —660 —858 | —594 
Vv I 5 30 70 112 42 | —180 495 715 275 
VI —l10 | —90 280 560 252, |\—1260 3960! 6435 2475 
VII 1 10 | —28 |—112/3 4.2 420 924) 1716 132 
Vill 10 | —56 —112 168 2100|—5544| —12012 | —924 
Ix 1 —27 | —210 495 1287 | —297 
xe —54 | —630| 1980 6435 |—1485 
XI 1 18} —66 —858 462 
XII : 18} —132 —858 462 
XIII 1 is —ll1 
XIV i —I] 


Wie man erkennt, ist die Auflésung des Gleichungssystems sukzessive méglich. Es kann sich 
im Verlauf der Rechnung héchstens um die Auflésung von zwei Gleichungen mit zwei Un- 


g* 


— 
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Tabelle 1b. Koeffizientenschema der rechten Seite des Gleichungssystems. 


Co Cy cs. Ge Op ce é, Oz (Os | 
I 15 15/2 9/4 27/28 13/28 | 25/112 | 31/224 31/352 | 157/2640 | 
I 60 75/2 27/2 27/4 26/7 | 225/112 | 155/112 31/32 | 157/220 
III —30/2 | 15/4 —9/4 Siy7 15/28 45/224 5/32 | 1/20 
IV 245 is —45/4 es 15/4 | —45/28 | 45/32 | =e 
Vv 15 15/2 15/4 15/4 —15/8. | —45/16 } 15/16 |—=3/16 @y 
VI 15/2 15/2 45/4 15/8. | —A5/8mi= 105/16 P= =3/2 
VII = 14 —3/4 1/2 5/4 15/16 21/16 
Vill —3/4 1 15/4 15/4 105/16 
1s —1/28 | —15/112 | —45/224| —15/32 3/16 
x —15/112 | —45/112| —45/32 3/4 
XI 1/224 1/32 | 1/20 
Us 1/32 | —1/10 
XII 1/880 | 
XIV | 
bekannten handeln. Es sind 14 Gleichungen mit 13 Unbekannten, so daB die iiberzahlige Glei- | 
chung als Rechenkontrolle dienen kann. Das Ergebnis ist in Tabelle 2 angegeben. | 
Tabelle 2. | 
Co ay Ce C3 C4 es Ce c; es 1 
a 18 7,142 857 2,761 904 1,2 0,485 714 0.164072 | —0,117396 | —0,364635| —0,658407 | 
by —30  |—10,142860 | —3,428571 | —1,285714 | —0,380952|} —0,012059 0,328 386 0,638 528 1,028 571 
a, Bais = 58 0,75 1 0,982143| 1,178571 1,458 333 1,933 333 
by a5 —1,5625 1,875 —1,6875 — 1,583 333 | —1,346726 | —1,508928 | —1,822916| —2,416667 
by 1,875 1,5 1,125 1 0,758 929 0,910 714 1,1875 1,70 
bs —0,281250 | —0,1875 —0,093 750 | —0,125 —0,0379 464 | —0,120536 | —0,234375| —0,425 
ag 0,0357143| 0,0558036 |—0,0133929 | —0,078125 | —0,175 
b, — 0,0251116 |—0,0167411 |—0,0150670 |—0,0200893 |—0,3597895 | 0,0193718 |—0,0376674|  0,0683036 
bs — 0.2558 035 |—0,0111607 |—0,0130 208 | —0,0208 333 | 
by 02465030 | 0,3558036 |—0,2651042 0,2104 167 
Aye — 0,2113 636 
oe 0,4268286| 0,4321944| 00,4626 002 0.2220 352 
Bea 0,4317066 | 0,4380479 | 00,8739 820 052260417 
3 3 5 lh 
a,= 7 4 + — bs, My = 596% to br, b,=— 15 ¢,. 
1 1 weet 
Cet oo 10 = 994 6? bs = 59 fa? 
3 5 1 9 i 
O; = x53 — 7 Os 4 = 355% — qq Oo» bio = — g50 C8: 


Aus der Anzahl] der méglichen Kombinationen sei hier nur 
noch ein instruktives Beispiel herausgegriffen. Die Belastung 
p sei in der Form gegeben (Abb. 10) 


p=m(l— 2). 


20 


Die Konstanten lassen sich aus Tabelle 2 berechnen, worauf hier verzichtet werden soll. Es _ 
soll hier vielmehr die Gleichgewichtsbedingung bzw. die Querkraftstetigkeit an den Diagonalen | 
gepriift werden. Fiir die Auflast p erhalt man allgemein das Integral 


+Qo 
Pi=2 | ople)de 


—@o 


und fiir das Integral der Querkrafte am Rand 0 = a, 


+@o 
g b(A +1 
Po= | le A «) dye Dialed pay on 
@ 80 js A v=l 
—eo 


Die Gleichgewichtsforderung lautet nun 
P4= Po. 
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Fir den oben angefithrten Fall erhalt man 
FO) BP, ae Pa = 0,105 py @2 . 


Die unwesentliche Differenz zwischen P, und Po laBt erwarten, daB auch die Querkrafte Ix 
und q, eine gute Ubereinstimmung aufweisen werden. In Abb. 11 sind beide Kurven gezeichnet 
und bestatigen die Vermutung. 


5. Schlu8betrachtung. Die zunichst durchgefithrten recht cinfachen Rechnungen an den 
schon bekannten Plattenproblemen haben gezeigt, da& man hierbei qualitativ gute Ergebnisse 
-erzielen kann. P,P, 

Das Problem besteht weiterhin darin, auch der # 


ee tzten physikalischen Bedingung der Stetigkeit der 


~Querkrafte an den Diagonalen Rechnung zu tragen. 
Im Abschn. 4 wurde gezeigt, daB es auch gewisse Last- 
falle gibt, die dieser Bedingung in sehr guter An- G6 


naherung geniigen. 06 


; Oo (A)v=1 “FC @ 3 
CC 


Man kann nun auch in allen anderen Fallen dieser gt 
Bedingung zumindestens punktweise oder im Mittel ga g 
-geniigen, wenn man dafiir die Lastfunktion durch ein 0 wld 


Polynom annahert. Alle anderen Bedingungen werden — -gz 
‘dabei beibehalten. Die bisher Erfolg versprechenden — -g¥ 
Ergebnisse dieser Rechnung sollen einer folgenden 
Arbeit vorbehalten bleiben. 


Es sei hier noch bemerkt, daB es auch miglich ist, Lastverteilungen der Form 


— al In 2 (n > 0) 


tv) 


Abb 11. Darstellung der Querkrifte q,., dy firv = 1 


zu rechnen, so daf man eine ganze Reihe von Uberlagerungsméglichkeiten zur Erfiillung der 
Stetigkeit der Querkrafte zur Verfiigung hatte. 

Was neben der leichten Auflésung der Gleichungssysteme noch besonders zu erwahnen ware, 
ist die Tatsache, da® eine Anderung der Belastungsfunktion lediglich die rechten Seiten des 
Gleichungssystems (Tab. Ib) andert. Wenn man also einmal den Rechengang zur Auflésung 
des Systems (Tab. la) schematisch festgelegt hat, liegt es auch fiir die rechten Seiten fest und 
kann unter Zuhilfenahme einer Rechenmaschine schnell erledigt werden. 

Diese Methode ist nun weiterhin auf eine unendlich ausgedehnte Scheibe mit quadratischem 
Loch und symmetrischer Belastung mit Erfolg angewendet worden, wofiir die Potential- und 
Bipotentialfunktionen verwandt wurden, die ihren Pol im Nullpunkt haben. Dieses Problem 
wird eine spatere Arbeit von H. Montag behandeln. 

Die vorliegende Arbeit ist 1951 am Institut fiir Mechanik der Techn. Hochschule Aachen 
als Habilitationsschrift entstanden. 


(Eingegangen am 26. September 1950.) 


Anschrift des Verfassers: Dozent Dr. Th. Schade, Wiirselen bei Aachen, Kaiserstr. 69. 
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Die achsensymmetrisch belastete dicke Kreisplatte 
auf elastischer Unterlage. 


Von I. Szabo. 


1, Einleitung, Uber das Problem der Bestimmung der Druckverteilung in und unter elasti- 
schen Tragern, insbesondere Balken und Platten, auf elastischer Unterlage ist seit Zimmermann a} 
eine groBe Anzahl von Publikationen? erschienen. Jeder dieser Arbeiten (soweit dem Verfasser| 
bekannt), enthalt mindestens eine der folgenden drei Kinschrankungen: 1. Der Trager ist ,,un-| 
endlich diinn‘‘, so daB seine Deformation beim Balken durch die Verbiegung der Achse, bei der: 
Platte durch die der Mittelebene beschrieben werden kann; 2. der Trager ist ,,unendlich lang‘¢, 
so daB man sich um die Erfiillung der Spannungs- und Deformationsbedingungen an den Randern 
nicht zu kiimmern braucht; 3. die sog. Zimmermannsche, in Wirklichkeit auf Winkler? zuriick-, 
gehende Annahme, die besagt, dafs die Hinsenkung der Unterlage und somit, bei Voraussetzung. 
stetiger Berithrung, auch die Durchbiegung des Tragers an der Berithrungsstelle dem dort herr-| 
schenden Druck proportional ist. Da® uns diese Annahmen mehr oder weniger von der Wirk- 
lichkeit entfernen, ist offensichtlich: Ein Schornsteinfundament von 2 m Dicke und 12 m Durch- 
messer ist weder ,,unendlich diinn‘‘ noch unendlich lang, und auch gegen die unter 3. ange-| 
fihrte Hypothese 1a8t sich der Einwand anfiihren, das die Einsenkung an einer Stelle nicht. 
nur von dem 6rtlichen Druck, sondern von allen von der Unterlage aufzunehmenden Lasten 
abhangen muf. 


2. Problemstellung, Eine Kreisplatte vom Halbmesser a, der Hohe h, Schubmodul G und der 
Poissonschen Zahl y ist in der Ebene z= 0 durch die achsensymmetrische Spannung o = f(r) 
belastet und liegt in der Ebene z= h auf dem elastischen Halbraum vom Schubmodul G, und 


Abb. 1. Rotationssymmetrisch belastete Kreisplatte auf Abb. 2. Volumenelement der Platte bzw. Unterlage. 
elastischem Halbraum, 


der Poissonschen Zahl y, (Abb. 1). Gesucht sind unter Zugrundelegung der Hookeschen Gesetze | 
die zur Charakterisierung des Formanderungs- und Belastungszustandes notwendigen Verschie-| 
bungen 0 = e(2,r), = (z, r) fiir die Platte bzw. 0, = 0,(z, r), ¢,—¢,(%, r) fiir die Unterlage 
(Abb. 1) und die Spannungen o,, 0,, 0, bzw. 612, 011, 01, (Abb. 2). Die nunmehr zu formulierenden 
Randbedingungen werden in der linearisierten Elastizitatslehre iiblicherweise fiir den unver- 
formten Zustand unter den Voraussetzungen aufgestellt, daB 1. die Platte in jedem Punkte 
ihrer Grundflache (z= h) mit der Unterlage in Kontakt bleibt, wodurch freilich der Lésung 
nur solange eine physikalische Realitét zukommt, als zwischen Platte und Unterlage Druck- 
krafte auftreten (was bei den iiblichen Belastungsarten meist der Fall sein wird), und 2. daB die | 
Oberflache des Halbraumes vollkommen glatt ist, so daB® sie dort keine Schubspannung auf- 
zunehmen vermag. (Die Beriicksichtigung derselben erfolgt in einer dieser bald folgenden 


1 H, Zimmermann,.Die Berechnung des Eisenbahnoberbaues. Berlin 1888. 
5S. z. B. E. Pflanz, Ing.-Arch. 12 (1941) S.201 und die dort angefiihrte Literatur. 


° E. Winkler, Die Lehre von Elastizitaét und Festigkeit. Prag 1867, 
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| Mitteilung.) Unter Zugrundelegung der in Abb. 1 gezeichneten Koordinatensysteme sind dann 
folgende Randbedingungen zu erfiillen: 


I. Fiir z=0 und alle r<a ist OS =f {r) aes 0. (0, r) > (1) 

| t= 0 = 7(0;7r) (2) 

t(h,r) = 0 fir’ sa, (3) 

. 7,(0,r) = 0 » aller, (4) 

Rebun’ <— kh. baw. == 0 ist Cr One =ve( is.) ape tte (5) 
Dy COIS 5) aes 2.0 

a0.) = 1 0 esi (6) 

C= 0 1 ZB, > a 

Greece anal ailens 0 leh tae eo) (7) 

| a) 6, (2, a) . (8) 


DaB die Randbedingung (5) im Sinne der linearisierten Elastizitatslehre zulassig ist, zeigt man 
-folgendermaBen: Streng genommen ist ¢ (h, r) = €,(0, r,) und r+ 9 = r,+ ¢,, woraus nach dem 
Taylorschen Satz 

OR) C00 — 0) = (0.7) 5 (0 0) 4 
folgt, und da 9 — 9, und d{,/dr kleine Gréfen erster Ordnung sind, ist die Randbedingung (5) 
bis auf kleine GréSen héherer Ordnung richtig. 
Mit der Abkiirzung fiir die Gesamtdehnung 


Q IGF) 06 
Jo eee (9) 


sind die Spannungen bekanntlich! durch 


_ OC ts VE Pe 00 oc 
la Cre area E era | 


do VE 0 VE 
== eS = | 
Cre G(s: en : Ce ie = 


(10) 


gegeben, wahrend die Gleichgewichtsbedingungen ohne die Beriicksichtigung der Massenkrafte 

(wortiber bei einer anderen Gelegenheit berichtet werden soll) das Bestehen der Differential- 
gleichungen 

00r Or —G Ot 00: Ot @ 

or ! r ‘ Pace Oz Ae arose qt) 

fordern. Durch Eliminationsprozesse? bekommt man aus (9) bis (11) die beiden Differential- 

gleichungen 


Ty (12) 
| und 
| (L—2n)(4o— 4) + = 9, (13) 


2 2 
| wobei A= = : - ; ee den Laplaceschen Operator bedeutet. 
| 


3. Bestimmung der Verschiebungen fiir die Platte, Aus (9) gewinnt man durch Differentiation 


und somit 


) ede = ide 06 0 ne 14. 
| AC=35, i (3: | Or Oz 33 or © 
Wir fihren ein 
2 
ga oe ee (15) 
Zz or 


1S. z. B. Handbuch der Physik, Bd. VI, S. 81. 
2 A, Nédai, Elastische Platten, 5. 310. 
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woraus 
Os ae "0 GG Os = a9 OG (16) 
Or oraz Or oz Oz? “ord, 
folgt. Setzt man (15) und (16) in (14) und den so erhaltenen Wert fiir 4¢ in (12) ein, so ergibt sich 
d« s Os ) 
20— N= A—29(F +5). (17) | 


Ebenfalls aus (9) flieBt die Beziehung 
OO Melee eke f Ot one e 
ar.) por “Or 72k “tords 


die mit (16) 


de 0 Os 
Ag— or a re 9 Oz 
liefert, womit aus (13) 
0 Os 
2(1L— n+ 1-25 =0 (18) 


hervorgeht. Mit Hilfe der aus (17) und (18) ermittelten d¢/dz und d¢/dr erhalten wir schlieBlich 
die bekannte Beziehung 
Ules), slace Oe 


22 ' r @r aggre Ae xl (19) | 
Errechnet man hingegen aus (17) und (18) ds/dz und @s/dr, so kommen wir zu 
Cia TREC NCEE ES Age biz») Ce sia 
an 5 er * or aera or ig © 
d.h. 
\ 
As — eral (20) 


Der erste Schritt zur Ermittlung der Verschiebungen ¢ und @ ist die Lésung der Differential- 
gleichungen (19) und (20). Der Produktansatz 


e= Ri(r) Z(z) baw. s= R(r) Z(z) 
fithrt aus (19) und (20) mit willkiirlichem 2? auf die gewéhnlichen Differentialgleichungen 
PRY (r) TR (rn) Are R or) = 0; ZG) 2Z (sz) = 0 
bzw. 
rR’ (r) + rR’ (r) + (a — 1)iR(r) = 055 22" (ze) Ze) 0; 


deren fiir r= 0 regular bleibende Lésungen J, (Ar), e=“* baw. J,(Ar) und e+’* sind, so daB 
wir fiir ¢ bzw. s die partikularen Lésungen 


= (pe **+ ge**) Jo(Ar) (21) 


bzw. ; 
—iz he 
shen Ree (22) 
bekommen; hierbei bedeuten p, q, K,, K, willkiirliche Konstanten und Jo(Ar) und J,(Ar) die 


Besselschen Funktionen erster Art. Das Bestehen der Gleichung (18) ergibt wegen der bekannten 
Beziehung 


dJ,(Ar , 1 
Gy i (Ar) =—— J, (Ar) + Ag (Ar), (23) 
dai 
_ 20—») . _ 2(1—») 
K\= ih =e ge Resi sere 7 


sein muB, und man zeigt leicht, das die an Stelle von (22) tretende Lésung 


7A ioe — Az 2 
am ae e— ge") Jy (Ar) (24) 


mit (21) der Gleichung (17) ebenfalls geniigt. Jetzt kénnen wir zur Bestimmung der Verschie- 
bungen iibergehen. | 
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Aus (12) folgt mit (21) die inhomogene Differentialgleichung 


Ab = — >" (—pe** + ge") Jy(Ar) (25) 


fiir €¢. Der homogene Teil 4¢=0 hat — entsprechend (19) und (21) — mit willkiirlichem 
P und Q die Lésung 


€= (Pe **+ Qe") J,(ar), 
wahrend der Ansatz 
Be Ae Ard oe Tar) Ace B 


fiir eine partikulare Lésung von (25) mit den neuen Konstanten A,, A,, 4 und B wegen J’ (Ar) 


= — J,(Ar) auf 
E Ae Mane diet, ds) ar) — pat ge), 
fiihrt, woraus 

Meee Pi eee en ts q 

WPT CTU Cy ak Oar MPT OREOE 


und somit nach Superposition der beiden Lésungen 


C= — say Pe + ae) Jo(Ar) + (Pe "+ Qe) I,(ar) + dz+ B (26) 


folgt. 
Die Ermittlung von @ gestaltet sich etwas miihsamer. Unter Beriicksichtigung von (15), 
(9), (18) und (24) erhalt man aus (13) zunachst 
Q I Os h —hz hz 
40 =a = eee Fae) Jar) (27) 
wieder eine inhomogene Differentialgleichung fiir 9. Die homogene Gleichung Ao — o/r? = 0 
hat mit den beliebigen Konstanten B, und B, analog zu (20) die Liésung 


Q EG oy ee gee B,e"*) Ji (Ar), (28) 
wahrend der Ansatz 
9 = (De “* + Dye") z J, (Ar) (29) 
in (27) gesetzt nach einigen Zwischenrechnungen zu 
ae 1 ee isa ehh) aS. 
ei 2(1 — 20)’ eee, 


fiihrt. Nach Einsetzen von D, und D, und Superposition der beiden Liésungen (28) und (29) 
entsteht 
: —hz hx —he Az . 
ast erases ae ‘ ge 4) Jy (nyse Bre + B,e**) J,(Ar). (30) 
Die Konstanten B, und B, lassen sich auf p und q bzw. P und Q zurickfiihren, wenn man 
beachtet, daB 9 so gebaut sein mu, daB die Gleichungen (9) und (15) bestehen: Aus (15) ergibt 
jsich mit (24), (26) und (30) 


(3 — 4») p (3 —4 9) q 
==, FoR ‘ se | a 
i i 2(1—2)A ’ Bs Om 2(1— 22’ 
bund man kann zeigen, dafB die mit diesen Werten aus (30) hervorgehende Lésung 
Sa ir, Me I 3—4 ay 
C= 59-95) 2 — qe"*) J, (Ar) |? a Oa Oeil | 


(31) 


(3 —4y eh 
+ | O+ 3a—2¢ nif? | Aa(2n) | 
‘zusammen mit € nach (26) der Gleichung (9) ebenfalls geniigt. 

_Damit sind gemaB (26) und (31) gewisse (spezielle) Lésungen gefunden, von denen aber von 
-vornherein noch keinesfalls erwartet werden kann, da®B man die in ihnen vorkommenden Kon- 
stanten p,q, P,Q, A, B und 1 so wahlen kann, daB die dadurch entstehenden Lésungen allen 
-Randbedingungen perecie werden; zur Klarung dieser Frage gehen wir jetzt dazu iiber, unter 
Heranziehung der Randbedingungen fiir die Platte, die Konstanten bzw. einige von ihnen zu 


bestimmen. 
i 
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| 
| 
i 
| 
| 


Die Forderung (7) (verschwindende Schubspannung auf dem freien Plattenmantel) fiihrt mit 
Hilfe von (26) und (31) aus (10) auf J,(Aa) = 0; also die Eigenwerte As Ay (k=, 2; 3.) sie 
durch die unendlich vielen, reellen und positiven Nullstellen von J,(t)—= 0 bestimmt: Aus} 


Dit es Ortolety A, = tk  wobei die t, z.B. aus dem Tafelwerk Jahnke-Emde entnommen werden) 
a 


kénnen. Durch Superposition der so entstehenden EKigenfunktionen treten an Stelle von (26) 

und (31) die Lésungen 

€ =— 5 (pe +e) Io(Aar) + > (Pee + Qe) So(Aer) + Az+B (32) 
k=1 | 


2(1 — 2») ao 


und 


z = A pag > 
See A fj, eM) J, (A, 1) > | Pet saan ° Ay, +| 


Q,, x =e ral on Jy (A r) . 

(Der in (33) auftretende physikalische Widerspruch, da8 namlich @ (a, z) =0 ist, wird noch eine | 
Erklarung dahingehend finden, da® diese Funktionen fiir €¢ und @ nicht mehr genug frei wahl-| 
bare Blomence fever um allen noch nicht erfiillten Randbedingungen gerecht zu werden; d. he 
(33) ist noch nicht die ,,richtige*‘ Liésung fiir 0. ) 
Mit (32) und (33) kénnen wir weitere Randbedingungen der Platte realisieren. Die Last- \ 


bedingung (1) liefert mit (10) 
26(5 to ee 
woraus mit (9), (32) und (33) 
S : Lew SEG) 
folgt. Multiplizieren wir diese Gleichung mit J)(A,1r)rdr (n= 1, 2,3...) und integrieren von 


r= 0 bis r= a, so erhalten wir wegen Ji (Ay a) = 0 und der gemaB dieser Gleichung giltigen 
Formeln 


(33) f 


a a ee 0 fiir k = Tt, | 
[To(Anr)rdr—= 0, | JoAnn) Jo Cen) rdrestne (34) 
0 0 | 5 Je (A; a) fiir cn. J 
— | f(r) Jo(A,r)rdr 
P, di jf pelted J (r) Jo (Agr) fo (35) 
D 9 ek ai “b k a2 GJ? (A;,@) ko 
wahrend sie mit rdr multipliziert und integriert 
Ym 2S : 
Aa [ forar (36) 
f 
liefert, womit A bekannt ist. 
Berechnen wir mit (32) und (33) gema® (10) die Schubspannung t = Tf (z, r), so bekommen | 


wir aus der Randbedingung (2) 


a0 


a > Ay (Py + Q;.) Ji Co sim Bes (Pe = Gr) Jy (Axr) aia 


Wiss (3 — 47) p, (3 — 4») q 
ios ay | Pe zm | ! a ee rae ae) = 0 ut 
woraus 
: ae 
Ay (Pit Qe) py, (Pe — ie) = 0 (37) 
folgt. Die Auflésung von (35) und (37) liefert 
(3 — 47) p, I 
Ia a =r sen aes Gg) 
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Pi _ B= 4) 4, 
2 Sty) BFS 2) 


2 i Qk = 


— Ch» (39) 


wobei die c, nach (35) definiert und somit durch die Belastung gegeben sind. 
Auf dem gleichen Wege liefert die Randbedingung (3) 


22 —Ayh Ak . a : 
am lee ge) Ah Si ar)— > (Pr2me + Q, 24,0) J (der) — 


k=1 ; k=1 


2(1—») ~ —Ayh Aygh\ 
— TF De ae) b= 0, 


woraus unter Benutzung von (38) und (39) 
pr(Ahe™ +. Sin Ah) + qn (Aphe** + Gin A,h)=2 (1 — 2) o Gin Ah (40) 


_hervorgeht. Damit kénnen aber alle in (32) und (33) vorkommenden Konstanten mit Ausnahme 
von B auf eine, etwa p,, zuriickgefiihrt werden, d.h. zur Beschreibung des Spannungs- und 
Formanderungszustandes der Platte fehlen uns noch die Koeffizienten p, und die Konstante B, 
Zu ibrer Bestimmung miissen wir die noch nicht realisierten Randbedingungen heranziehen und 
dazu benétigen wir die Verschiebungen der Unterlage. 


4, Die Verschiebungen in der Unterlage (Halbraum), Da fiir die Verschiebungen und Span- 
-nungen des Halbraumes dieselben Differentialgleichungen gelten wie fiir die Platte, sind die 
Verschiebungen ¢, und g, fiir die Unterlage unter Beachtung von (32) und (33) sofort hinzu- 
schreiben: An Stelle der Summe treten im Sinne des Laplaceschen Gedankens Integrale, wahrend 
die Konstanten durch unbekannte Funktionen zu ersetzen sind; bedenken wir noch, dab fiir 
2, —> co sowohl ¢, wie 0, verschwinden miissen, so haben wir 


co fee) 


= — aq tay | PM Sone art [KAI (ane al, (41) 
0 0 
ie f 3— 42 \D A) Pmtics" 
= sq O(a) J, (Arye dA + i wa) Rtas 2D di, (42) 
0 0 


-wobei (A) und ¥(/) aus den Randbedingungen zu ermittelnde Funktionen sind. Die Rand- 
bedingung (4) ergibt 


f paw) oF sa ee P(A ) J, (Ar) di = 0, 


woraus die Beziehung 


(=e Dan, 
P(A) + él =e) POW eae v (43) 


zwischen Y(A) und (A) folgt. 


5. Die Ubergangsbedingungen zwischen Platte und Unterlage und die Aufstellung eines un- 
endlichen Gleichungssystems fiir die Koeffizienten p,. Die Ubergangsbedingung (6) (Gleichheit 
der Normalspannungen) lautet mit Riicksicht auf (10) 


| Vy ey = OG. Be 
ag Ge ETA 26 (5: | ee 


wobei ¢ bzw. €, durch (9) definiert sind, so da®B wir durch Heranziehung von (32) und (33) (mit 
Benutzung von (38) und (39)) bzw. (41) und (42) [mit Benutzung von (43)] mit der Abkiirzung 


= (44) 


a 
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nach laingeren Rechnungen wegen o(h,r)= 0 fiir r > a zu dem Ergebnis 


[[PA+ 9, PO] Solar) das sug | MAAS (Ar) da = 
g 0 
Oar re 
zi, “= : 5 45 
| ae | De I) (A, he~“#* — Gin Ah) — (45) 


k 


a yey (A, he” — Gin Ah) + e, Co} uh] Jo (Ann) | fix r<a 


gelangen. Multiplizieren wir diese Gleichung mit J,(#r)rdr und integrieren von r = 0 bis r = 00 
(bzw. rechts von r= 0 bis r= a), so bekommen wir unter Benutzung der Hankelschen Umkehr- - 
formel 


| rJy(0r) AA pe TO f J, (Or) F(8) dO = g(r) 


und der Integralbezichungen 


a 


| TolOryrdr = Na), | Toler) To(Dr) dr P20 eee) 
é ae Ie 


0 


wenn wir nach Ausfiihrung der Integration fiir wieder A schreiben 


Pay Jaden _ (1—») Aa J, (Aa) 
ee i faah va > One ee ee (46) 


wobei 


A,he “*" — Gin A,h) — 4 


2 (I say (4s he’*” _ Gin Ayh) + c4Cof A,h 


oe Pk 
ua (1 — 2») ( 
oder mit (40) 
(A; h)? — Gin? A, h Lo, wht e*k* Sin Ah 
k le 
(ao) ane Gin d,h] | "A,he** + Sin a,h 


bedeutet, also lediglich von p, abhangt. 
Die Gleichheit der senkrechten Verschiebungen (Randbedingung (5)) liefert die fiir r<a 
giltige Gleichung 


le - 
> aso 4y+2A,h)e er oe 


C, = p (47) 


q Be | | 
i a (3 —4y—2 Ane + "| 1 26, Sin a ete ae (48) | 


ie ARE ale Ov aya 


die mit r dr multipliziert und zwischen r= 0 und r= a integriert zu der Beziehung 


fv.) ) F(a) A = * (4h + B) (49) 


fiihrt. Multiplizieren wir (46) mit J,(Aa) d//A und integrieren von r—0 bis r= oo, so ergibt 


sich wegen (40) 
> Gh wa ty, ae [ve i(Aa) a). (50) 


6 
wahrend die Gleichung (48) mit Jy (A,r) rdr multipliziert und von r= 0 bis r—a integriert 


Ah-+ B=—4y(1— ») 


J WHA) J Solar) Jy(Anryrdr daa | (ay Pon) SAO) gy a8 Ona) py 
4=0 0 0 ne 4) 


Ay 
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oder 
caprovel ANU W9) a J,(2a) 
[va pee ig. Vo (51) 


liefert, wobei (nach Zuriickfiihrung aller Koeffizienten auf Pn) 


4(1 —») lens 
ices | | a (7. h+ > Gin 2 Anh) + c,e*® Sin Aah (52) 


A,he’n® 1 Gin d,h 


ist, also bis auf p, ebenfalls als bekannt anzusehen ist. 


Und zum Schlu8 multiplizierten wir (46) mit ce dj und integrieren von A = 0 bis A= co 
und bekommen unter Beriicksichtigung von (51) : 
Jona) py hikes! r Bie da (I=) A f JR(Aa) dA 
4H, "1 > ae 2) (Ba) | 1— 20 P72 
0 0 


und nach Einsetzen von (47) und (52) 


(1—y») Jo(4,,4) Pp 2 1 me 
Be h | 2) 9 » Anh 
A, (A,he*n® + Gin A,h) | 1— 39 ( nh 3 S24, h)+ c,e2""” Sind, h| = 


co 


| 
ky eo Gey > [ps (1,h)? — Sin? Ah Lo, et pe : | 
| 
J 


1G hd — 2) ayhe®* + Sindh) | “Ake + Sin Ah 


| : 2? J? (Aa) di ea Pca I5 1) 
oe aot (#28) (#— at) © 12» jas an 
6 6 

Non stellt diese Beziehung fiir n = 1, 2, 3, .. . ein lineares, unendliches Gleichungssystem fiir die 
unbekannten Koeffizienten p, dar. (Auf die uncigentlichen Integrale in (53), die sich als konver- 
gent, sogar geschlossen darstellbar erweisen werden, kommen wir noch zuriick.) Die Auflésbar- 
keit dieses Gleichungssystems vorausgesetzt (die Frage der Konvergenz wird in einer Sonder- 
betrachtung im bejahenden Sinne beantwortet), kommen wir zu dem Schlu8, daB nach Ermitt- 
tung der p, riicklaufend alle anderen Koeffizienten (die a aus (40), Baus (50), P, und Q, aus (38) 
und (39)) und unbekannte Funktionen ('¥(A) aus (46), ® (A) aus (43)) bestimmt werden kénnen, 
wodurch Verformungs- und Spannungszustand in der Platte und in der Unterlage bekannt sind; 
und das bedeutet wiederum, da wir keine Méglichkeit mehr haben, die noch offene Rant 
bedingung (8) (Verschwinden der Radialspannung am Plattenmantel) zu erfiillen. Die Erklarung 
hierfiir ist: Die durch den Produktansatz gewonnenen Liésungen der Form e © 47 J)(Ar), e* 47 J,(Ar) 
baw. ze+4* J, (Ar) und ze+7’* J,(Ar) und insbesondere die Wahl des Eigenwertes sind nicht all- 
zemein genug, um alle Randbedingungen zu erfiillen. AnschlieBend zeigen wir, wie eine ,,Kor- 
rektur“ fiir die Verletzung dieser Randbedingung im Sinne des de Saint- Venantschen Prinzips an- 
gebracht werden kann. 


6, Radialspannungen und Radialmoment am Plattenmantel, Da die Radialspannung am 
Plattenmantel o, (a, z) nicht fiir alle z verschwindet, wird die ,,resultierende Spannung“ 


Ql 
>| 
——s 
a 
e 
ee 
QQ 
we 


| 0 

und auch das auf die Langeneinheit bezogene Radialmoment 
h 

| mr = | 0,(4, z) 2dz 

| 0 

von Null verschieden sein. Aus (10) erhalten wir mit (32) und (33) 


> & dy (pre “#*— que’**) ae iG (Pye “ ee Q,e"**) ze 


6. (4,2) == 26 Pe a1 — 2») 


3—4 vA 


ry (54) 
v (ppe *#* que’**) TROD ; | ; | 
2(1 — 2») 1— J | 


| 
| 
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und somit ergibt sich fiir o, und mg, indem wir vermége (38), (39) und (40) alle Koetiionte 
durch p, ausdriicken, nach recht mithsamen Zwischenrechnungen 


on = oe y) (55) 
2 2 = 61h") (Ah = Crt ayh) 
2 | Pp(Aph? — Stn? pe) c (Ue, eth ly A; | 
le DoE 20) he (A,he*® + Gin d,h) | ABA? (A,he™* 4+ Gin A,h) oes sl 
be } ( ( ) 
Liat — 29){ ; | 


Um die Randbedingung am Plattenmantel beziiglich der Resultierenden zu erfiillen, belasten 
wir den Rand durch die gleichmaBig verteilte Spannung —o, auf Zug oder Druck und durch; 


a 2 
das Kraftepaar — mp = — mp -- uae auf reine Biegung, wodurch die Verschiebungen + 
he a Pcs Le) Gy 
eS anela ‘). Serie 
h\? 57) 
ma |2o[2— 3] +a—n-] 7 ~67ig(l—»)(2— 5) ‘ 
ote (1 +») Gh ey Reet +») Gh j 
und somit nach (10) die Spannungen 
C08 Ca C.= 0, == <6 | 
i yt 1 ut 12 Mr h (58) 
G, = 0, a == CO, == 0; = hs (: =) | 


hervorgerufen werden. Auf diese Weise wird man der Randbedingung (8) nach de Saint-Venan 
gerecht, aber durch diese ,,zusatzlichen Belastungen** wird das Spannungsfeld in der Platte und| 
wegen des Kontaktes zwischen Platte und Unterlage auch das des Halbraumes geandert und es 
ist jetzt zu priifen, wieweit die Randbedingungen bzw. die aus ihnen flieBenden Gleichungen zx 
modifizieren sind. 


7. Die durch achsensymmetrische Lasten, reine Biegung und zur Achse orthogonale (gleich- 
mafige) Randkrafte beanspruchte Platte. Der neue yecenaary ee ist durch die Super-' 
position der Verschiebungen C, @ bzw. ¢’, 0’ und €’’, 0” gegeben: 


Q 


2b bo anc. pee Mie 
Da nach (58) die von —o, und — mg herriihrenden Normal- und Schubspannungen identisch ver- 
schwinden, ist es offensichtlich, da die aus den Randbedingungen (1), (2) und (3) hervorgehender 
Gleichungen (36), (38), (39) und (40) bestehen bleiben; dasselbe gilt fiir die aus (4) und (6) flieBen- 


den Gleichungen (43) und (46). Eine Anderung wird dagegen die aus (5) entstandene Gleichung 
(48) erfahren: Aus 


Z(h, r) = C(h, v) + "(hy r) + 0°"(h, r) = ¢, (0, 7) 
tritt an Stelle von (48) die Gleichung 


2 yer: 1! Jo A r | 
oe (3 4 v 2 Ah) ehh Lo EB Sin a ot 4 Ah | B - (59) 
Gh 3mp [vh? +2(1—v) 27] faz 
a 2(1 + ») Gh =| P(A) Sy (Ar)da 


6 
Mit dieser Gleichung verfahren wir genau so wie mit (48): Zuerst multiplizieren wir sie mit rdr' 
und integrieren von r= 0 bis r= a und erhalten 


fieac ) I(aa) S = Ns Lepo rh _ 3 mp [yh + (1 —») a?) (60) 


FCO) 4 2(1 + ») Gh3 (? 


ey) 
sl 


1 A, Clebsch, Elastizitat fester Kérper (1862), S. 181. 
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womit (50) durch 


orh 3m [v h? + (1 —7)a?] 

eet 2G(1—») 2(1 + v) Gh3 
os. anes eI (adi (=) 4 ¢ Jala) di (61) 
=—4y (1 pS CoTo(des) | pag + Tay | ae | | 


zu ersetzen ist. Multiplizieren wir (59) mit rJ, (A,r) dr und integrieren zwischen r = 0 und r= a, 
so bekommen wir wegen 

Me ae 

[PJ nr) dr = Fe Jo(Ana) 
0 n 
anstatt (51) 


Pw (ay 2h Ga) gg aJo(2na) py | 6 mn (1—>) a 
J poe di, 4 An D, + (lL -+ »)G RB 22’ (62) 


»wodurch an Stelle von (53) nach Einfiihrung von (55) und (56) fiir die p, das Gleichungssystem 


(1 —»v Jo) (Ana) Anh cs Pn in 
n } | 2 

Van (nhe'n* + Gin Anh) ewe Se Ta (2 Pg et 2 uh) 

a 12 (1 — 7) ll p, (Az h? — Sin? A, h) plies oe) (A, — Sin A,h) TO nee 
Anh (L +») | (1 — 2) agh?(A,he“k*+ Gind,h) AB R2(A,he“e* + Gin d,h) |~O 7 
= (2h? — Gin? A,h 4 bCin Ah 
fe 2y a—»)> Pr ; ae 47) ! eum e Sin A, . (63) 

i | d—2) (A he" + Sin A,h) A,hetk™ 4+ Gin Ah 
; f ABB (Aa) dA 24 (1—»,)(1—»)A f[ JB (Aa) dd 
x Jo (A; @) | (2? — 22) (A? — 22) earn haa Dy Ae Ae J 
0 0 
» hervorgeht. 


) Fiir die rechnerische Handhabung dieses Gleichungssystems ist, vom Standpunkt der Mathe- 
»matik, zunachst die Existenz der in (50) und (63) vorkommenden uneigentlichen Integrale und 
| dann die Existenz der Liésungen p,, py,... zu zeigen +. Wir wenden uns jetzt der ersten Frage zu. 


I(t) dt fait) at 
HUES In= | 


9 2 Peed 
t? —- t 


8, Berechnung der Integrale: J = f 
0 


at) ) cee 
Tin= J (@—t) (P—f) Wl bey Bsn 


[Durch die Substitution Aa = t lassen sich die in (63) bzw. (50) vorkommenden uneigentlichen 
{ Integrale auf diese zuriickfiihren: 


me (Je da) dh f  j2J2(Aa) da 
[22 aa=at, [08 = ol, [ pes = tl 64) 
0 0 0 


— 72) (2 — 72) 


' Hierbei sind t, bzw. t, Nullstellen von J,(t), d.h. es ist J, (¢,) = J, (t,) = 9. 
| Um die Auswertung von I brauchen wir uns nicht zu bemiihen. Es ist ein Sonderfall des sog. 
» Schafheitlinschen Integrals”, und es gilt 


[Roe e (65) 
0 


4, 
i east 
Die Bestimmung von I, erfordert mehr Aufwand an Uberlegung und Rechnung. Durch Partial- 
bruchzerlegung und vermége der Substitution t—t,— x bzw. t-+-t,= x erhalten wir zunachst 


(sid 1 f 1 i re eal eae eas PIa—t,) 
z,— | 2 =z,{ Hol )at a! erp J ie 


?—t,? Zt tb 
ib 


. (66) 


_ 1 Bei der Klaérung dieser Fragen habe ich meinem verehrten Kollegen, Herrn Schmeidler, fiir wertvolle 
} Unterhaltungen zu danken; insbesondere riihrt der in Ziff. 9 gegebene Konvergenzbeweis und die Naherungs- 
-formel (76) von ihm her. 

| 2 Watson, Theory of Bessel Functions, 5. 403. 
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Unter Beriicksichtigung der Additionstheoreme der Besselfunktionen+ und Beachtung vo 
Jy(t:) = 0 und J_,(x) = (—1)* J,(x) — fir a > 0, ganz — ergeben sich 


Tet) = S (Ig (2) (Sana Oo) — Ines > | 


K=1 


(67 
Hla = SI) [Ina — Jers Ol] | 


und damit 


FIR (a+ : P Jo (x) Ja(x) dee 
Ay = f BEE) gp — ST (18 +8 LS, 3s) — Fas 1(60)] [pa ts) — Jp galt] | 22 
—! Ko ce 
und 
:— ( Ja(x) Ja (x) de 
iti jes Ne = DS) a-att) Su salt)] Yp-al)~ Inatoy] J EP | 
tn 
Wegen a = a i und ia ey, bekommen wir mit 
8 : Oe ] 
Jo(%) Jp (x) Fp ie | Ja G8) Jp) C= ijn PORE 
ak i, ; i, 
A;— Ay= 1)*F? — [Tua (ta) — Fa (ta) |[Jp—1 (te) — Tpa(ta)] x | 
a, p=1 
x des ede —2 DS) [Funalts) = Snsr t)|lJo-a(6) = Jonat6]% 5 (6g) 
ee grea 
«| Jo.(x) ECBHIC) 


Andererseits ist 2 


(2) 


; Tt 
f net a. 2 sin (« — B) = 
Ld (ted) 


0 
und da im vorliegenden Falle a -+ f = 2 m- 1 (ungerade) ist, erhalten wir 


a—B-1 
2 


sin (a — f) > = sin (2 m+ 1—2f) 5 = cos (m—f)n=(—1)"?=(—1) 7 | 


und somit nach Einsetzen von (68) in (66) 


he, ; ‘ a 
=| peta a 2 tg (Fa-1 (te) Jaa ta)]LSp—1 (bn) — Spar lta) (69) 
a+B rai ' 


Die Kenntnis von I, erméglicht es, vermége einer identischen Umformung I,,, ohne jede weitere | 
Rechnung ldnmechreneal 


9 ? 2 2 
ae PERO GE age F(t) tr UR dt 
ee Pee R) BS ee rae lene 
6 0 


und das bedeutet, daB 


ty = t2 f,, (70) 


 Jahnke-Emde, Funktionentafeln (1933), S. 213. 
2 Watson, Theory of Bessel Functions, S. 404. 
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ist, wobei I, bzw. I, durch (69) gegeben sind. Der Fall k= n muf gesondert behandelt werden. 
Setzen wir I, und I, nach (69) in (70) ein, so ergibt sich zunachst 


a—B-1 
| 2 =< Gye 
he oe pay (8? — a®) ( — 2) {te [Fn —1 (te) = Joos (te) [Ja—1 (te) — Spr (te)] — 
ie o+p=2m+1 


ae tn tit) a J 44.1 (ta) ] [Je—1 (tn) = Jp i(t,)]} o 


Fiir k = n erscheint diese Gleichung in der unbestimmten Form 0/0; die Bernoullische Regel an- 


; 5 9 " 
gewendet, erhalten wir unter Heranziehung von Funktionsgleichungen fiir Besselsche Funk- 
tionen nach langerer Zwischenrechnung 


2 
se 4 a—p—1 (2) a= | 
BEG) ais 1 4 = ae NGS 
Ti, == | (2 — Rp aay fi, iz pa () : rer J (ty) [Jp—1 (te) -— Jp4i(t)]. (71) 
0 Palas 


a&+B=2m+1 


Damit sind die fraglichen uneigentlichen Integrale numerisch auswertbar. Fiir k, n= 1, 2, 3,4 
sind die Zahlenwerte!: 


n 1 2 3 | 4 
es —0,04455282 —0,01707376 —0,00917584 —0,00583368 
Lh —0.08610633 —0,00539201 —0,00332785 —(0,00234266 
Ths —0,00539201 —0,05418659 —0,00201452 | —0,00152184 
Ie —0,00332785 —0,00201452 —0,03905402 | —0,00116055 
{Bar —(0),00234266 —0,00152184 —0,00116055 | —0,03013391 


9. Beweis der Auflésbarkeit des unendlichen Gleichungssystems, Mit Riicksicht auf (64) 
erscheint das Gleichungssystem (63) in folgender Form: 
(1 —») Jo Gn@) 
An Unhe’n™ + Sin Ah) 


« ooh Gin Ah — Pe (Jai ++ Gin2 Aah + 


i W(il—r) < Py (zh? — Sin? A,h) | (1 —e*k") (1h —Gind,h) reaea 
SS + ¢ 7 k 2) = 
Ah (1 + 2) ee (1—2v) Azh2(A,he“k* + Gin A,h) |” ABR? (Ayheée* + Sin Aph) | ~° (79) 
* } 
Sy) py, (Azh? — Gin? A, h) Tes pae etkh Sin Ah | 
a 2 1 (1 — e 2 = SS J, (A, a) aL», — 
a aS aaa (ahete® + Sindh) ' "7 hete* + Sin Ah 0 (Aca) 


ALAN Den eat ay pe | 
I —2y 6 


jUm die Existenz einer Lisung p,, po... . fiir dieses Gleichungssystems nachzuweisen, dividieren 
wir (72) durch 


1 
1—v) Jo (Ano) Anh + = Gin 2 1,h) 
eee. : (73) 
(1—2r)A, (A,he4n"+ Gin A,,h) 


Dann wird das erste Glied der linken Seite 


Anh ei 
rn” Sind, h 
te eee Da ciel 2) po (1 =P) e a a 
(A,h + 5 Sin 24,h) A,h + 5Gin 2Ayh 
ind somit geht nach Einfithrung von 
Pr=2¢,(1—27)+ 4, (74) 


a1 Bei diesen und den anderen numerischen Rechnungen haben mich meine AY die Herree 
Jr. Hoffmann und Dipl.-Ing. Reckling wirksam unterstiitzt; ersterer hatte auch die oft langwierigen 
ligemeinen Rechnungen kontrolliert. 
10 
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(72) wber in 
Oy Akh? — Gin? Ah) | 


baa aul ae) ce 
e oie Da see nT Wee tae 
Ah + 5 Gin 2A, h ” fai CL 2 0) AG be et Sine 


— 6,+(1— 2 r)¢, 


| 2 (Azh? — Gin? J,h) (1 —e*k") (1,h—Gind,h) \y oe 

L_¢ fe Wels I — 

* | jen? (A,he'k” + Sindh) | A2h2(A,he*e* + Sindh) || °° 

died Bt) SS cL ra (75)} 
wT: N pat | (1—2p) (A, he**h + Gin A, h) 


Jo (Ax) AT pn — | 


| Aah beteh Sin yh ee 
“lA, hen’ + Gind,h | A,he“k*+ Gina,h |] 
2 u(1—y»,) =P) 

Se Ae 

Mit wachsendem n(n—> oc) nehmen auch die Higenwerte 7, zu (A,—> oo), so daB das zweite> 
Glied auf der linken Seite fiir n—> oo verschwindet, da die eckige Klammer gegen Null strebt. | 
In den weiteren Gliedern streben die eckigen Klammern mit wachsendem k ebenfalls gegen Null! 
und dasselbe gilt fiir das letzte Glied der rechten Seite (da fiir n—> co nach (73) N ier alle; 
Grenzen wichst). Wegen des mit n wachsenden Nenners N kénnen daher fiir groBes n alle Glieder| 
mit ¢, und ¢, wie auch das letzte Glied fortgelassen werden, so da ein homogenes Gleichungs- ! 
system iibrigbleibt, dessen Glieder fiir k= 1,2,... n— 1 wegen des grofBen N ebenfalls un-; 
endlich klein sind. Das nun iibrigbleibende Gileichungssystem fiir 6,,6,,,,... kann durch} 


O¢ = Ony1 = +++ =O geldst werden. Damit ist nicht nur die Lésbarkeit des Gleichungssystems 
(72) bewiesen, sondern wegen (74) zugleich gezeigt, daB fiir groBe n(> 4) die Naherungsformel 
Pa 2 ¢,(1— 2 v) (76) 


gilt. In den meisten Fallen wird es also geniigen, p,, Py, P3, Pa aus dem fiir n = 1, 2, 3, 4 aus (72) 
flieBenden linearen Gleichungssystem zu berechnen und die weiteren p, (wenn ndtig) nach (76) 
zu ermitteln. Die Herleitung eines solchen Gleichungssystems fiir die ersten vier Koeffizienten 
erfolgt in der bei einem unendlichen Gleichungssystem iiblichen Weise: Man setzt in (73) n= 
1, 2, 3, 4 und geht jedesmal in den Summen mit k ebenfalls bis 4, wodurch man zu einem linearen 
inhomogenen Gleichungssystem kommt. 


10. Ein Beispiel. Eine mehr den Bediirfnissen der Praxis angemessene numerische Aus- 
schépfung der hier gewonnenen allgemeinen Ergebnisse, insbesondere Vergleiche mit den Ergeb- 
nissen der auf der Kirchhoffschen Plattentheorie und der Zimmermannschen Hypothese auf- 
gebauten Untersuchungen sollen an einer anderen Stelle erfolgen; hier sei lediglich ein kurzes 
Beispiel angefiihrt. 

Fundamentplatte aus Stahlbeton: a=6m; h=2m; G=15:104kg/em™?; y= 1/6. | 
Die Last von P= 500t sei innerhalb eines Kreises vom Halbmesser r, = 2,5 m gleichmahig 
verteilt (Abb. 3). Fir die Unterlage (Mittelsand) nehmen wir an G, = 413 kg/em~? und », = 0,45. 


Nach (35) erhalten wir mit o) = —— 


Tp ot 
To 
ere eC Lee i OT) Jy(Ay To) 
“b= a GF3 (Aya) | Jol ie i,a® G J3(4,a) ° (77) 


0 


Fir diesen Belastungsfall ergibt die numerische Rechnung fiir die ersten vier Koeffizienten 
Py = — 22,3272-10-*; p, = —8,8040-10-*; ps = 2,7557-10-*; p, = 4,9260-10-*, wahrend die 
Naherungsformel (76) py = 5,019-10~®liefert. Mit Hilfe dieser Koeffizienten bzw. der zugehérigen 
Kigenfunktionen sind die wegen der nicht streng erfiillten Mantelbedingung in erster Linie inter- 
essierenden Radialspannungen am Plattenumfang und in der Mitte berechnet und in Abb. 3 mit 
und ohne Randkorrektur aufgetragen. Man ersieht aus diesem Bild, daB die Randbedingung am 
Mantel (0, (a, z) = 0) in ausreichender Weise erfiillt ist. Die maximale Druckspannung zwischen 
Platte und Unterlage ergab sich fiir r— 0 und z= h zu 0,70 kg/em~? und die Einsenkung der- 
selben Stelle betragt 0,32 cm. 


>a 
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11. Eine verbesserte Niherungsformel fiir die Koeffizienten Px. Die Approximation (76) stellt 
}natiirlich fiir die ersten Koeffizienten noch eine recht rohe Naherung dar, aber ihre Kenntnis 
verhilft uns zu einer SchluBweise, die auf eine bessere Naherungsformel fiihrt. 

Das vorangehend angefiihrte Beispiel zeigte bei der numerischen Behandlung, daB die Auf- 
lésungen der Gleichungen (53) ohne Randkorrektur bzw. (63) mit Randkorrektur sich nicht 
wesentlich unterscheiden, so da® wir zur Herleitung einer Naherungsformel von den ohne die 
/Randkorrektur gewonnenen Gleichungen ausgehen kénnen. Wir kénnen (48) mit 
Jo (Ax) 


2 Mi, 


on Ah+ B= g(r) (78) 
k=1 


in der Form 


2 


| | 
20-40 0 40 20kG/ew 
WZ 


Abb. 3. Platte mit gleichmiBiger Last. Radialspannungen am Rand und in der Mitte. 


— — — — —_ — ohne Randkorrektur, mit Randkorrektur. 


ischreiben. Dehnen wir die Giiltigkeit dieser nur fiir r < a bestehenden Beziehung fiir alle r aus, 
so erhalten wir nach der Hankelschen Umkehrformel 


=e) uss Pane) rdr= J Jo(Ar) g(r) rdr+ igor g(r) rdr 


oder mit 
[ Jo(An) g(r) rdr = y(2) (79) 
MO) LJg(Ar) ale) rdr + yl). (80) 
fNach Einfithrung von (78) und Ausfiihrung der Integration ergibt sich aus (80) 
we = (Ah+ B)* J, (Aa) Be Mi = D, + yd), 


wworaus fiir A= J, (da auf der rechten Seite nur k = n einen Beitrag liefert) 


Pan) a2 Dy, | = 
Go) Pe 54,0) + yn) (81) 

wnervorgeht. Auf demselben Wege bekommen wir aus (46) fiir A= A, 
VY (An) = — w(1 — »,) a2 J2 (A, 4) C, , (82) 


10* 


| 
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und aus (81) und (82) flieBt 
4 An Y (An) | 
D, + 4u(1— »,) G, = — | 


a Ji (Ana) 


Setzen wir hier (47) und (52) ein, so erhalten wir 


(1 —»)e*2" Cin hs — o,) M(A,h + e An® Gin An h) | | 
—>) (2A,h + Sin 2/,h) — 2 w (1 —»,) (2h? — Sin? 7h) (83)) 


1g see red yg 


A Hey Onh + Gin Anh) ) p(An) 


| ; 
 @ J2(A,a) (1 —r) 24, h + Sin 2/,h) — 2 w(1 — »,) 02h? — Sin? 7,h) | | 
Da das erste Glied in (83) fiir n—>co gegen 2 (1— 2¥)c, strebt und im zweiten der Faktor’ 


von w(A,) gegen Null, folgern wir mit Riicksicht auf (76), daB 
(1—»)e4n® Gin Ah + (lL —»,) M(Ayh + e’m* Gin 1h) 
"GW —v) Q4,h + Gin2/,h) —2 (1 —»,) w(2 hk? — Gin? A, h) 


Pa 2(1— 2) ¢ (84)) 


eine Naherungsformel fiir die Koeffizienten p, bzw. p, darstellt. Fiir n—> co geht aus (84) wie-} 
der (76) hervor. | 


Aus (63) _ | nach (84) | macni(de) Die Gitte der Formel (84) zeigt die nebenstehends} 
p,: 106 | —22,3272 | —27,840 | —8,643 Tabelle, in der die aus (63) ermittelten bzw. nach (76) | 
po LOe | == 8.8040) 297,275) 5,482 und (84) berechneten ersten vier Koeffizienten fiir den} 
P3° 10° 2,7557 2,400 2,264 in Ziff.10 behandelten Belastungsfall nebeneinander-; 
Daclos 4,9260 5,082 5,019 


gestellt sind. | 


12, Zusammenfassung, Behandelt wurde das Problem der auf dem Halbraum aufliegenden, , 
achsensymmetrisch belasteten dicken Kreisplatte. Unter Zugrundelegung der Hockeschen Ge- 
setze und der Annahme des stetigen Kontaktes bei vollkommen glatter Beriithrungsflache wut 
den fiir den Spannungs- und Verformungszustand Lésungen hergeleitet, die mit Ausnahme dea 

| 


Radialspannung am Plattenmantel simtlichen Randbedingungen geniigen. Dieser Randbedin- 
gung wurde nach dem de Saint- Venantschen Prinzip in der Weise gerecht, indem das aus den Rest- | 
spannungen hervorgehende Moment und die Einzelkraft am Plattenmantel im umgekehrten 
Sinne als Last angebracht wurde. Die Entwicklungskoeffizienten der Lésungen fiir die Platte er- 
geben sich aus einem unendlichen, linearen Gleichungssystem, dessen Aufliésbarkeit nachgewiesen 
und fiir die héheren Koeffizienten eine Naherungsauflésung gegeben wurde. 


(Hingegangen am 6. Oktober 1950.) 


Anschrift des Verfassers: Prof. Dr.-Ing. I. Szabé, Berlin-Charlottenburg 2, Technische Universitat, 
Hardenbergstr, 34. 
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Lésung einer Aufgabe tuber Stromung in Leitungssystemen *, 
Von W, Richter, 


1, Einleitung, Im Folgenden soll cine Aufgabe behandelt werden, die fiir Pumpenleitungen> 
Tunnelbeliiftung, Grubenbewetterung, Saugzug- und Unterwindanlagen in Kraftwerken, ver- 
schiedene industrielle Ventilationszwecke vorwiegend in der chemischen Industrie und fiir andere 
Gasleitungssysteme von Bedeutung ist. 

An der Eintrittsstelle E stréme das Férdermittel in den zu betrachtenden Teil eines Leitungs- 
systems ein und verlasse diesen an der Austrittsstelle A. Das Leitungssystem sei iiberall voll- 
kommen dicht. Die Unterschiede der geodatischen Héhen sollen vernachliassigt werden. Die 
Dichte des Férdermittels werde im betrachteten Teil des Leitungssystems als konstant voraus- 
gesetzt, und es gilt daher in einer iiberall dichten einfachen Leitung fiir den Férderstrom Q die 
Kontinuitatsgleichung 


| Q = konst [Q]= m?/s . 


Ein Leitungssystem soll passives System genannt werden, wenn in ihm dem Férdermittel 
keine Energie zugefithrt wird, in ihm also insbesondere kein Liifter bzw. keine Pumpe arbeitet. 
In einem passiven System wird dem Férdermittel itberall durch Wandreibung und an bestimmten 
Stellen des Leitungssystems durch Querschnittsinderungen, Umlenkungen und andere értliche 
Widerstande Energie entzogen. Es sei P;, = P;— P, die Gesamtdruckdifferenz zwischen den 
Stellen i und k gemessen in kg/m?. 

Fir stationare turbulente Strémungen in einer iiberall dichten passiven einfachen Leitung 
besteht dann die Beziehung 

Ri, = Piz - (1) 
Der darin enthaltene R-Wert ist von den Abmessungen, der Form und der Rauhigkeit der Lei- 
tung, von den in ihr liegenden ortlichen Widerstaénden, von der Dichte des Férdermittels und 
auBerdem noch vom Sinn, in welchem die Leitung durchflossen wird, abhangig. Bei der Um- 
kehrung des Sinnes werden ja bisherige Querschnittserweiterungen zu Querschnittsverengungen 
und umgekehrt, und auch manche anderen értlichen Widerstande ergeben dabei verschiedene 
Druckverluste. Im allgemeinen ist daher Rj, R,;. In jedem Falle haben R und P;;, gleiche 


Vorzeichen. 


; In zwei Sonderfallen kann man den R-Wert fiir ein bpe Nigel 7a aN Ln 
Leitungssystem einfach angeben, wenn man die R,- Werte OE alk aes a a 
aller seiner Einzelleitungen kennt. Abb. 1. Hintereinander geschaltete Leitungen. 
) 1. Fiir hintereinander geschaltete Leitungen (Abb. 1) 
gilt iy 
nm 
R= 5 K;. 
t=] 
| 2. Fir parallel geschaltete Leitungen (Abb. 2) ergibt — 9 ae 
sich 
n 

ly om 1 

VR oe VR, 5 , Abb. 2. Parallel geschaltete Leitungen. 

ea aN 


Im Folgenden soll fiir ein Leitungssystem, das sich nicht aus diesen beiden Grundformen 
aufbauen JaBt, eine graphische Lésung mit nomographischen Hilfsmitteln angegeben werden. 


2. Zwei parallel geschaltete Leitungen mit einer Querverbindung als passives System, Ein 
iiberall dichtes passives System habe die in Abb. 3 dargestellte Form. Im Punkte E teilt sich 
dic Zustrémleitung in zwei Teilleitungen, die sich im Punkte A wieder zur Abstrémleitung ver- 
einigen. Vom Punkte 1 der einen dieser Tcilleitungen fiihrt eine Querverbindung zum Punkte 2 
der anderen Teilleitung. Das so beschriebene Leitungssystem besteht also aus fiinf einfachen 


* Uber einen Teil der hier gegebenen Entwicklungen wurde auf der Tagung der Gesellschaft fiir ange- 
wandte Mathematik und Mechanik im April 1950 in Darmstadt berichtet. 
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Leitungen, fiir welche fiinf Bezichungen von der Form (1) aufgestellt werden kénnen. Die fiint 
entsprechenden R-Werte scien bekannt. Mit P soll die positiv vorausgesetzte Gesamtdruck; 
differenz P — Pry= Py — P, bezeichnet werden. Aus P > 0 folgt fiir ein iiberall dichteg 
passives System Pr, > 0, Py, [> 0, Pr, [0, Pos = 0. | 
Dadurch ist in beiden Teilleitungen des Systems der in Abb.3 eingezeichnete DurchfluBsing 
festgelegt, und die entsprechenden R-Werte fiir (1) sind ebenfalls eindeutig R,,,R, R,, .R. Ent 
sprechend Vie 0 sind in der Quer- 
verbindung | bis 2 zwei verschieden af 
Sinne fiir den Durchflu8B méglich und 
dazu gehéren i. a. zwei verschiedenex 
R-Werte Ry. baw. Ry. 
Die in E einstrémende Menge v4 

; teilt sich dort in Q, und Q,. Wird diey 
Ay é of in der Querverbindung strémendey 

Abb. 3. Passives System. Menge Q, als positiv bzw. negatiy 

angenommen, je nachdem ob sie vont 

I nach 2 oder von 2 nach | flieBt, so gehért in beiden Fallen zum Teilstrang 1A die Menges 
Q,— Q; und zum Teilstrang 2.A die Menge Q,-++- Q,. Man erhalt daher folgendes Gleichungssystem 


Oi VO =—-05 ems | 


R, 


R, OF = Pr, ° 
Re O = Pr ° 
Rk, Q: ae P,, ° 
1R (Q, — Qs)? = Pua, | (2 
2h (Qe. Qs)? 2 P,4 ? | 
Pr + 1A i > | 
Pra Py a1 
Pre eepeeceRs, = P, 


Dabei bezeichnet R, je nach der zugehérigen Strémungsrichtung einen der beiden Werte R,, 
bzw. R,,. 


Das Gleichungssystem (2) legt die neuen Gréfen Q, Q,, Qo, Q3, Per, Pro, Pra, Pos, Pi fir 
jeden Wert des Parameters P fest. Mit Ausnahme vonQ, und P,, sind diese GréBen nur dann Lé 
sungen der vorgelegten Aufgabe, wenn sie nicht negativ sind. Die Werte P,, und Q, miissen gleiche | 
Vorzeichen haben, und dieses Vorzeichen muB dem unter den beiden fiir R, méglichen Werten i 
gewahlten Wert entsprechen. AuBerdem gelten die B-dingungen |Q,;| <Q und |P,.| < P. Ist 
aber zu einem bestimmten Wert von P eine Gruppe der gesuchten neun GréBen so bestimmt, 
daf fiir sie das Gleichungssystem erfiillt ist, so gilt dasselbe fiir jede Gruppe, die man erhilt. | 
wenn man alle darin enthaltenen Férderstréme mit der beliebig gewahlten Konstanten k und!) 
alle darin enthaltenen Druckdifferenzen mit k? multipliziert. Daraus folgt fiir das betrachtete | 
Leitungssystem die Beziehung 


RG— P, (3) 


wo R den zu bestimmenden R-Wert des ganzen Leitungssystems bezeichnet. Aus den Glei- | 
chungen (2) ergibt sich | 


1R, OF + RQ = 2,R0, 0.4 Bs | 
oR. Q, + ok Q? = — 2,RQ,0,+ P, f (4) | 
RO} — RO} Rs. 02 —— 0, J | 


worin ,R, = R,-+- ,R der (1) entsprechende R-Wert fiir die ganze erste Teilleitung E1A und 
2R,= R, + 2R jener fiir die ganze zweite Teilleitung LH 2 A bedeutet. SchlieSlich erhalt man aus 
(4) das Gleichungssystem 


Qi = a, Q, Q3 + 8, Q.Q3-+ ¢ P, (5a) | 
Q3 = 42.Q; Q3 + by Qe Qs + co P, (5b) 
G2 = a,Q, 03+ bs Q2 Qs ++ e, P, (5¢) 
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wenn die Bezeichnungen 
_ 2 (12Ro3 + 2k.) 


a, = ae : b, = =< ge SY oR) ; 
eS A = aR 
a, = N : b) ay ve ? 
maa a (i = 1, 2, 3) 
9 eh 2R 9 9 
und 
aft OTR 
N= 0 eh, oh 
Ae R, Ral 


eingefiihrt werden, N+ 0 ist und die Symbole 
| oki = otk (.R-+ R,) | Senge | 
ik Rim = iR raat R, Re | 


beniitzt werden. Die Gleichungen des Systems (4) sind aber fiir N= 0 nur dann gleichzeitig 
fiir nicht negative Werte von Q,, Q,,,R und ,R méglich, wenn Q, = Q, Q; P= (ist, dh. 
wenn das Leitungssystem nicht durchstrémt wird. 


(i,k. L m= I, 2, 3) 


3. Graphische Lésung. Zur Lésung des Gleichungssystems (5) sollen nomographische Hilfs- 
mittel verwendet werden. Die Gleichung (5a) ergibt ein Fluchtliniennomogramm mit drei 
Leitern, die in rechtwinkligen Koordinaten durch folgende Parameterdarstellungen bestimmt 
sind: 


F b Q? —c,P 
-Leiter: x= — —+, eee SU 
@ a, Q; 2 a, Q, 
| Q,-Leiter : eee ia UR 
| 0, 
| Oretter: «= 0, Vers 


Die Q,-Leiter ist eine Hyperbel mit den Asymptoten «= 0 (y-Achse) und y= (c, P/b,) x, 
deren Mittelpunkt der Koordinatenursprung ist. Ebenso folgt fiir (5b) ein Fluchtliniennomo- 
gramm mit drei Leitern: 


, Uh 
Q,-Leiter: += 0,” y= 9; 
— Q2 —c,P 
Reg ° rt = a, ; ze SON a Se 
Q,-Leiter: x oro: ¥ bao: 
Q,-Leiter: x= 0, va», 


Hier ist die Q,-Leiter eine Hyperbel mit den Asymptoten x— 0 (y-Achse) und y = (ce, P/a,) x, 
deren Mittelpunkt der Koordinatenursprung ist. 
Endlich erhalt man fiir (5c) ein Fluchtliniennomogramm mit drei Leitern: 


. b. I 
Q;-Leiter: 4 = a0; > a= Q,’ 
Q,-Leiter: «= a > pe 

a, Q. 
Q,-Leiter: x= 0, si Gi =aP ; 


Diese drei Leiter sind Gerade durch den Koordinatenursprung, von denen zwei mit den Ko- 


ordinatenachsen zusammenfallen. 

Es empfiehlt sich, wie in Abb.4 alle drei Fluchtliniennomogramme zu vereinigen, weil sie 
eine Leiter gemeinsam haben. 

Nun kann das Gleichungssystem (5) graphisch gelést werden. Wird fiir Q, ein beliebiger 
fester Wert Q' gewahlt, so liefern die zu den Gleichungen (5 a) und (5c) gehérigen Teile des 
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Nomogramms Abb. 4 fiir jeden gewahlten Wert Q,— Q* zwei im allgemeinen verschiedene 
Werte Q" bzw. Q*’ fiir Q,. Wird fiir Q, eine Reihe verschiedener Werte Qk gewahlt, so kann man 
die in Abb. 5 angedeuteten Kurven [Q*’] bzw. [Q/’], namlich eine Kurve zweiter Ordnung und 
eine Gerade, zeichnen. Zu ihrem Schnittpunkt P, fiir den QY” = QY’ = Qi ist, gehért mit 
Q; = Qi der Punkt P’. Wird dieser Vorgang fiir eine Reihe von Werten Q,= Q; durchgefihrt, 


A 
1 
7 & 
le A Ee: A 


(0) 
(2) 


Abb. 4. Nomogramm fiir das passive System Abb. 6. Sonderfille: einer der R-Werte gleich Null, 
4, g : 


Q Q 


@ @) Hy 


Abb. 5. Teilforderstréme fiir das passive System. Abb. 7. Sonderfalle: einer der R-Werte gleich Null. 


so erhalt man die in Abb. 5 angedeuteten Kurven [P] bzw. [P’]. Die zu den Punkten der Kurve 
[P] gehérigen Wertepaare Q”, Q™ liefern in dem zu (5b) gehérigen Teil des Nomogramms 
Abb. 4 jeweils einen Wert Q™ fiir Q,, der im allgemeinen nicht gleich Q: ist. Der zugehérige 
Punkt P” ist in Abb. 5 angedeutet. Die Punkte P’ liegen auf einer Kurve [P”], deren Schnitt- 
punkte x mit der Kurve [P’] als Koordinaten Lésungen des Gleichungssystems (5) haben, 
unter denen sich die fiir die gestellte Aufgabe geeigneten befinden. Fiir ein solches Lésungs- 
system ist R— P/Q? der*R-Wert des betrachteten Systems. 
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| Die Abb. 6 und 7 stellen Sonderfalle des betrachteten Systems dar, in denen die R-Werte 


Beevelncr Teilleitungen Null werden. 


4, Beispiel fiir ein passives Lei- 
|} tungssystem, Fiir ein Beispiel wur- 


den folgende Annahmen getroffen: 


R,= 0,0033, R= 0,0066, 
R,= 0,0051, R= 0,0053, 
R, = 00,188. 


| Abb.8 zeigt das fiir P = 100kg/m? 
daraus abgeleitete Nomogramm. 
Die Lage der Koordinatenachsen 
, und die auf ihnen gewahlten Zeichen- 
einheiten zeigt die Tabelle. 


Abb. 8. Nomogramm fiir das passive System: 
Anwendungsbeispiel 


Aus der in Abb. 9 dargestellten 
graphischen Lisung, die mit Hilfe 
j}des Nomogramms der Abb.8 in der 
jin Abschnitt 3 angegebenen Weise 
) bestimmt wurde, folgt 


Q, = 108 m?/s, 0;— 90,5 m/s, 
QO; == 13,1 m*/s 


und daher der gesamte Férderstrom 
des betrachteten Systems 


0= 0,+ Q,= 198,5 m3/s , 
ifiir welchen sich zu P = 100 kg/m? 
jder R-Wert des passiven Systems 
) R= 0,00254 ergibt. Weil Q, positiv 
hist, strémt das stromende Mittel in 
)der Querverbindung vom Punkt ] 
»zum Punkt 2. In Abb. 10 wurden fiir 
}das ganze passive System und fiir 
‘seine fiinf Teilleitungen die Wider- 


o9 
10 ¢; @ 
YAR 
109 i) 
igs v7 me ¥ 
~008 a 
el Acegonss 7 te = 00053 
106 Se e+ Ks 
SS 
105 XN 
9 NN 
4 \ pce 
Ss 
MI-\\|-93 Soo 
SS ~ 
102 
l Se 
107 fs is 
100 Pat, 
Los SS 
9 | Ss 
Na “XN 
98 1 zm 
0 600 500 400 300 200 700 Wb — IO 
- 
oo 
—s 
5 “3 ae 
pee 
TE 
=a 
ee 
mm 
oat 
an 
oma 
Nomo- Abszissen + x Ordinaten + y 
ere Richtung Richtung 
Glevehuns horizontal Einheit vertikal Einheit 
eee nach nach 
5a rechts oben 2 cm 
b Sb 40 m 
DEC links unten 40 m 
Qs; 
20 


Lue 


! 


200 Ws 


700 


-100 


Ye, 


105 


ms 70 


Abb. 9, Teilférderstréme fiir das passive System: Anwendungsbeispiel. 
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standsparabeln und auf ihnen die zu den ermittelten Lésungen gchiérenden Betriebspunkte | 
dargestellt. 


5. Aktives System. Ein Leitungssystem soll aktives System genannt werden, wenn in ihm | 
dem strémenden Mittel Energie zugefiithrt wird. Das kann z. B. so geschehen, dafi in ihm Liifter 


bzw. Pumpen arbeiten. Bezeichnet ,P4 den Gesamtdruck unmittelbar vor und ,P, den 


AP 
100 iT 
kg/m f OL | karl fe ee 
80 | a | 


60 
iW 
| | 
| | 

( 
| 

40 + fe + 

VAR 
| 
i} 
| 


20 ~ y a 


Be H ! Q 
0 700 m/s 200 


Abb. 10. Widerstandsparabeln fiir das passive System und seine Teilleitungen. 


Gesamtdruck unmittelbar hinter einem Lifter bzw. einer Pumpe, so stellt ;Pe4= 1Pr— 1P4 
die erzielte Gesamtdruckerhéhung dar und Q:;P 4 ist ein Maf fiir die dabei dem strémenden 
Mittel zugefiihrte Leistung. Zwischen Q und ,P;4 besteht ein dem betreffenden Liifter (bzw. 
ri der Pumpe) ecigentiimlicher Zu- 
Y 

a ; sammenhang, dessen Darstellung 
gj? Le & ly | in einem rechtwinkligen Koordi- 
natensystem mit Q als Abszissen 
und ,;Pr, als Ordinaten die 
Kennlinie des Liifters bzw. der 
Pumpe heiBt. Die Kurve K in 
Abb. 11 ist eine solche Kennlinie 
fir einen Axiallifter Bauart 
Schicht (D = 3,40m, n= 440 
U/min, y= 1,25 kg/m?). Arbeitet 
dieser Liifter in einer einfachen 
Leitung, deren R-Wert gegeben 
ist, und in der dem strémenden 
Mittel auBer durch diesen Lifter 
keine Energie zugefiihrt wird, so 
erhalt man fiir diese, nachdem 
man sich den Liifter aus der 
Abb. 11. Aktives System: Liifterkennlinie HOPES ae ee el 

oe arid i Puz = RQ?, wofiir man die 

Kurve W ins Koordinatensystem der Liifterkennlinie in Abb, 11 einzeichnet. Wird zu jedem Q 
als Abszisse sPe4 = 1 Pra— Par als Ordinate eingezeichnet, so erhalt man die Kurve S. Ihre 
Punkte bestimmen fiir das betreffende aktive System, das aus der betrachteten einfachen Lei- 
tung und dem in ihr arbeitenden Liifter besteht, zu jedem Q die zugehérige Gesamtdruckdifferenz 


100-200 
| 


: | 100 
at 40 
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sPre4= fs (Q) als Funktion von Q. Die Beziehung 
Rie sPra 


Q? 

ordnet jedem Wertepaar Q, 5Pr4 eindeutig einen R-Wert zu. Zu jedem Punkt der Kurve S 
gehért daher cin ganz bestimmter R-Wert. In Abb. 11 ist die Kurve R eingezeichnet, die aus 
den Punkten der Kurve S so abgeleitet wurde. Ein aktives System unterscheidet sich daher 
von einem passiven dadurch, daB sein R-Wert nicht konstant ist, sondern vom Q abhangt. 

Wird ein solches aktives System mit einem passiven System vereinigt, wie dasin Abb. 12 gezeich- 
net ist, so daS der Endpunkt ,E des aktiven Systems mit dem Anfangspunkt ,A des passiven 
Systems und der Anfangspunkt ,A des aktiven Systems mit dem Endpunkt ,E des passiven 


Systems zusammenfallen, so erhalt man ein geschlossenes aktives System. In Abb.11 ist zur 
Liifterkennlinie K und zu den fiir 


den aktiven Teil des Systems be- gr a) g; 
stimmten Kurven S und R noch die ote 
die zum passiven Teil des Systems 

gehérige Kurve G gezeichnet. Der fig 
Schnittpunkt X der Kurven S und 

G bestimmt den Forderstrom Q des (Of 
Systems und ,P4n x = «Pri_x die 
Gesamtdruckdifferenz zwischen den 
Vereinigungspunkten des aktiven 


—e 
(2) 
@) (4) 
Abb, 12. Geschlossenes aktives System. Abb, 13. Nomogramm fiir das aktive System: Liifter in der Querverbindung 1—2. 


und des passiven Systems. Zum Punkt X gehért fiir dasselbe Q auf der Kurve K der Arbeits- 
punkt Y des Lifters, der die vom Liifter erzeugte Gesamtdruckdifferenz ;P,4 festlegt. Der 
Abb. 11 liegt als passives System jenes zugrunde, das in Abschnitt 4 behandelt wurde. Fir 
die zum aktiven System gehérende Kurve W wurde R= 0,0035 gewahlt und der dort ar- 
beitende Liifter war so ausgelegt worden, da zu seinem Arbeitspunkt Y der Bestwert des 
Liifterwirkungsgrades 7 = 85°, gehért. 


6. Zwei parallel geschaltete Leitungen mit einer Querverbindung als aktives System. Arbeitet 
in einer der Teilleitungen des in Abb.3 dargestellten Leitungssystems ein Liifter bzw. eine Pumpe, 
so liegt ein aktives System vor. Es soll zuerst der Fall besprochen werden, daf} in der Querver- 
bindung 1—2 ein Liifter arbeitet und alle anderen Teilleitungen des Leitungssystems passiv sind. 
P fiir das Leitungssystem sei konstant und gegeben. Je nach der Art des Liifters ist jetzt 


Pie, a es ’ Pin 0) P,42 


2A = ? 


und der Sinn des Durchflusses in den Teilleitungen ist daher nicht mehr von vornherein bekannt. 
Wenn nun fiir Q, und Q, auch negative Werte und Null zugelassen werden, wobei jeweils posi- 
tiven Werten der in Abschnitt 2 angegebene Durchflu®sinn entsprechen soll, so gelten alle 
Gleichungen des Systems (2), nur ist hier R, nicht mehr konstant, sondern eine Funktion von Q;, 
die wie in Abb. 8 aus der Liifterkennlinie und dem R-Wert der Teilleitung 1—2 ableitbar ist. 
Die Determinante N ist eine lineare Funktion von R, und die a;, 6;, ¢; (t==1, 2, 3) sind linear 
gebrochene Funktionen von R;. In dem zu (5a) entwickelten ersten Teil des Nomogramms 
(Abb. 4) erhalt man anstatt der Q,-Leiter, und in seinem zu (5 b) gehérenden zweiten Teil anstatt 
der Q,-Leiter nun nicht mehr Kurven, sondern Kurvenscharen. 

Im ersten Nomogrammteil (Abb. 13) ergibt jedes Q, zu dem ihm entsprechenden R,-Wert 
cine mit den Werten von Q; zu beziffernde Hyperbel. Fiir alle Hyperbeln dieser Schar ist der 
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Koordinatenursprung Mittelpunkt und die Y-Achse eine Asymptote. Werden alle zu gleichen 
Q,-Werten gehérenden Punkte verbunden, so erhalt man eine zweite mit Q, zu beziffernde Kur- 
venschar. Analog findet man fiir den zweiten Nomogrammteil ein Kurvenfeld, bestehend aus 
einer mit Q, bezifferten Kurvenschar und aus einer mit Q; bezifferten Hyperbelschar. Der dritte 
Nomogrammteil unterscheidet sich von dem der Abb.4 nur durch die andere Teilung der Q;- 
Leiter. Die Benutzung des so entwickelten 
(@) @) Nomogramms geschieht ganz dahnlich, wie 

of das in Abschnitt 3 fiir ein passives System 
beschricben wurde. 
Wenn der Liifter in der Teilleitung E 1 _ 
arbeitet, sind N und die a;,, b;, ¢; (t= 1, 2, 3) 
Funktionen des nunmehr von Q, abhangigen 
R,-Wertes. Im ersten Nomogrammteil erhalt 
man anstatt der mit Q, bezifferten urspriing- 
lichen Hyperbel eine andere Kurve. Im zwei- 
ten Nomogrammteil ergibt sich ein mit Q, 
bzw. Q, beziffertes der Abb. 14 ihnelndes 
Kurvenfeld, und im dritten Nomogrammteil 
ist die Q,-Leiter eine Kurve und nicht mehr 
einc Gerade, wahrend die Q,-Leiter zu einer 
Schar von je zu einem Q,-Wert gehérenden 
auf der Y-Achse liegenden Leitern wird. 
Dieser Nachtcil wird durch die Transformation 


f=x+taQ,y, («4 0, willkirliche | 
(ay n= y Konstante) 


Abb. 14. Nomogramm fiir das aktive System: Lifter in der beseitigt wodurch die QO - Leiter in eine andere 
Teilleitung El. oe % : 1 : < : : 
Form ibergefiihrt wird und die Q,-Leitern in 
ein mit Q, beziffertes Geradenbiischel mit dem Koordinatenursprung als Scheitel und eine mit Q, 
bezifferte Kurvenschar itbergehen. Alle anderen Aufstellungsméglichkeiten des Liifters konnen 
Joe (V2 auf die eben besprochene 
gu 200 zurickgefiihrt werden. 


7. Beispiel fiir ein aktives 
Leitungssystem, In der Teil- 
leitung 1— 2 soll ein Liifter 
eee ve| L’ arbeiten, dessen Kenn- 
linie die Kurve K in Abb. 15 
ist. DL’ ist ein Axialliifter 
Bauart Schicht mit einem 
Laufradau8endurchmesser 
von 1400 mm und einer Dreh- 
zahl von 720 U/min. Unter 
Beriicksichtigung derin Ab- 
ae schnitt 5 fiir diese Teillei- 
. tung angegebenen Konstan- 
ten R, erhalt man daraus 
% die KurveS, fiir deren Punk- 
g, te die jeweiligen R-Werte 

bs 0 25 ey m/s ow berechnet und in Kurve R 

Abb. 15. Kennlinie des Lifters in der Querverbindung. in Abhiangigkeit vom Foér- 

derstrom aufgetragen wur- 

den. Wenn der Lifter L’ in der Richtung vom Punkt 1 nach Punkt 2 férdert, so sind diese 

R-Werte im dargestellten Bereich des Férderstroms negativ bezogen auf P,.. Dem Nomogramm 

in Abb. 16 liegen dieselben Annahmen iber die Koordinatenachsen und die MaBstabcinheiten 
zugrunde wie jenem in Abb. 8. Fiir P= 100 kg/m? ergab die graphische Lisung 

C= 120 mers 03-80%? /s 5 Q; == 31,5 m3’s 


und daher als Férderstrom fiir das gesamte betrachtete System 


Q= 0+ Q,= 200 mi/s , 
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also einen vom Beispiel in Abschnitt 4 wenig verschiedenen Wert. Wird daher das betrachtete 
System, wie das in Abb. 17 dargestellt ist, wieder mit demselben aktiven System wie in Ab- 
-schnitt 4 zu einem Kreislauf vereinigt, so erhilt man fiir dessen Liifter L denselben Arbeitspunkt. 
Fir den Lifter L’ findet man den in Abb. 15 angegebenen Arbeitspunkt Y,,. Dadurch, dah 

diesmal in der Teilleitung 12 der Liifter L’ arbeitete, wurde fiir den gesamten Kreislauf nichts 
_ geandert, in den fiinf Teilleitungen des betrachteten Systems wurde jedoch eine andere Verteilung 
der Férderstréme und der Gesamtdruckdifferenzen als im Abschnitt 4 erzielt. 

Nimmt man jedoch an, daB derselbe bs 
| ‘Hilfsliifter L so in die Querverbindung 
-eingebaut wurde, da®B er im umgekehr- 
ten Sinn, also von 2 nach 1 blast, wah- 
-rend der R-Wert der passiven Quer- 
_verbindung 1—2, also ohne eingebauten 
Lifter, fiir beide Durchstrémrichtungen 
gleich vorausgesetzt werden soll, so muB 
man die R-Werte in Abb.15 mit —1 
multiplizieren; sie sind dann positiv. 
Das der Abb. 16 fiir diesen Fall ent- 
sprechende Nomogramm ergibt zu 


P= 100 kg/m? 


die Lisung 

Q, = 80 m/s, 
OF== 111 dims , 
Q,= 28,15 m/s . 


Der Gesamtforderstrom 


O02 O1-- Oo — 1915 m/s 


unterscheidet sich dann von dem in den OM oan ae wie a tr 
friiheren Beispiclen erhaltenen. Damit g 500280 tO 180 3000 ae a 
verschiebt sich aber der Arbeitspunkt Y 2 ood 

in Abb.1]1. Das angegebene graphische oe 

Verfahren mu® fiir verschiedene Werte ve eee 

von P wiederholt werden, wobei sich pet 

im Nomogramm eee 

Abb.16 in beiden ae 

Kurvenfeldern je- oe 2 ‘ 


weils die Ordinaten 
‘der Netzpunkte 
im. wahrend 
die Abszissen un- 


Abb. 16. Nomogramm fiir das aktive System: Anwendungsbeispiel. 


_ verandert bleiben. 

‘Von den Einzelleitern andert sich nur die Q.-Leiter des dritten Nomogrammteiles. Die so er- 
‘mittelten Wertepaare P,Q kénnen als Kennlinie des betrachteten aktiven Systems in Abb. 11 
‘eingezeichnet werden, und es kann so der neue Arbeitspunkt des Liifters L bestimmt werden. 
Im vorliegenden Beispiel erhalt man als Lésung 


0; == 82;5'm*/s , O,= 114, m?/s, Q, = —28.7 m,/s , 


find zwar fiir die zum Gesamtforderstrom 
Q= Q, + Q,= 196,5 m/s 


gehérige Gesamtdruckdifferenz P= 105 kg/m”. . 

Fiir den Hilfsliifter L’ wurde der diesen Werten entsprechende Arbeitspunkt Y,, in Abb. 15 
cingezeichnet. pat 

In Abb.18 sind die Férderstréme in den Teilleitungen und die ihnen entsprechenden Gesamt- 
druckdifferenzen fiir die drei angegebenen Beispiele miteinander verglichen. Man erkennt, da® 
bei diesen Beispielen der Hilfsliifter L’ im wesentlichen die Verteilung der Forderstréme auf 
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die cinzelnen Teilleitungen beeinfluBt, gerade das ist aber oft bei den praktischen Anwendungen 
die Aufgabe eines solchen zusitzlichen Liifters L’, 

Es mu stets beachtet werden, ob der Liifter L’ irgendwo im System innere Kreislaufe ver- 
ursacht, oder ob durch seinen Einsatz in einzelnen Teilleitungen iiberhaupt keine Strémung mehr 
stattfindet. In den meisten Anwen- 
dungsfallen sind innere Kreislaufe und 


des Systems unerwiinscht, oft sogar 


des Liifters L durch Einsatz des Hilfs- 
lifters L’ stark beeinfluBt wird, so ist es 


neue Liifter L’ den urspriinglichen Liif- 


Lifter L hindurch seiner urspriinglichen 


Arbeitsrichtung entgegengesetzt ver- 
Abb, 17. Geschlossenes aktives System mit 2 Liiftern. laiuft. 


8. Zusammenfassung. Fir ein iiberall dichtes Leitungssystem, das aus zwei parallel geschal- 


teten Leitungen mit einer Querverbindung besteht, wurden mit nomographischen Hilfsmitteln | 


zuerst fiir den Fall, daB in ihm der Strémung keine Energie zugefithrt wird, und dann fiir den 


Qg 
i g CO) ohne Hilfslitter in der Querverbindung 12 
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Abb. 18. Vergleich der Foérderstréme und der Druckdifferenzen fiir die besprochenen Beispiele. 


Fall, daB in einer seiner Teilleitungen cin Liifter bzw. eine Pumpe mit bekannter Kennlinie 
arbeitet, die Verteilung der Férderstréme und die zugehérigen Gesamtdruckdifferenzen zwischen 
je zwei Verzweigungspunkten bestimmt. Es wurden Sonderfialle besprochen und Beispiele fiir 
das Verfahren durchgefiihrt. Die Verwendung nomographischer Hilfsmittel liegt nahe, weil die 
Liifter- bzw. Pumpenkennlinie nicht in analytischer Form festgelegt werden kann und die Glei- 
chungen des Systems (5) quadratisch und nicht homogen sind. 


Stockungen der Strémung in Teilen. 


sehr gefahrlich. Wenn der Arbeitspunkt | 


im Extremfall sogar méglich, daB der | 


ter L iiberspiilt, so daB dieser, anstatt | 
dem Kreislauf Energie zuzufithren, nur | 
ein Strémungshindernis darstellt, oder | 
daB sogar die Strémung durch den | 
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Es soll darauf besonders hingewiesen werden, daB jedoch keinesfalls durch eine solche Lisungs- 
midglichkeit Erfahrungen auf den Gebieten, in die das betreffende Leitungssystem gehért, und 
die Erfahrungen des Liifter- bzw. Pumpenbaues entbehrlich werden. Im Gegenteil kann nur aus 
diesen Quellen eine verniinftige Wahl der erforderlichen Konstanten und der einzusetzenden 
Lifter bzw. Pumpen stammen. Die praktische Anwendbarkeit ist zudem so schwierig und eine 
Fehllésung verursacht stets, vor allem auf dem Gebiete der Grubenbewetterung, Betriebszu- 
stande, die nur der véllig mit dem betreffenden Sondergebiet Vertraute iiberblicken und verant- 
worten kann, so dafB die evtl. praktische Anwendung des entwickelten Verfahrens ausschlieBlich 
fiir diese Fachleute bestimmt sein kann. 

Eine Analogie mit Fragen auf dem Gebiete der Elektrotechnik liegt nahe, wo iibrigens auch 
eine experimentelle Behandlung der Aufgabe verhaltnismabig einfach méglich ist. Das ist auch 
bekannt und ein Vergleich der experimentell gewonnenen Erkenntnisse mit dem entwickelten 
Verfahren ware wiinschenswert. Im Falle eines aktiven Systems diirfte die experimentelle Dar- 
stellung des elektrischen Analogons fiir den in einer der Teilleitungen arbeitenden Liifter mit 
als bekannt vorausgesetzter Kennlinie einige Schwierigkeiten bereiten. 


(Eingegangen am 9, Oktober 1950.) 
Anschrift des Verfassers: Dr. W. Richter, Frankenthal, Beindersheimer Strafe 13. 
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Berichtigung 


zu meinem Aufsatz in Band XVIII, S. 394 des Ingenieur-Archivs ,,Beitrage zur 
Theorie der durch gleichférmigen Schub beanspruchten Platte. II. Mitteilung™. 


Von W. Kucharski. 


In der genannten Arbeit habe ich unter (24) auf $.397 irrtiimlicherweise angegeben, da8 die 
Bewegung stabil, und zwar periodisch verlauft, wenn die aus der Fundamentalgleichung (23) 
der maBgebenden linearen Substitutionen folgende Konstante s = + | ausfallt. Dies trifft nur 
dann zu, wenn s nicht gleichzeitig eine Doppelwurzel der Fundamentalgleichung ist. Im Falle 
der Doppelwurzel artet die Substitution in bekannter Weise aus, so da dann zu der periodi- 
schen Liésung eine nichtperiodische hinzutritt, die mit der Zeit anwachst. Hier liegt gerade 
dieser Fal] vor; in (24) ist die auf s= + 1 beziigliche Aussage im angegebenen Sinne zu korri- 
gieren, und in dem Stabilitatskriterium (28) auf S.398 ist das Gleichheitszeichen zu streichen. 
Die den spater diskutierten und in Abb.2 skizzierten Grenzkurven entsprechenden Bewegungen 
sind daher im allgemeinen nicht mehr stabil. 


Wenn auch die Lésungen der hier behandelten partiellen Differentialgleichung naturgemaf | 


keine Matthieuschen Funktionen sind, wie sie aus der gewohnlichen Differentialgleichung (18) 
in bekannter Weise entstehen, sondern sich von diesen in vieler Beziehung grundsatzlich unter- 
scheiden (man beachte z. B. die starken Abweichungen von Abb.2 gegeniiber der bekannten 
sog. Struttschen ,,Karte‘‘), so haben sie doch diejenigen Eigenschaften gemeinsam, die aus dem 
Auftreten der gleichartigen linearen Substitutionen folgen. Zu diesen gehéren die Besonder- 
heiten, welche im Falle einer Doppelwurzel fiir die Fundamentalgleichung auftreten. 


(Eingegangen am 10. Mai 1951.) 
Anschrift des Verfassers: Prof. W. Kucharski, Berlin-Charlottenburg, Hardenbergstrabe 34. 
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